World in motion

Professor Ed Galea CEng, MIFireE provides a welcome to Pedestrian and Evacuation Dynamics 2003, (PED 2003) to be held in London on 20-22 August 2003.

Predicting the past

H. Jiang, S. Gwynne, E.R. Galea, P. Lawrence, F. Jia and H. Ingason model a disco fire in Gothenburg, Sweden to compare the simulation’s predictions with actual events

The airEXODUS evacuation model and its application to aircraft safety

Computer based mathematical models describing the aircraft evacuation process have a vital role to play in the design and development of safer aircraft, the implementation of safer and more rigorous certification criteria, in cabin crew training and post-mortem accident investigation. As the risk of personal injury and the costs involved in performing full-scale certification trials are high, the development and use of these evacuation modelling tools are essential. Furthermore, evacuation models provide insight into the evacuation process that is impossible to derive from a single certification trial. The airEXODUS evacuation model has been under development since 1989 with support from the UK CAA and the aviation industry. In addition to describing the capabilities of the airEXODUS evacuation model, this paper describes the findings of a recent CAA project aimed at investigating model accuracy in predicting past certification trials. Furthermore, airEXODUS is used to examine issues related to the Blended Wing Body (BWB) and Very Large Transport Aircraft (VLTA). These radical new aircraft concepts pose considerable challenges to designers, operators and certification authorities. BWB concepts involving one or two decks with possibly four or more aisles offer even greater challenges. Can the largest exits currently available cope with passenger flow arising from four or five aisles? Do we need to consider new concepts in exit design? Should the main aisle be made wider to accommodate more passengers? In this paper we discuss various issues evacuation related issues associated VLTA and BWB aircraft and demonstrate how computer based evacuation models can be used to investigage these issues through examination of aisle/exit configurations for BWB cabin layouts.

The AASK Database V3.0: a database of human experience during aircraft evacuation incidents

This paper describes recent developments with the Aircraft Accident Statistics and Knowledge (AASK) database. The AASK database is a repository of survivor accounts from aviation accidents developed by the Fire Safety Engineering Group of the University of Greenwich with support from the UK CAA. Its main purpose is to store observational and anecdotal data from the actual interviews of the occupants involved in aircraft accidents. Access to the latest version of the database (AASK V3.0) is available over the Internet. AASK consists of information derived from both passenger and cabin crew interviews, information concerning fatalities and basic accident details. Also provided with AASK is the Seat Plan Viewer that graphically displays the starting locations of all the passengers - both survivors and fatalities - as well as the exits used by the survivors. Data entered into the AASK database is extracted from the transcripts supplied by the National Transportation Safety Board in the US and the Air Accident Investigation Branch in the UK. The quality and quantity of the data was very variable ranging from short summary reports of the accidents to boxes of individual accounts from passengers, crew and investigators. Data imported into AASK V3.0 includes information from 55 accidents and individual accounts from 1295 passengers and 110 crew.

The collection of pre-evacuation times from evacuation trials involving a Hospital Outpatient area and a University Library facility

This paper presents data relating to occupant pre-evacuation times from a university and a hospital outpatient facility. Although the two structures are entirely different they do employ relatively similar procedures: members of staff sweeping areas of the structure to encourage individuals to evacuate. However, the manner in which the dependent population reacts to these procedures is quite different. In the hospital case the patients only evacuated once a member of the nursing staff had instructed them to do so while in the university evacuation the students were less dependent upon the actions of the staff with over 50% of them evacuating with no prior prompting. Although this data may be useful in a variety of areas, it was collected primarily for use within evacuation models.

A generic time and space accurate numerical approach to closely coupled fluid-structure interaction problems

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, flow in elastic pipes and blood vessels and extrusion of metals through dies. However a comprehensive computational model of these multi-physics phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply even to the extent in metal forming, for example, that the deformation of the die is totally ignored. More recently, strategies for solving the full coupling between the fluid and soild mechanics behaviour have developed. Conventionally, the computational modelling of fluid structure interaction is problematical since computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. In the past the concurrent, but rather disparate, development paths for the finite element and finite volume methods have resulted in numerical software tools for CFD and CSM that are different in almost every respect. Hence, progress is frustrated in modelling the emerging multi-physics problem of fluid structure interaction in a consistent manner. Unless the fluid-structure coupling is either one way, very weak or both, transferring and filtering data from one mesh and solution procedure to another may lead to significant problems in computational convergence. Using a novel three phase technique the full interaction between the fluid and the dynamic structural response are represented. The procedure is demonstrated on some challenging applications in complex three dimensional geometries involving aircraft flutter, metal forming and blood flow in arteries.

Multi-physics-multi-scale simulation and optimisation: the next generation

FEA and CFD analysis is becoming ever more complex with an emerging demand for simulation software technologies that can address ranges of problems that involve combinations of interactions amongst varying physical phenomena over a variety of time and length scales. Computation modelling of such problems requires software technologies that enable the representation of these complex suites of 'physical' interactions. This functionality requires the structuring of simulation modules for specific physical phemonmena so that the coupling can be effectiely represented. These 'multi-physics' and 'multi-scale' computations are very compute intensive and so the simulation software must operate effectively in parallel if it is to be used in this context. Of course the objective of 'multi-physics' and 'multi-scale' simulation is the optimal design of engineered systems so optimistation is an important feature of such classes of simulation. In this presentation, a multi-disciplinary approach to simulation based optimisation is described with some key examples of application to challenging engineering problems.

Application of CFD in vacuum dezincing process

Most lead bullion is refined by pyrometallurgical methods - this involves a serics of processes that remove the antimony (softening) silver (Parkes process), zinc (vacuum dezincing) and if need be, bismuth (Betterton-Kroll process). The first step, softening, removes the antimony, arsenic and tin by air oxidation in a furnace or by the Harris process. Next, in the Parkes process, zinc is added to the melt to remove the silver and gold. Insoluble zinc, silver and gold compounds are skimmed off from the melt surface. Excess zinc added during desilvering is removed from lead bullion using one of ghree methods: * Vacuum dezincing; * Chlorine dezincing; or * Harris dezincing. The present study concentrates on the Vacuum dezincing process for lead refining. The main aims of the research are to develop mathematical model(s), using Computational Fluid Dyanmics (CFD) a Surface Averaged Model (SAM), to predict the process behaviour under various operating conditions, thus providing detailed information of the process - insight into its reaction to changes of key operating parameters. Finally, the model will be used to optimise the process in terms of initial feed concentration, temperature, vacuum height cooling rate, etc.

A multilevel algorithm for force-directed graph-drawing

We describe a heuristic method for drawing graphs which uses a multilevel framework combined with a force-directed placement algorithm. The multilevel technique matches and coalesces pairs of adjacent vertices to define a new graph and is repeated recursively to create a hierarchy of increasingly coarse graphs, G0, G1, …, GL. The coarsest graph, GL, is then given an initial layout and the layout is refined and extended to all the graphs starting with the coarsest and ending with the original. At each successive change of level, l, the initial layout for Gl is taken from its coarser and smaller child graph, Gl+1, and refined using force-directed placement. In this way the multilevel framework both accelerates and appears to give a more global quality to the drawing. The algorithm can compute both 2 & 3 dimensional layouts and we demonstrate it on examples ranging in size from 10 to 225,000 vertices. It is also very fast and can compute a 2D layout of a sparse graph in around 12 seconds for a 10,000 vertex graph to around 5-7 minutes for the largest graphs. This is an order of magnitude faster than recent implementations of force-directed placement algorithms.

Combined vertex-based - cell-centred finite volume method for flows in complex geometries

CFD modelling of 'real-life' processes often requires solutions in complex three dimensional geometries, which can often result in meshes where aspects of it are badly distorted. Cell-centred finite volume methods, typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes which feature cells with high non-orthogonal shapes. The vertex-based finite volume method handles distorted meshes with relative ease, but is computationally expensive. A combined vertex-based - cell-centred (VB-CC) technique, detailed in this paper, allows solutions on distorted meshes that defeat purely cell-centred physical models to be employed in the solution of other transported quantities. The VB-CC method is validated with benchmark solutions for thermally driven flow and turbulent flow. An early application of this hybrid technique is to three-dimensional flow over an aircraft wing, although it is planned to use it in a wide variety of processing applications in the future.

Computational modelling of bubbles, droplets and particles in metals reduction and refining

The CFD modelling of metals reduction processes particularly always seems to involve the interaction of liquid metals, a gas (often air) top space, liquid droplets in the top space and injection of both solid particles and gaseous bubbles into the bath. These phases all interact and exhange mass, momentum and energy. Often it is the extent to which these multi-phase phemomena can be effectively captured within the CFD model which determines whether or not a tool of genuine use to the target industry sector can constructed. In this paper we discuss these issues in the context of two problems - one involving the injection of sparging gases into a steel continuous caster and the other based on the development of a novel process for aluminium electrolysis.

Examining the effect of exit separation on aircraft evacuation performance using evacuation modelling techniques: "Is the 60 foot rule relevant?"

This paper examines the influence of exit separation, exit availability and seating configuration on aircraft evacuation efficiency and evacuation time. The purpose of this analysis is to explore how these parameters influence the 60 foot exit separation requirement found in aircraft certification rules. The analysis makes use of the airEXODUS evacuation model and is based on a typical wide-body aircraft cabin section involving two pairs of Type-A exits located at either end of the section with a maximum permissible loading of 220 passengers located between the exits. The analysis reveals that there is a complex relationship between exit separation and evacuation efficiency. Indeed, other factors such as exit flow rate and exit availability are shown to exert a strong influence on critical exit separations. A main finding of this work is that for the cabin section examined under certification conditions, exit separations up to 170 feet will result in approximately constant total evacuation times and average personal evacuation times. This practical exit separation threshold is decreased to 114 feet if another combination of exits is selected. While other factors must also be considered when determining maximum allowable exit separations, these results suggest it is not possible to mandate a maximum exit separation without taking into consideration exit type, exit availability and aircraft configuration. This has implications when determining maximum allowable exit separations for wide and narrow body aircraft. It is also relevant when considering the maximum allowable separation between different exit types on a given aircraft configuration.