Numerical stability of the predictor-corrector method in Monte Carlo burnup calculations of critical reactors

Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the simplest methods, such as the beginning-of-step and middle-of-step constant flux approximations, are numerically unstable in fuel cycle calculations of critical reactors. Here we show that even the predictor-corrector methods that are implemented in established Monte Carlo burnup codes can be numerically unstable in cycle calculations of large systems. © 2013 Elsevier Ltd. All rights reserved.

Coupled neutronic thermo-hydraulic analysis of full PWR core with Monte-Carlo based BGCore system

BGCore reactor analysis system was recently developed at Ben-Gurion University for calculating in-core fuel composition and spent fuel emissions following discharge. It couples the Monte Carlo transport code MCNP with an independently developed burnup and decay module SARAF. Most of the existing MCNP based depletion codes (e.g. MOCUP, Monteburns, MCODE) tally directly the one-group fluxes and reaction rates in order to prepare one-group cross sections necessary for the fuel depletion analysis. BGCore, on the other hand, uses a multi-group (MG) approach for generation of one group cross-sections. This coupling approach significantly reduces the code execution time without compromising the accuracy of the results. Substantial reduction in the BGCore code execution time allows consideration of problems with much higher degree of complexity, such as introduction of thermal hydraulic (TH) feedback into the calculation scheme. Recently, a simplified TH feedback module, THERMO, was developed and integrated into the BGCore system. To demonstrate the capabilities of the upgraded BGCore system, a coupled neutronic TH analysis of a full PWR core was performed. The BGCore results were compared with those of the state of the art 3D deterministic nodal diffusion code DYN3D (Grundmann et al.; 2000). Very good agreement in major core operational parameters including k-eff eigenvalue, axial and radial power profiles, and temperature distributions between the BGCore and DYN3D results was observed. This agreement confirms the consistency of the implementation of the TH feedback module. Although the upgraded BGCore system is capable of performing both, depletion and TH analyses, the calculations in this study were performed for the beginning of cycle state with pre-generated fuel compositions. © 2011 Published by Elsevier B.V.

The stochastic implicit Euler method - A stable coupling scheme for Monte Carlo burnup calculations

Existing Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the coupling schemes of the existing Monte Carlo burnup codes can be numerically unstable. Here we develop the Stochastic Implicit Euler method - a stable and efficient new coupling scheme. The implicit solution is obtained by the stochastic approximation at each time step. Our test calculations demonstrate that the Stochastic Implicit Euler method can provide an accurate solution to problems where the methods in the existing Monte Carlo burnup codes fail. © 2013 Elsevier Ltd. All rights reserved.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

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.

Spin-dependent transport and spin polarization in coupled quantum wells

We theoretically investigate the electron transport and spin polarization of two coupled quantum wells with Dresselhaus spin-orbit interaction. In analogy with the optical dual-channel directional coupler, the resonant tunneling effect is treated by the coupled-mode equations. We demonstrate that spin-up and -down electrons can be completely separated from each other for the system with an appropriate system geometry and a controllable barrier. Our result provides a new approach to construct spin-switching devices without containing any magnetic materials or applying a magnetic field. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2981204]

Multi-channel micro neural probe fabricated with SOI

Silicon-on-insulator (SOI) substrate is widely used in micro-electro-mechanical systems (MEMS). With the buried oxide layer of SOI acting as an etching stop, silicon based micro neural probe can be fabricated with improved uniformity and manufacturability. A seven-record-site neural probe was formed by inductive-coupled plasma (ICP) dry etching of an SOI substrate. The thickness of the probe is 15 mu m. The shaft of the probe has dimensions of 3 mmx100 mu mx15 mu m with typical area of the record site of 78.5 mu m(2). The impedance of the record site was measured in-vitro. The typical impedance characteristics of the record sites are around 2 M Omega at 1 kHz. The performance of the neural probe in-vivo was tested on anesthetic rat. The recorded neural spike was typically around 140 mu V. Spike from individual site could exceed 700 mu V. The average signal noise ratio was 7 or more.

Design and Fabrication of a Photonic Crystal Channel Drop Filter Based on an Asymmetric Silicon-on-Insulator Slab

We report on the design and fabrication of a photonic crystal (PC) channel drop filter based on an asymmetric silicon-on-insulator (SOI) slab. The filter is composed of two symmetric stick-shape micro-cavities between two single-line-defect (W1) waveguides in a triangular lattice, and the phase matching condition for the filter to improve the drop efficiency is satisfied by modifying the positions and radii of the air holes around the micro-cavities. A sample is then fabricated by using electron beam lithography (EBL) and inductively coupled plasma (ICP) etching processes. The measured 0 factor of the filter is about 1140, and the drop efficiency is estimated to be 73% +/- 5% by fitting the transmission spectrum.

Prediction and suppression of strong dispersive coupling in microracetrack channel drop filters

For a four-port microracetrack channel drop filter, unexpected transmission characteristics due to strong dispersive coupling are demonstrated by the light tunneling between the input-output waveguides and the resonator, where a large dropping transmission at off-resonance wavelengths is observed by finite-difference time-domain simulation. It causes a severe decline of the extinction ratio and finesse. An appropriate decrease of the coupling strength is found to suppress the dispersive coupling and greately increase the extinction ratio and finesse, a decreased coupling strength can be realized by the application of an asymmetrical coupling waveguide structure. In addition, the profile of the coupling dispersion in the transmission spectra can be predicted based on a coupled mode theory analysis of an equivalent system consisting of two coupling straight waveguides. The effects of structure parameters on the transmission spectra obtained by this method agree well with the numerical results. It is useful to avoid the strong dispersive coupling region in the filter design. (c) 2007 Optical Society of America.

Design and Fabrication of a Photonic Crystal Channel Drop Filter Based on an Asymmetric Silicon-on-Insulator Slab

We report on the design and fabrication of a photonic crystal (PC) channel drop filter based on an asymmetric silicon-on-insulator (SOI) slab. The filter is composed of two symmetric stick-shape micro-cavities between two single-line-defect (W1) waveguides in a triangular lattice, and the phase matching condition for the filter to improve the drop efficiency is satisfied by modifying the positions and radii of the air holes around the micro-cavities. A sample is then fabricated by using electron beam lithography (EBL) and inductively coupled plasma (ICP) etching processes. The measured 0 factor of the filter is about 1140, and the drop efficiency is estimated to be 73% +/- 5% by fitting the transmission spectrum.

The Fabrication of Eight-Channel DFB Laser Array Using Sampled Gratings

An eight-channel monolithically integrated complex-coupled distributed-feedback laser array based on sampled gratings has been designed and fabricated. Selective lasing at different wavelengths is obtained. The frequency separation between each adjacent channel is about 200 GHz. The typical threshold current is between 30 and 40 mA. The optical output power of each channel is about 10 mW at an injection current of 100 mA. The continuous tuning of emission wavelength with injected currents is also demonstrated.

Effective-mass theory for InAs/GaAs strained coupled quantum dots

In the framework of effective-mass envelope-function theory, the optical transitions of InAs/GaAs strained coupled quantum dots grown on GaAs (100) oriented substrates are studied. At the Gamma point, the electron and hole energy levels, the distribution of electron and hole wave functions along the growth and parallel directions, the optical transition-matrix elements, the exciton states, and absorption spectra are calculated. In calculations, the effects due to the different effective masses of electrons and holes in different materials are included. Our theoretical results are in good agreement with the available experimental data.

Quantum-mechanical calculations of vibrationally resolved cross sections for non-dissociative charge transfer of O3+ with H-2

Charge transfer due to collisions of ground state O3+ (2s(2)2p P-2) ions with molecular hydrogen is investigated using the quantum-mechanical molecular-orbital close-coupling (MOCC) method, and electronic and vibrational state-selective cross sections along with the corresponding differential cross sections are calculated for projectile energies of 100, 500, 1000 and 5000 eV/u at the orientation angles of 25 degrees,45 degrees and 89 degrees. The adiabatic potentials and radial coupling matrix elements utilized in the QMOCC calculations were obtained with the spin-coupled valence-bond approach. The infinite order sudden approximation (IOSA) and the vibrational sudden approximation (VSA) are utilized to deal with the rotation of H-2 and the coupling between the electron and the vibration of H-2. It is found that the distribution of vibrationally resolved cross sections with the vibrational quantum number upsilon' of H-2(+) (upsilon') varies with the increment of the projectile energy; and the electronic and vibrational stateselective differential cross sections show similar behaviors: there is a highest platform within a very small scattering angle, beyond which the differential cross sections decrease as the scattering angle increases and lots of oscillating structures appear, where the scattering angle of the first structure decreases as E-P(-1/2) with the increment of the projectile energy E-P; and the structure and amplitude of the differential cross sections are sensitive to the orientation of molecule H-2, which provides a possibility to identify the orientations of molecule H-2 by the vibrational state-selective differential scattering processes.

Distribution of uranium isotopes in the main channel of Yellow river (Huanghe), China

Uranium isotopes were measured in waters and suspended particulate matters (SPM) of the main channel of Yellow River, China that were sampled during four field trips between August 2005 and July 2006. The results show that the concentration of dissolved U (2.04-7.83 mu g/l) and the activity ratio of U-234/U-238 (1.36-1.67) are much higher than the average U concentrations and activity ratios of global major rivers. Mass balance calculations using the results of simulated experiments and measurement data show that the section of the Yellow River between Lanzhou and Sanmenxia has its dissolved U derived from two sources: suspended sediments (68%) and groundwater/runoff from loess deposits (32%). Both sources are related to the heavy erosion of the Chinese Loess Plateau. (C) 2008 Elsevier Ltd. All rights reserved.

Interface of chip-based capillary electrophoresis-inductively coupled plasma-atomic emission spectrometry (CE-ICP-AES)

An interface of chip-based capillary electrophoresis (CE)-inductively coupled plasma-atomic emission spectrometry (ICP-AES) that is based on cross-flow nebulization has been developed. A polydimethylsiloxane (PDMS) CE-chip with conventional cross channel layout was used. A stainless steel tube was placed orthogonal to the exit of the CE separation channel for cross flow nebulization. A supplementary flow of buffer solution at the channel exit was used to improve nebulization efficiency. Two capillaries were inserted into the CE chip near the inlet of the separation channel for sample and buffer solution injection. Syringe pumps were used to manipulate the flow rate and flow direction of the sample, buffer, and supplementary buffer solution. Peak broadening due to the shape (bulb and tube-shaped) and size of the spray chambers was studied. The smaller tube-shaped spray chamber was used because of smaller peak broadening effect due to aerosol transport. The nebulization and transport efficiency of the CE-ICP interface was approximately 10%. Ba2+ and Mg2+ ions were eluted from the CE-chip within 30 s. Resolution of the Ba2+ and Mg2+ peaks was 0.7 using the chip-based CE-ICP-AES system.