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.

A new approach to the analytical and numerical form-finding of tensegrity structures

We develop a new formulation for the form-finding of tensegrity structures in which the primary variables are the Cartesian components of element lengths. Both an analytical and a numerical implementation of the formulation are described; each require a description of the connectivity of the tensegrity, with the iterative numerical method also requiring a random starting vector of member force densities. The analytical and numerical form-finding of tensegrity structures is demonstrated through six examples, and the results obtained are compared and contrasted with those available in the literature to verify the accuracy and viability of the suggested methods. © 2013 Elsevier Ltd. All rights reserved.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

An investigation into numerical analysis alternatives for predicting reingestion in turbine disc rim cavities

Reliable means of predicting ingestion in cavities adjacent to the main gas path are increasingly being sought by engineers involved in the design of gas turbines. In this paper, analysis is to be presented that results from an extended research programme, MAGPI, sponsored by the EU and several leading gas turbine manufactures and universities. Extensive use is made of CFD modelling techniques to understand the aerodynamic behaviour of a turbine stator well cavity, focusing on the interaction of cooling air supply with the main annulus gas. The objective of the study has been to benchmark a number of CFD codes and numerical techniques covering RANS and URANS calculations with different turbulence models in order to assess the suitability of the standard settings used in the industry for calculating the mechanics of the flow travelling between cavities in a turbine through the main gas path. The modelling methods employed have been compared making use of experimental data gathered from a dedicated two-stage turbine rig, running at engine representative conditions. Extensive measurements are available for a range of flow conditions and alternative cooling arrangements. The limitations of the numerical methods in calculating the interaction of the cooling flow egress and the main stream gas, and subsequent ingestion into downstream cavities in the engine (i.e. re-ingestion), have been exposed. This has been done without losing sight of the validation of the CFD for its use for predicting heat transfer, which was the main objective of the partners of the MAGPI Work- Package 1 consortium. Copyright © 2012 by ASME.

Efficient numerical models for investigating the statistics of random dynamic systems

The study of random dynamic systems usually requires the definition of an ensemble of structures and the solution of the eigenproblem for each member of the ensemble. If the process is carried out using a conventional numerical approach, the computational cost becomes prohibitive for complex systems. In this work, an alternative numerical method is proposed. The results for the response statistics are compared with values obtained from a detailed stochastic FE analysis of plates. The proposed method seems to capture the statistical behaviour of the response with a reduced computational cost.

A numerical study into the effects of Bio-synthetic Paraffinic Kerosine blends with Jet-A fuel for civil aircraft engine

Growing concerns regarding fluctuating fuel costs and pollution targets for gas emissions, have led the aviation industry to seek alternative technologies to reduce its dependency on crude oil, and its net emissions. Recently blends of bio-fuel with kerosine, have become an alternative solution as they offer "greener" aircraft and reduce demand on crude oil. Interestingly, this technique is able to be implemented in current aircraft as it does not require any modification to the engine. Therefore, the present study investigates the effect of blends of bio-synthetic paraffinic kerosine with Jet-A in a civil aircraft engine, focusing on its performance and exhaust emissions. Two bio-fuels are considered: Jatropha Bio-synthetic Paraffinic Kerosine (JSPK) and Camelina Bio-synthetic Paraffinic Kerosine (CSPK); there are evaluated as pure fuels, and as 10% and 50% blend with Jet-A. Results obtained show improvement in thrust, fuel flow and SFC as composition of bio-fuel in the blend increases. At design point condition, results on engine emissions show reduction in NO x, and CO, but increases of CO is observed at fixed fuel condition, as the composition of bio-fuel in the mixture increases. Copyright © 2012 by ASME.

Numerical analysis of shallow foundations under combined horizontal and torsional loading

Foundations of subsea infrastructure in deep water subjected to asymmetric environmental loads have underscored the importance of combined torsional and horizontal loading effects on the bearing capacity of rectangular shallow foundations. The purpose of this study is to investigate the undrained sliding and torsional bearing capacity of rectangular and square shallow foundations together with the interaction response under combined loading using three-dimensional finite element (3D-FE) analysis. Upper bound plastic limit analysis is employed to establish a reference value for horizontal and torsional bearing capacity, and an interaction relationship for the combined loading condition. Satisfactory agreement of plastic limit analysis (PLA) and 3D-FE results for ultimate capacity and interaction curves ensures that simple PLA solution could be used to evaluate the bearing capacity problem of foundation under combined sliding and torsion.

Numerical modelling of chirality-Induced bi-Directional swimming of artificial flagella

Biomimetic micro-swimmers can be used for various medical applications, such as targeted drug delivery and micro-object (e.g. biological cells) manipulation, in lab-on-a-chip devices. Bacteria swim using a bundle of flagella (flexible hair-like structures) that form a rotating cork-screw of chiral shape. To mimic bacterial swimming, we employ a computational approach to design a bacterial (chirality-induced) swimmer whose chiral shape and rotational velocity can be controlled by an external magnetic field. In our model, we numerically solve the coupled governing equations that describe the system dynamics (i.e. solid mechanics, fluid dynamics and magnetostatics). We explore the swimming response as a function of the characteristic dimensionless parameters and put special emphasis on controlling the swimming direction. Our results provide fundamental physical insight on the chirality-induced propulsion, and it provides guidelines for the design of magnetic bi-directional micro-swimmers. © 2013 The Author(s) Published by the Royal Society. All rights reserved.

Numerical simulations supporting the process design of ring rolling processes

In conventional Finite Element Analysis (FEA) of radial-axial ring rolling (RAR) the motions of all tools are usually defined prior to simulation in the preprocessing step. However, the real process holds up to 8 degrees of freedom (DOF) that are controlled by industrial control systems according to actual sensor values and preselected control strategies. Since the histories of the motions are unknown before the experiment and are dependent on sensor data, the conventional FEA cannot represent the process before experiment. In order to enable the usage of FEA in the process design stage, this approach integrates the industrially applied control algorithms of the real process including all relevant sensors and actuators into the FE model of ring rolling. Additionally, the process design of a novel process 'the axial profiling', in which a profiled roll is used for rolling axially profiled rings, is supported by FEA. Using this approach suitable control strategies can be tested in virtual environment before processing. © 2013 AIP Publishing LLC.

Numerical study of the impact of process variations on the motional resistance of weakly coupled MEMS resonators

This paper presents a numerical study of the impact of process-induced variations on the achievable motional resistance Rx of one-dimensional, cyclic and cross-coupled architectures of electrostatically transduced MEMS resonators operating in the 250 kHz range. Monte Carlo numerical simulations which accounted for up to 0.75% variation in critical resonator feature sizes were initiated on 1, 2, 3, 4, 5 and 9 coupled MEMS resonators for three distinct coupling architectures. Improvements of 100X in the spread of Rx and 2.7X in mean achievable Rx are reported for the case of 9 resonators when implemented in the cross-coupled topology, as opposed to the traditional one-dimensional chain. © 2013 IEEE.