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.

Numerical analysis of the settlement induced damage to Palazzo Loggia in Brescia

One of the main causes of failure of historic buildings is represented by the differential settlements of foundations. Finite element analysis provides a useful tool for predicting the consequences of given ground displacements in terms of structural damage and also assesses the need of strengthening techniques. The actual damage classification for buildings subject to settlement bases the assessment of the potential damage on the expected crack pattern of the structure. In this paper, the correlation between the physical description of the damage in terms of crack width and the interpretation of the finite element analysis output is analyzed. Different discrete and continuum crack models are applied to simulate an experiment carried on a scale model of a masonry historical building, the Loggia Palace in Brescia (Italy). Results are discussed and a modified version of the fixed total strain smeared crack model is evaluated, in order to solve the problem related to the calculation of the exact crack width.

Numerical analyses of tunnel-induced settlement damage to a masonry wall

Settlements due to underground construction represent a risk for the architectural heritage, especially in The Netherlands, because of the combination of soft soil, fragile pile foundation and brittle, un-reinforced masonry façade. Modelling of soil-structure interaction is fundamental to assess the risk of building damage due to tunnelling. This paper presents results of finite element analyses carried out with different models for a simple masonry wall. Focus is paid on the comparison between coupled, uncoupled and semi-coupled analyses, in which the soil-structure interaction is represented in different ways. In particular, the implementation of a soil-structure interface model in the numerical analyses is analysed, in order to asses its validity. The aim of the research project is the development of a damage classification system for different building typologies.

NUMERICAL ANALYSIS OF A NOVEL 980nm MICROCAVITY SEMICONDUCTOR OPTICAL AMPLIFIER WITH HIGH GAIN AND LOW NOISE

A novel microcavity semiconductor optical amplifier ( MCSOA) was proposed by incorporating top and bottom distributed Bragg reflectors ( DBRs) into the waveguide structure of conventional traveling-wave semiconductor optical amplifiers(TW-SOAs). The incoming( outgoing) light beam incidented onto (escaped from) the waveguide structure at a oblique angle through two optical windows, where the top DBR was etched away, and anti-reflection coating was deposited. The light beams inside the optical cavity were reflected repeatedly between two DBRs and propagated along waveguide in a zigzag optical path. The performance of the MCSOA was systematically investigated by extensive numerical simulation based on a traveling-wave model by taking into account the comprehensive effects of DBRs on both the amplification of signals and the filtering of spontaneous emission( SE). Our results show that the MCSOA is capable of achieving a fiber-to-fiber gain as high as 40dB and a low noise figure is less than 3.5dB.

Numerical analysis of gain saturation, noise figure, and carrier distribution for quantum-dot semiconductor-optical amplifiers

The gain saturation behaviors and noise figure are numerically analyzed for quantum-dot semiconductor optical amplifiers (QD-SOAs). The carrier and photon distributions in the longitudinal direction as well as the photon energy dependent facet reflectivity are accounted in the rate equations, which are solved with output amplified spontaneous emission spectrum as iterative variables. The longitudinal distributions of the occupation probabilities and spectral-hole burning are presented for electrons in the excited and ground states of quantum dots. The saturation output power 19.7 dBm and device gain 20.6 dB are obtained for a QD-SOA with the cavity length of 6 rum at the bias current of 500 mA. The influences of them electron intradot relaxation time and the QD capture time on the gain spectrum are simulated with the relaxation time of 1, 30, and 60 ps and capture time of 1, 5, and 10 ps. The noise figure as low as 3.5 dB is expected due to the strong polarization sensitive spontaneous emission. The characteristics of gain saturation and noise figure versus input signal power for QD-SOAs are similar to that of semiconductor. linear optical amplifiers with gain clamping by vertical laser fields.