USE OF TSIM IN DIGITAL SIMULATION.

In this paper we shall discuss the use of the TSIM simulation software for modelling large-scale industrial processes. The discussion draws on our recent experience of modelling a large plant in the food-processing industry. We shall focus on those features of software organization and software engineering which proved to be particularly necessary for the execution of this project, and illustrate the extent to which the use of TISM facilitated the implementation of these features. We shall also make some general remarks about the 'life-cycle' of models resulting from projects of this kind.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (3rd report, direct numerical simulation on wall blocking effects for anisotropic turbulence)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method proposed in the second report using Townsend's "simple model eddies" for DNS was extended to generate axisymmetric anisotropic turbulence. A force field is obtained in real space by sprinkling many space-filling "simple model eddies" whose centers are randomly but uniformly distributed in space. The axes of rotation are controlled in this study to generate axisymmetric anisotropic turbulence. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects for anisotropic turbulence are investigated. The results show that stationary axisymmetric anisotropic turbulence is generated using the present method. Turbulence intensities near the wall showed good agreements with the rapid distortion theory (RDT) for small t (t ≪ TL), where TL. is the eddy turnover time. The splat effect (i. e. turbulence intensities of the components parallel to the surface are amplified) occurs near the boundary and the viscous effect attenuates the splat effect at the quasi steady state at low Reynolds number as for Isotropic turbulence. Prandtl's secondary flow of the second kind does not occur for low Reynolds number flows, which qualitatively agrees with previous observetion in a mixing-box.

Using Discrete Event Simulation (DES) to manage theatre operations in healthcare: An audit-based case study

This paper discusses the application of Discrete Event Simulation (DES) in modelling the complex relationship between patient types, case-mix and operating theatre allocation in a large National Health Service (NHS) Trust in London. The simulation model that was constructed described the main features of nine theatres, focusing on operational processes and patient throughput times. The model was used to test three scenarios of case-mix and to demonstrate the potential of using simulation modelling as a cost effective method for understanding the issues of healthcare operations management and the role of simulation techniques in problem solving. The results indicated that removing all day cases will reduce patient throughput by 23.3% and the utilization of the orthopaedic theatre in particular by 6.5%. This represents a case example of how DES can be used by healthcare managers to inform decision making. © 2008 IEEE.

Fast and compact simulation models for a variety of FET nano devices by the CMOS EKV equations

In this paper we explore the possibility of using the equations of a well known compact model for CMOS transistors as a parameterized compact model for a variety of FET based nano-technology devices. This can turn out to be a practical preliminary solution for system level architectural researchers, who could simulate behaviourally large scale systems, while more physically based models become available for each new device. We have used a four parameter version of the EKV model equations and verified that fitting errors are similar to those when using them for standard CMOS FET transistors. The model has been used for fitting measured data from three types of FET nano-technology devices obeying different physics, for different fabrication steps, and under different programming conditions. © 2009 IEEE NANO Organizers.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Turn-off failure mechanism in large area IGCTs

The destruction mechanism in large area IGCTs (Integrated Gate Commutated Thyristors) under inductive switching conditions is analyzed in detail. The three-dimensional nature of the turn-off process in a 91mm diameter wafer is simulated with a two-dimensional representation. Simulation results show that the final destruction is caused by the uneven dynamic avalanche current distribution across the wafer. © 2011 IEEE.

Statistical user simulation with a hidden agenda

Recent work in the area of probabilistic user simulation for training statistical dialogue managers has investigated a new agenda-based user model and presented preliminary experiments with a handcrafted model parameter set. Training the model on dialogue data is an important next step, but non-trivial since the user agenda states are not observable in data and the space of possible states and state transitions is intractably large. This paper presents a summary-space mapping which greatly reduces the number of state transitions and introduces a tree-based method for representing the space of possible agenda state sequences. Treating the user agenda as a hidden variable, the forward/backward algorithm can then be successfully applied to iteratively estimate the model parameters on dialogue data. © 2007 Association for Computational Linguistics.

Simulation studies on the magnetization of (RE)BCO bulk superconductors using various split-coil arrangements

Simulation studies were conducted on the magnetization of (RE)BCO (RE-Ba-Cu-O, where RE represents a rare earth element) bulk superconductors using various split-coil arrangements by solving the critical state equation using the commercial software FlexPDE. A pair of coaxial coils of identical size is identified as an optimum arrangement for practical magnetization at 77K by the zero-field cooling technique. In general, the magnetization process is likely to be most effective when the outer radius of the coils lies between 100% and 50% of the sample radius. A relatively large coil pair is necessary for samples with either a smaller aspect ratio or larger values of J c0. Two different regimes of flux penetration are found to be involved in the magnetization process. For a sufficiently small sample, the penetration field is determined by flux propagation from beneath the coil to the centre of the sample; for a sufficiently large sample, the definitive propagation route is from beneath the coil to the periphery of the sample. Effective split-coil magnetization occurs only in the former regime, and both penetration regimes are completely different from that involved in the solenoidal-coil magnetization process. © 2012 IOP Publishing Ltd.