A natural extension of the conventional finite volume method into polygonal unstructured meshes for CFD application.

A new general cell-centered solution procedure based upon the conventional control or finite volume (CV or FV) approach has been developed for numerical heat transfer and fluid flow which encompasses both structured and unstructured meshes for any kind of mixed polygon cell. Unlike conventional FV methods for structured and block structured meshes and both FV and FE methods for unstructured meshes, the irregular control volume (ICV) method does not require the shape of the element or cell to be predefined because it simply exploits the concept of fluxes across cell faces. That is, the ICV method enables meshes employing mixtures of triangular, quadrilateral, and any other higher order polygonal cells to be exploited using a single solution procedure. The ICV approach otherwise preserves all the desirable features of conventional FV procedures for a structured mesh; in the current implementation, collocation of variables at cell centers is used with a Rhie and Chow interpolation (to suppress pressure oscillation in the flow field) in the context of the SIMPLE pressure correction solution procedure. In fact all other FV structured mesh-based methods may be perceived as a subset of the ICV formulation. The new ICV formulation is benchmarked using two standard computational fluid dynamics (CFD) problems i.e., the moving lid cavity and the natural convection driven cavity. Both cases were solved with a variety of structured and unstructured meshes, the latter exploiting mixed polygonal cell meshes. The polygonal mesh experiments show a higher degree of accuracy for equivalent meshes (in nodal density terms) using triangular or quadrilateral cells; these results may be interpreted in a manner similar to the CUPID scheme used in structured meshes for reducing numerical diffusion for flows with changing direction.

Modelling of industrial multi-physics processes—a key role for computational mechanics.

In the analysis of industrial processes, there is an increasing emphasis on systems governed by interacting continuum phenomena. Mathematical models of such multi-physics processes can only be achieved for practical simulations through computational solution procedures—computational mechanics. Examples of such multi-physics systems in the context of metals processing are used to explore some of the key issues. Finite-volume methods on unstructured meshes are proposed as a means to achieve efficient rapid solutions to such systems. Issues associated with the software design, the exploitation of high performance computers, and the concept of the virtual computational-mechanics modelling laboratory are also addressed in this context.

Computer modelling of human behaviour in aircraft fire accidents

The mathematical simulation of the evacuation process has a wide and largely untapped scope of application within the aircraft industry. The function of the mathematical model is to provide insight into complex behaviour by allowing designers, legislators, and investigators to ask ‘what if’ questions. Such a model, EXODUS, is currently under development, and this paper describes its evolution and potential applications. EXODUS is an egress model designed to simulate the evacuation of large numbers of individuals from an enclosure, such as an aircraft. The model tracks the trajectory of each individual as they make their way out of the enclosure or are overcome by fire hazards, such as heat and toxic gases. The software is expert system-based, the progressive motion and behaviour of each individual being determined by a set of heuristics or rules. EXODUS comprises five core interacting components: (i) the Movement Submodel — controls the physical movement of individual passengers from their current position to the most suitable neighbouring location; (ii) the Behaviour Submodel — determines an individual's response to the current prevailing situation; (iii) the Passenger Submodel — describes an individual as a collection of 22 defining attributes and variables; (iv) the Hazard Submodel — controls the atmospheric and physical environment; and (v) the Toxicity Submodel — determines the effects on an individual exposed to the fire products, heat, and narcotic gases through the Fractional Effective Dose calculations. These components are briefly described and their capabilities and limitations are demonstrated through comparison with experimental data and several hypothetical evacuation scenarios.

The prediction of fire propagation in enclosure fires

In this paper we present some early work concerned with the development of a simple solid fuel combustion model incorporated within a Computational Fluid Dynamics (CFD) framework. The model is intended for use in engineering applications of fire field modelling and represents an extension of this technique to situations involving the combustion of solid cellulosic hels A simple solid &el combustion model consisting of a thermal pyrolysis model, a six flux radiation model and an eddydissipation model for gaseous combustion have been developed and implemented within the CFD code CFDS-FLOW3D The model is briefly described and demonstrated through two applications involving fire spread in a compartment with a plywood lined ceiling. The two scenarios considered involve a fire in an open and closed compartment The model is shown to be able to qualitatively predict behaviours similar to flashover - in the case of the open room - and backdrafl - in the case of the initially closed room.

Numerical optimisation of water curtain barriers design

The interactions between water curtain protective barriers and impinging cold gas clouds released during industrial accidents are simulated by a 2D finite difference code using a PSI-cell of the heat, mass and momentum transfer processes between the phases. A consistent derivation of the continuous phase source terms is presented. The results of early simulations of an existing 1:10 experimental model of a water curtain impacted by a cold gas cloud are presented for two typical curtain configurations.

A novel neural/symbolic hybrid approach to heuristically optimized fuel allocation and automated revision of heuristics in nuclear engineering

FUELCON is an expert system for optimized refueling design in nuclear engineering. This task is crucial for keeping down operating costs at a plant without compromising safety. FUELCON proposes sets of alternative configurations of allocation of fuel assemblies that are each positioned in the planar grid of a horizontal section of a reactor core. Results are simulated, and an expert user can also use FUELCON to revise rulesets and improve on his or her heuristics. The successful completion of FUELCON led this research team into undertaking a panoply of sequel projects, of which we provide a meta-architectural comparative formal discussion. In this paper, we demonstrate a novel adaptive technique that learns the optimal allocation heuristic for the various cores. The algorithm is a hybrid of a fine-grained neural network and symbolic computation components. This hybrid architecture is sensitive enough to learn the particular characteristics of the ‘in-core fuel management problem’ at hand, and is powerful enough to use this information fully to automatically revise heuristics, thus improving upon those provided by a human expert.

A two-dimensional numerical investigation of the oscillatory flow behaviour in rectangular fire compartments with a single horizontal ceiling vent

This paper presents a comparison of fire field model predictions with experiment for the case of a fire within a compartment which is vented (buoyancydriven) to the outside by a single horizontal ceiling vent. Unlike previous work, the mathematical model does not employ a mixing ratio to represent vent temperatures but allows the model to predict vent temperatures a priori. The experiment suggests that the flow through the vent produces oscillatory behaviour in vent temperatures with puffs of smoke emerging from the fire compartment. This type of flow is also predicted by the fire field model. While the numerical predictions are in good qualitative agreement with observations, they overpredict the amplitudes of the temperature oscillations within the vent and also the compartment temperatures. The discrepancies are thought to be due to three-dimensional effects not accounted for in this model as well as using standard ‘practices’ normally used by the community with regards to discretization and turbulence models. Furthermore, it is important to note that the use of the k–ε turbulence model in a transient mode, as is used here, may have a significant effect on the results. The numerical results also suggest that a linear relationship exists between the frequency of vent temperature oscillation (n) and the heat release rate (Q0) of the type n∝Q0.290, similar to that observed for compartments with two horizontal vents. This relationship is predicted to occur only for heat release rates below a critical value. Furthermore, the vent discharge coefficient is found to vary in an oscillatory fashion with a mean value of 0.58. Below the critical heat release rate the mean discharge coefficient is found to be insensitive to fire size.

An investigation of the aspects of occupant behavior required for evacuation modeling

This article examines occupant behavior exhibited during evacuation conditions. This is based on a review of a wide range of published literature concerned with evacuation. Factors influencing evacuation performance can be categorized into four broad areas, namely, configurational, environmental, procedural, and, most importantly, behavioral. The contributory factors associated with each of the four influencing categories are examined in detail and it is suggested that these factors should be represented within evacuation models.

Multiphysics modelling of the metals casting process

Metals casting is a process governed by the interaction of a range of physical phenomena. Most computational models of this process address only what are conventionally regarded as the primary phenomena-heat conduction and solidification. However, to predict the formation of porosity (a factor of key importance in cast quality) requires the modelling of the interaction of the fluid flow, heat transfer, solidification and the development of stress-deformation in the solidified part of a component. In this paper, a model of the casting process is described which addresses all the main continuum phenomena involved in a coupled manner. The model is solved numerically using novel finite volume unstructured mesh techniques, and then applied to both the prediction of shape deformation (plus the subsequent formation of a gap at the metal-mould interface and its impact on the heat transfer behaviour) and porosity formation in solidifying metal components. Although the porosity prediction model is phenomenologically simplistic it is based on the interaction of the continuum phenomena and yields good comparisons with available experimental results. This work represents the first of the next generation of casting simulation tools to predict aspects of the structure of cast components.

The integration of an intelligent knowledge-based system into engineering software using the blackboard structure

Over recent years there has been an increase in the use of generic Computational Fluid Dynamics (CFD) software packages spread across various application fields. This has created the need for the integration of expertise into CFD software. Expertise can be integrated into CFD software in the form of an Intelligent Knowledge-Based System (IKBS). The advantages of integrating intelligence into generic engineering software are discussed with a special view to software engineering considerations. The software modelling cycle of a typical engineering problem is identified and the respective expertise and user control needed for each modelling phase is shown. The requirements of an IKBS for CFD software are discussed and compared to current practice. The blackboard software architecture is presented. This is shown to be appropriate for the integration of an IKBS into an engineering software package. This is demonstrated through the presentation of the prototype CFD software package FLOWES.

Magnetic modelling and magneto-hydro-dynamic simulations of an aluminium production electrolytic cell

A simulation of the motion of molten aluminium inside an electrolytic cell is presented. Since the driving term of the aluminium motion is the Lorentz (j × B) body force acting within the fluid,this problem involves the solution of the magneto-hydro-dynamic equations. Different solver modules for the magnetic field computation and for the fluid motion simulation are coupled together. The interactions of all these are presented and discussed.

Estimating the flow rate capacity of an overturned rail carriage end exit in the presence of smoke

While incidents requiring the rapid egress of passengers from trains are infrequent, perhaps the most challenging scenario for passengers involves the evacuation from an overturned carriage subjected to fire. In this paper we attempt to estimate the flow rate capacity of an overturned rail carriage end exit. This was achieved through two full-scale evacuation experiments, in one of which the participants were subjected to non-toxic smoke. The experiments were conducted as part of a pilot study into evacuation from rail carriages. In reviewing the experimental results, it should be noted that only a single run of each trial was undertaken with a limited — though varied — population. As a result it is not possible to test the statistical significance of the evacuation times quoted and so the results should be treated as indicative rather than definitive. The carriage used in the experiments was a standard class Mark IID which, while an old carriage design, shares many features with those carriages commonly found on the British rail network. In the evacuation involving smoke, the carriage end exit was found to achieve an average flow rate capacity of approximately 5.0 persons/min. The average flow rate capacity of the exit without smoke was found to be approximately 9.2 persons/min. It was noted that the presence of smoke tended to reduce significantly the exit flow rate. Due to the nature of the experimental conditions, these flow rates are considered optimistic. Finally, the authors make several recommendations for improving survivability in rail accidents. Copyright © 2000 John Wiley & Sons, Ltd.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

The effect of curing on the performance of ACF bonded chip-on-flex assemblies after thermal ageing

Purpose – Anisotropic conductive film (ACF) is now an attractive technology for direct mounting of chips onto the substrate as an alternative to lead-free solders. However, despite its various advantages over other technologies, it also has many unresolved reliability issues. For instance, the performance of ACF assembly in high temperature applications is questionable. The purpose of this paper is to study the effect of bonding temperatures on the curing of ACFs, and their mechanical and electrical performance after high temperature ageing. Design/methodology/approach – In the work presented in this paper, the curing degree of an ACF at different bonding temperatures was measured using a differential scanning calorimeter. The adhesion strength and the contact resistance of ACF bonded chip-on-flex assembly were measured before and after thermal ageing and the results were correlated with the curing degree of ACF. The ACF was an epoxy-based adhesive in which Au-Ni coated polymer particles were randomly dispersed. Findings – The results showed that higher bonding temperatures had resulted in better ACF curing and stronger adhesion. After ageing, the adhesion strength increased for the samples bonded at lower temperatures and decreased for the samples bonded at higher temperatures. ACF assemblies with higher degrees of curing showed smaller increases in contact resistance after ageing. Conduction gaps at the bump-particle and/or particle-pad interfaces were found with the help of scanning electron microscopy and are thought to be the root cause of the increase in contact resistance. Originality/value – The present study focuses on the effect of bonding temperatures on the curing of ACFs, and their adhesion strength and electrical performances after high temperature ageing. The results of this study may help the development of ACFs with higher heat resistance, so that ACFs can be considered as an alternative to lead-free solders.