SBLOCK: A framework for efficient stencil-based PDE solvers on multi-core platforms

We present a new software framework for the implementation of applications that use stencil computations on block-structured grids to solve partial differential equations. A key feature of the framework is the extensive use of automatic source code generation which is used to achieve high performance on a range of leading multi-core processors. Results are presented for a simple model stencil running on Intel and AMD CPUs as well as the NVIDIA GT200 GPU. The generality of the framework is demonstrated through the implementation of a complete application consisting of many different stencil computations, taken from the field of computational fluid dynamics. © 2010 IEEE.

Efficient generation of one-group cross sections for coupled Monte Carlo depletion calculations

Coupled Monte Carlo depletion systems provide a versatile and an accurate tool for analyzing advanced thermal and fast reactor designs for a variety of fuel compositions and geometries. The main drawback of Monte Carlo-based systems is a long calculation time imposing significant restrictions on the complexity and amount of design-oriented calculations. This paper presents an alternative approach to interfacing the Monte Carlo and depletion modules aimed at addressing this problem. The main idea is to calculate the one-group cross sections for all relevant isotopes required by the depletion module in a separate module external to Monte Carlo calculations. Thus, the Monte Carlo module will produce the criticality and neutron spectrum only, without tallying of the individual isotope reaction rates. The onegroup cross section for all isotopes will be generated in a separate module by collapsing a universal multigroup (MG) cross-section library using the Monte Carlo calculated flux. Here, the term "universal" means that a single MG cross-section set will be applicable for all reactor systems and is independent of reactor characteristics such as a neutron spectrum; fuel composition; and fuel cell, assembly, and core geometries. This approach was originally proposed by Haeck et al. and implemented in the ALEPH code. Implementation of the proposed approach to Monte Carlo burnup interfacing was carried out through the BGCORE system. One-group cross sections generated by the BGCORE system were compared with those tallied directly by the MCNP code. Analysis of this comparison was carried out and led to the conclusion that in order to achieve the accuracy required for a reliable core and fuel cycle analysis, accounting for the background cross section (σ0) in the unresolved resonance energy region is essential. An extension of the one-group cross-section generation model was implemented and tested by tabulating and interpolating by a simplified σ0 model. A significant improvement of the one-group cross-section accuracy was demonstrated.

The rhythm of fountains: The length and time scales of rise height fluctuations at low and high Froude numbers

The magnitude and frequency of vertical fluctuations of the top of an axisymmetric miscible Boussinesq fountain forms the focus of this work. We present measurements of these quantities for saline-aqueous fountains in uniform quiescent surroundings. Our results span source Froude numbers 0.3 ≤ Fr 0 ≤ 40 and, thereby, encompass very weak, weak, intermediate and forced classes of fountain. We identify distinct scalings, based on known quantities at the fountain source, for the frequency of fountain height fluctuations which collapse our data within bands of Fr0. Notably, our scalings reveal that the (dimensionless) frequency takes a constant value within each band. These results highlight characteristic time scales for the fluctuations which we decompose into a single, physically apparent, length scale and velocity scale within each band. Moreover, within one particular band, spanning source Froude numbers towards the lower end of the full range considered, we identify unexpectedly long-period fluctuations indicating a near balance of inertia and (opposing) buoyancy at the source. Our analysis identifies four distinct classes of fluctuation behaviour (four bands of Fr 0) and this classification matches well with existing classifications of fountains based on rise heights. As such, we show that an analysis of the behaviour of the fountain top alone, rather than the entire fountain, provides an alternative approach to classifying fountains. The similarity of classifications based on the two different methods confirms that the boundaries between classes mark tangible changes in the physics of fountains. For high Fr0 we show that the dominant fluctuations occur at the scale of the largest eddies which can be contained within the fountain near its top. Extending this, we develop a Strouhal number, Strtop, based on experimental measures of the fountain top, defined such that Strtop = 1 would suggest the dominant fluctuations are caused by a continual cycle of eddies forming and collapsing at this largest physical scale. For high- Fr 0 fountains we find Strtop ≈ 0. 9. © 2013 Cambridge University Press.

The effect of beam interruptions on the integrity of ADSR fuel pin cladding: A thermo-mechanical analysis

During its lifetime in the core, the cladding of an Accelerator Driven Subcritical Reactor (ADSR) fuel pin is expected to experience variable stresses due to frequent interruptions in the accelerator proton beam. This paper investigates the thermal fatigue damage in the cladding due to repetitive and unplanned beam interruptions under certain operational conditions. Beam trip data was obtained for four operating high power proton accelerators, among which the Spallation Neutron Source (SNS) superconducting accelerator was selected for further analysis. 9Cr-1Mo-Nb-V (T91) steel was selected as the cladding material because of its proven compatibility with proposed ADSR design concepts. The neutronic, thermal and stress analyses were performed using the PTS-ADS, a code that has been specifically developed for studying the dynamic response to beam-induced transients in accelerator driven subcritical systems. The lifetime of the fuel cladding in the core was estimated for three levels of allowed pin power and specific operating conditions. © 2012 Elsevier Ltd. All rights reserved.

Achieving Csiszár's source-channel coding exponent with product distributions

We derive a random-coding upper bound on the average probability of error of joint source-channel coding that recovers Csiszár's error exponent when used with product distributions over the channel inputs. Our proof technique for the error probability analysis employs a code construction for which source messages are assigned to subsets and codewords are generated with a distribution that depends on the subset. © 2012 IEEE.

An extension of slow feature analysis for nonlinear blind source separation

We present and test an extension of slow feature analysis as a novel approach to nonlinear blind source separation. The algorithm relies on temporal correlations and iteratively reconstructs a set of statistically independent sources from arbitrary nonlinear instantaneous mixtures. Simulations show that it is able to invert a complicated nonlinear mixture of two audio signals with a high reliability. The algorithm is based on a mathematical analysis of slow feature analysis for the case of input data that are generated from statistically independent sources. © 2014 Henning Sprekeler, Tiziano Zito and Laurenz Wiskott.

Properties of amorphous carbon-silicon alloys deposited by a high plasma density source

The addition of silicon to hydrogenated amorphous carbon can have the advantageous effect of lowering the compressive stress, improving the thermal stability of its hydrogen, and maintaining a low friction coefficient up to high humidity. Most experiments to date have been on hydrogenated amorphous carbon-silicon alloys (a-C1-xSix:H) deposited by rf plasma enhanced chemical vapor deposition. This method gives alloys with sizeable hydrogen content and only moderate hardness. Here we use a high plasma density source known as the electron cyclotron wave resonance source to prepare films with higher sp3 content and lower hydrogen content. The composition and bonding in the alloys is determined by x-ray photoelectron spectroscopy, Rutherford backscattering, elastic recoil detection analysis, visible and ultraviolet (UV) Raman spectroscopy, infrared spectroscopy, and x-ray reflectivity. We find that it is possible to produce relatively hard, low stress, low friction, almost humidity insensitive a-C1-xSix:H alloys with a good optical transparency and a band gap well over 2.5 eV. The friction behavior and friction mechanism of these alloys are studied and compared with that of a-C:H, ta-C:H, and ta-C. We show how UV Raman spectroscopy allows the direct detection of Si-C, Si-Hx, and C-Hx vibrations, not seen in visible Raman spectra. © 2001 American Institute of Physics.

Coupled deformation-flow analysis for methane hydrate production by depressurized wells

A formulation for coupled flow-deformation analysis of methane-hydrate extraction problems is presented. By assuming that the hydrate does not flow, a two phase flow formulation is considered, based on Darcy's law and capillary pressure relation. The formulation is implemented in the finite difference code FLAC. The code was used to investigate the stability of a methane extraction well by depressurizing the well. © 2005 Taylor & Francis Group, London.

On-line topographic analysis in the SEM

A novel method for on-line topographic analysis of rough surfaces in the SEM has been investigated. It utilises a digital minicomputer configured to act as a programmable scan generator and automatic focusing unit. The computer is coupled to the microscope through digital-to-analogue converters which enable it to generate ramp waveforms allowing the beam to be scanned over a small sub-region of the field under program control. A further digital-to-analogue converter regulates the current supply to the objective lens of the microscope. The video signal is sampled by means of an analogue-to-digital converter and the resultant binary code stored in the computer's memory as an array of numbers describing relative image intensity. Computations based on the intensity gradient of the image allow the objective lens current to be found for the in-focus condition, which may be related to the working distance through a previous calibration experiment. The sensitivity of the method for detecting small height changes is theoretically of the order of 1 μm. In practice the operator specifies features of interest by means of a mobile spot cursor injected into the SEM display screen, or he may scan the specimen at sub-regions corresponding to pre-determined points on a regular grid defined by him. The operation then proceeds under program control. | A novel method for on-line topographic analysis of rough surfaces in the SEM has been investigated. It utilizes a digital minicomputer configured to act as a programmable scan generator and automatic focusing unit. A further digital-to-analog converter regulates the current supply to the objective lens of the microscope. The video signal is sampled by means of an analog-to-digital converter and the resultant binary code stored in the computer's memory as an array of numbers describing relative image intensity. The sensitivity of the method for detecting small height changes is theroretically of the order of 1 mu m.

Advanced DC link capacitor technology application for a stiff voltage-source inverter

Problem of DC link size in a stiff voltage-source inverter for electric drive is described in the paper. Advantages of advanced film capacitor technology over conventional one for DC link application are reviewed. Conventional DC link capacitor selection methods are questioned in view of advanced capacitor technology utilization in stiff voltage-source inverter. For capacitor selection maximum ripple rms current point is shown. DC link ripple current spectrum analysis under modern PWM techniques is presented. Some capacitor selection recommendations are given. The analysis has been aided greatly by computer modeling in PSpice. ©2005 IEEE.

Eccentricity effects on acoustic radiation from a spherical source suspended within a thermoviscous fluid sphere.

Acoustic radiation from a spherical source undergoing angularly periodic axisymmetric harmonic surface vibrations while eccentrically suspended within a thermoviscous fluid sphere, which is immersed in a viscous thermally conducting unbounded fluid medium, is analyzed in an exact fashion. The formulation uses the appropriate wave-harmonic field expansions along with the translational addition theorem for spherical wave functions and the relevant boundary conditions to develop a closed-form solution in form of infinite series. The analytical results are illustrated with a numerical example in which the vibrating source is eccentrically positioned within a chemical fluid sphere submerged in water. The modal acoustic radiation impedance load on the source and the radiated far-field pressure are evaluated and discussed for representative values of the parameters characterizing the system. The proposed model can lead to a better understanding of dynamic response of an underwater acoustic lens. It is equally applicable in miniature transducer analysis and design with applications in medical ultrasonics.

Dynamic analysis of the brushless doubly-fed induction generator during symmetrical three-phase voltage dips

The Brushless Doubly-Fed Induction Generator (BDFIG) shows commercial promise as replacement for doublyfed slip-ring generators for wind power applications by offering reduced capital and operational costs due to its brushless operation. In order to facilitate its commercial deployment, the capabilities of the BDFIG system to comply with grid code requirements have to be assessed. This paper, for the first time, studies the performance of the BDFIG under grid fault ride-through and presents the dynamic behaviour of the machine during three-phase symmetrical voltage dips. Both full and partial voltage dips are studied using a vector model. Simulation and experimental results are provided for a 180 frame BDFIG.

Model structure and uncertainty for stochastic non-point source modelling applications

The uncertainty associated with a rainfall-runoff and non-point source loading (NPS) model can be attributed to both the parameterization and model structure. An interesting implication of the areal nature of NPS models is the direct relationship between model structure (i.e. sub-watershed size) and sample size for the parameterization of spatial data. The approach of this research is to find structural limitations in scale for the use of the conceptual NPS model, then examine the scales at which suitable stochastic depictions of key parameter sets can be generated. The overlapping regions are optimal (and possibly the only suitable regions) for conducting meaningful stochastic analysis with a given NPS model. Previous work has sought to find optimal scales for deterministic analysis (where, in fact, calibration can be adjusted to compensate for sub-optimal scale selection); however, analysis of stochastic suitability and uncertainty associated with both the conceptual model and the parameter set, as presented here, is novel; as is the strategy of delineating a watershed based on the uncertainty distribution. The results of this paper demonstrate a narrow range of acceptable model structure for stochastic analysis in the chosen NPS model. In the case examined, the uncertainties associated with parameterization and parameter sensitivity are shown to be outweighed in significance by those resulting from structural and conceptual decisions. © 2011 Copyright IAHS Press.