Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.

Well-Aligned Zn-Doped InN Nanorods Grown by Metal-Organic Chemical Vapor Deposition and the Dopant Distribution

We report the synthesis and characterization of Zn-doped InN nanorods by metal-organic chemical vapor deposition. Electron microscopy images show that the InN nanorods are single-crystalline structures and vertically well-aligned. Energy-dispersive X-ray spectroscopy analyses suggest that Zn ions are distributed nonhomogenously in InN nanorods. Simulations based on diffusion model show that the doping concentration along the radial direction of InN nanorod is bowl-like from the exterior to the interior, the doping concentration decreases, and Such dopant distribution result in a bimodal EDXS spectrum of Zn across the nanorod. The study of the mechanism of doping effect is useful for the design of InN-based nanometer devices. Also, high-quality Zn-doped InN nanorods will be very attractive as building blocks for nano-optoelectronic devices.'

Stone-Wales defects created by low energy recoils in single-walled silicon carbide nanotubes

The defect creation at low energy events was studied using density functional theory molecular dynamics simulations in silicon carbide nanotubes, and the displacement threshold energies determined exhibit a dependence on sizes, which decrease with decreasing diameter of the nanotubes. The Stone-Wales (SW) defect, which is a common defect configurations induced through irradiation in nanotubes, has also been investigated, and the formation energies of the SW defects increase with increasing diameter of the nanotubes. The mean threshold energies were found to be 23 and 18 eV for Si and C in armchair (5,5) nanotubes. (C) 2009 American Institute of Physics. [doi: 10.1063/1.3238307]

Experimental study of the high-voltage breakdown and simulations for the RFQ cooler and buncher RFQ1L

The RFQ cooler and buncher RFQ1L is one of the key parts of the being-built super-heavy nuclide research spectrometer. In order to understand the high-voltage breakdown phenomenon, the voltages between electrodes have been measured. In addition, more extensive simulations have been performed for better understanding and optimizing the RFQ1L work points.

Isobaric yield ratios and the symmetry energy in heavy-ion reactions near the Fermi energy

The relative isobaric yields of fragments produced in a series of heavy-ion-induced multifragmentation reactions have been analyzed in the framework of a modified Fisher model, primarily to determine the ratio of the symmetry energy coefficient to the temperature, a(sym)/T, as a function of fragment mass A. The extracted values increase from 5 to similar to 16 as A increases from 9 to 37. These values have been compared to the results of calculations using the antisymmetrized molecular dynamics (AMD) model together with the statistical decay code GEMINI. The calculated ratios are in good agreement with those extracted from the experiment. In contrast, the values extracted from the ratios of the primary isobars from the AMD model calculation are similar to 4 to 5 and show little variation with A. This observation indicates that the value of the symmetry energy coefficient derived from final fragment observables may be significantly different than the actual value at the time of fragment formation. The experimentally observed pairing effect is also studied within the same simulations. The Coulomb coefficient is also discussed.

Fluorene-Based Copolymers Containing Dinaphtho-s-indacene as New Building Blocks for High-Efficiency and Color-Stable Blue LEDs

By incorporating a new building block, 7,7,15,15-tetraoctyldinaphtho-s-indacene (NSI), into the backbone of poly(9,9-dioctylfluorene) (PFO), a novel series of blue light-emitting copolymers (PFO-NSI) have been developed. The insertion of the NSI unit into the PFO backbone leads to the increase of local effective conjugation length, to form low-energy fluorene-NSI-fluorene (FNF) segments that serve as exciton trapping sites, to which the energy transfers from the high-energy PFO segments. This causes these copolymers to show red-shifted emissions compared with PFO, with a high efficiency and good color stability and purity. The best device performance with a luminance efficiency of 3.43 cd . A(-1), a maximum brightness of 6 539 cd . m(-2) and CIE coordinates of (0.152, 0.164) was achieved.

Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

Diffusion and single molecule dynamics on biomolecular interface binding energy landscape

We propose a new approach to study the diffusion dynamics on biomolecular interface binding energy landscape. The resulting mean first passage time (MFPT) has 'U'curve dependence on the temperature. It is shown that the large specificity ratio of gap to roughness of the underlying binding energy landscape not only guarantees the thermodynamic stability and the specificity [P.A. Rejto, G.M. Verkhivker, in: Proc. Natl. Acad. Sci. 93 (1996) 8945; C.J. Tsai, S. Kumar, B. Ma, R. Nussinov, Protein Sci. 8 (1999) 1181; G.A. Papoian, P.G. Wolynes, Biopolymers 68 (2003) 333; J. Wang, G.M. Verkhivker, Phys. Rev. Lett. 90 (2003) 198101] but also the kinetic accessibility. The complex kinetics and the associated fluctuations reflecting the structures of the binding energy landscape emerge upon temperature changes. The theory suggests a way of connecting the models/simulations with single molecule experiments by analysing the kinetic trajectories.

Temperature dependence of the distribution of the first passage time: Results from discontinuous molecular dynamics simulations of an all-atom model of the second beta-hairpin fragment of protein G

More than 22 000 folding kinetic simulations were performed to study the temperature dependence of the distribution of first passage time (FPT) for the folding of an all-atom Go-like model of the second beta-hairpin fragment of protein G. We find that the mean FPT (MFPT) for folding has a U (or V)-shaped dependence on the temperature with a minimum at a characteristic optimal folding temperature T-opt*. The optimal folding temperature T-opt* is located between the thermodynamic folding transition temperature and the solidification temperature based on the Lindemann criterion for the solid. Both the T-opt* and the MFPT decrease when the energy bias gap against nonnative contacts increases. The high-order moments are nearly constant when the temperature is higher than T-opt* and start to diverge when the temperature is lower than T-opt*. The distribution of FPT is close to a log-normal-like distribution at T* greater than or equal to T-opt*. At even lower temperatures, the distribution starts to develop long power-law-like tails, indicating the non-self-averaging intermittent behavior of the folding dynamics. It is demonstrated that the distribution of FPT can also be calculated reliably from the derivative of the fraction not folded (or fraction folded), a measurable quantity by routine ensemble-averaged experimental techniques at dilute protein concentrations.

Building Responsive Systems from Physically-correct Specifications

Predictability - the ability to foretell that an implementation will not violate a set of specified reliability and timeliness requirements - is a crucial, highly desirable property of responsive embedded systems. This paper overviews a development methodology for responsive systems, which enhances predictability by eliminating potential hazards resulting from physically-unsound specifications. The backbone of our methodology is the Time-constrained Reactive Automaton (TRA) formalism, which adopts a fundamental notion of space and time that restricts expressiveness in a way that allows the specification of only reactive, spontaneous, and causal computation. Using the TRA model, unrealistic systems - possessing properties such as clairvoyance, caprice, in finite capacity, or perfect timing - cannot even be specified. We argue that this "ounce of prevention" at the specification level is likely to spare a lot of time and energy in the development cycle of responsive systems - not to mention the elimination of potential hazards that would have gone, otherwise, unnoticed. The TRA model is presented to system developers through the CLEOPATRA programming language. CLEOPATRA features a C-like imperative syntax for the description of computation, which makes it easier to incorporate in applications already using C. It is event-driven, and thus appropriate for embedded process control applications. It is object-oriented and compositional, thus advocating modularity and reusability. CLEOPATRA is semantically sound; its objects can be transformed, mechanically and unambiguously, into formal TRA automata for verification purposes, which can be pursued using model-checking or theorem proving techniques. Since 1989, an ancestor of CLEOPATRA has been in use as a specification and simulation language for embedded time-critical robotic processes.

An Energy-conscious Transport Protocol for Multi-hop Wireless Networks

We present a transport protocol whose goal is to reduce power consumption without compromising delivery requirements of applications. To meet its goal of energy efficiency, our transport protocol (1) contains mechanisms to balance end-to-end vs. local retransmissions; (2) minimizes acknowledgment traffic using receiver regulated rate-based flow control combined with selected acknowledgements and in-network caching of packets; and (3) aggressively seeks to avoid any congestion-based packet loss. Within a recently developed ultra low-power multi-hop wireless network system, extensive simulations and experimental results demonstrate that our transport protocol meets its goal of preserving the energy efficiency of the underlying network.

JTP: An Energy-conscious Transport Protocol for Wireless Ad Hoc Networks

Within a recently developed low-power ad hoc network system, we present a transport protocol (JTP) whose goal is to reduce power consumption without trading off delivery requirements of applications. JTP has the following features: it is lightweight whereby end-nodes control in-network actions by encoding delivery requirements in packet headers; JTP enables applications to specify a range of reliability requirements, thus allocating the right energy budget to packets; JTP minimizes feedback control traffic from the destination by varying its frequency based on delivery requirements and stability of the network; JTP minimizes energy consumption by implementing in-network caching and increasing the chances that data retransmission requests from destinations "hit" these caches, thus avoiding costly source retransmissions; and JTP fairly allocates bandwidth among flows by backing off the sending rate of a source to account for in-network retransmissions on its behalf. Analysis and extensive simulations demonstrate the energy gains of JTP over one-size-fits-all transport protocols.

Electrical machine characterisation and analysis for renewable energy applications

There has been an increased use of the Doubly-Fed Induction Machine (DFIM) in ac drive applications in recent times, particularly in the field of renewable energy systems and other high power variable-speed drives. The DFIM is widely regarded as the optimal generation system for both onshore and offshore wind turbines and has also been considered in wave power applications. Wind power generation is the most mature renewable technology. However, wave energy has attracted a large interest recently as the potential for power extraction is very significant. Various wave energy converter (WEC) technologies currently exist with the oscillating water column (OWC) type converter being one of the most advanced. There are fundemental differences in the power profile of the pneumatic power supplied by the OWC WEC and that of a wind turbine and this causes significant challenges in the selection and rating of electrical generators for the OWC devises. The thesis initially aims to provide an accurate per-phase equivalent circuit model of the DFIM by investigating various characterisation testing procedures. Novel testing methodologies based on the series-coupling tests is employed and is found to provide a more accurate representation of the DFIM than the standard IEEE testing methods because the series-coupling tests provide a direct method of determining the equivalent-circuit resistances and inductances of the machine. A second novel method known as the extended short-circuit test is also presented and investigated as an alternative characterisation method. Experimental results on a 1.1 kW DFIM and a 30 kW DFIM utilising the various characterisation procedures are presented in the thesis. The various test methods are analysed and validated through comparison of model predictions and torque-versus-speed curves for each induction machine. Sensitivity analysis is also used as a means of quantifying the effect of experimental error on the results taken from each of the testing procedures and is used to determine the suitability of the test procedures for characterising each of the devices. The series-coupling differential test is demonstrated to be the optimum test. The research then focuses on the OWC WEC and the modelling of this device. A software model is implemented based on data obtained from a scaled prototype device situated at the Irish test site. Test data from the electrical system of the device is analysed and this data is used to develop a performance curve for the air turbine utilised in the WEC. This performance curve was applied in a software model to represent the turbine in the electro-mechanical system and the software results are validated by the measured electrical output data from the prototype test device. Finally, once both the DFIM and OWC WEC power take-off system have been modeled succesfully, an investigation of the application of the DFIM to the OWC WEC model is carried out to determine the electrical machine rating required for the pulsating power derived from OWC WEC device. Thermal analysis of a 30 kW induction machine is carried out using a first-order thermal model. The simulations quantify the limits of operation of the machine and enable thedevelopment of rating requirements for the electrical generation system of the OWC WEC. The thesis can be considered to have three sections. The first section of the thesis contains Chapters 2 and 3 and focuses on the accurate characterisation of the doubly-fed induction machine using various testing procedures. The second section, containing Chapter 4, concentrates on the modelling of the OWC WEC power-takeoff with particular focus on the Wells turbine. Validation of this model is carried out through comparision of simulations and experimental measurements. The third section of the thesis utilises the OWC WEC model from Chapter 4 with a 30 kW induction machine model to determine the optimum device rating for the specified machine. Simulations are carried out to perform thermal analysis of the machine to give a general insight into electrical machine rating for an OWC WEC device.

Engineering of metal oxide interfaces for renewable energy applications

Diminishing non-renewable energy resources and planet-wide de-pollution on our planet are among the major problems which mankind faces into the future. To solve these problems, renewable energy sources such as readily available and inexhaustible sunlight will have to be used. There are however no readily available photocatalysts that are photocatalytically active under visible light; it is well established that the band gap of the prototypical photocatalyst, titanium dioxide, is the UV region with the consequence that only 4% of sun light is utilized. For this reason, this PhD project focused on developing new materials, based on titanium dioxide, which can be used in visible light activated photocatalytic hydrogen production and destruction of pollutant molecules. The main goal of this project is to use simulations based on first principles to engineer and understand rationally, materials based on modifying TiO2 that will have the following properties: (1) a suitable band gap in order to increase the efficiency of visible light absorption, with a gap around 2 – 2.5 eV considered optimum. (2). The second key aspect in the photocatalytic process is electron and hole separation after photoexcitation, which enable oxidation/reduction reactions necessary to i.e. decompose pollutants. (3) Enhanced activity over unmodified TiO2. In this thesis I present results on new materials based on modifying TiO2 with supported metal oxide nanoclusters, from two classes, namely: transition metal oxides (Ti, Ni, Cu) and p-block metal oxides (Sn, Pb, Bi). We find that the deposited metal oxide nanoclusters are stable at rutile and anatase TiO2 surfaces and present an analysis of changes to the band gap of TiO2, identifying those modifiers that can change the band gap to the desirable range and the origin of this. A successful collaboration with experimental researchers in Japan confirms many of the simulation results where the origin of improved visible light photocatalytic activity of oxide nanocluster-modified TiO2 is now well understood. The work presented in this thesis, creates a road map for the design of materials with desired photocatalytic properties and contributes to better understanding these properties which are of great application in renewable energy utilization.