Load balancing for aspect ratio in adaptive finite element simulations

In parallel adaptive finite element simulations the work load on the individual processors may change frequently. To (re)distribute the load evenly over the processors a load balancing heuristic is needed. Common strategies try to minimise subdomain dependencies by optimising the cutsize of the partitioning. However for certain solvers cutsize only plays a minor role, and their convergence is highly dependent on the subdomain shapes. Degenerated subdomain shapes cause them to need significantly more iterations to converge. In this work a new parallel load balancing strategy is introduced which directly addresses the problem of generating and conserving reasonably good subdomain shapes in a dynamically changing Finite Element Simulation. Geometric data is used to formulate several cost functions to rate elements in terms of their suitability to be migrated. The well known diffusive method which calculates the necessary load flow is enhanced by weighting the subdomain edges with the help of these cost functions. The proposed methods have been tested and results are presented.

SMARTFIRE: an integrated computational fluid dynamics code and expert system for fire field modelling

SMARTFIRE, an open architecture integrated CFD code and knowledge based system attempts to make fire field modeling accessible to non-experts in Computational Fluid Dynamics (CFD) such as fire fighters, architects and fire safety engineers. This is achieved by embedding expert knowledge into CFD software. This enables the 'black-art' associated with the CFD analysis such as selection of solvers, relaxation parameters, convergence criteria, time steps, grid and boundary condition specification to be guided by expert advice from the software. The user is however given the option of overriding these decisions, thus retaining ultimate control. SMARTFIRE also makes use of recent developments in CFD technology such as unstructured meshes and group solvers in order to make the CFD analysis more efficient. This paper describes the incorporation within SMARTFIRE of the expert fire modeling knowledge required for automatic problem setup and mesh generation as well as the concept and use of group solvers for automatic and manual dynamic control of the CFD code.

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.

Unstructured staggered mesh methods for fluid flow, heat transfer and phase change

Unstructured grid meshes used in most commercial CFD codes inevitably adopt collocated variable solution schemes. These schemes have several shortcomings, mainly due to the interpolation of the pressure gradient, that lead to slow convergence. In this publication we show how it is possible to use a much more stable staggered mesh arrangement in an unstructured code. Several alternative groupings of variables are investigated in a search for the optimum scheme.

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

A multi-phase framework is typically required for the CFD modelling of metals reduction processes. Such processes typically 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. The exchange of mass, momentum and energy between the phases is fundamental to these processes. Multi-phase algorithms are complex and can be unreliable in terms of either or both convergence behaviour or in the extent to which the physics is captured. In this contribution, we discuss these multi-phase flow issues and describe an example of each of the main “single phase” approaches to modelling this class of problems (i.e., Eulerian–Lagrangian and Eulerian–Eulerian). Their utility is illustrated in the context of two problems – one involving the injection of sparging gases into a steel continuous slab caster and the other based on the development of a novel process for aluminium electrolysis. In the steel caster, the coupling of the Lagrangian tracking of the gas phase with the continuum enables the simulation of the transient motion of the metal–flux interface. The model of the electrolysis process employs a novel method for the calculation of slip velocities of oxygen bubbles, resulting from the dissolution of alumina, which allows the efficiency of the process to be predicted.

Scalable and efficient graph colouring in 3 dimensions using emergence engineering principles

This paper describes ways in which emergence engineering principles can be applied to the development of distributed applications. A distributed solution to the graph-colouring problem is used as a vehicle to illustrate some novel techniques. Each node acts autonomously to colour itself based only on its local view of its neighbourhood, and following a simple set of carefully tuned rules. Randomness breaks symmetry and thus enhances stability. The algorithm has been developed to enable self-configuration in wireless sensor networks, and to reflect real-world configurations the algorithm operates with 3 dimensional topologies (reflecting the propagation of radio waves and the placement of sensors in buildings, bridge structures etc.). The algorithm’s performance is evaluated and results presented. It is shown to be simultaneously highly stable and scalable whilst achieving low convergence times. The use of eavesdropping gives rise to low interaction complexity and high efficiency in terms of the communication overheads.

Numerical and experimental studies of a substructural identification strategy

For structural health monitoring it is impractical to identify a large structure with complete measurement due to limited number of sensors and difficulty in field instrumentation. Furthermore, it is not desirable to identify a large number of unknown parameters in a full system because of numerical difficulty in convergence. A novel substructural strategy was presented for identification of stiffness matrices and damage assessment with incomplete measurement. The substructural approach was employed to identify large systems in a divide-and-conquer manner. In addition, the concept of model condensation was invoked to avoid the need for complete measurement, and the recovery process to obtain the full set of parameters was formulated. The efficiency of the proposed method is demonstrated numerically through multi-storey shear buildings subjected to random force. A fairly large structural system with 50 DOFs was identified with good results, taking into consideration the effects of noisy signals and the limited number of sensors. Two variations of the method were applied, depending on whether the sensor could be repositioned. The proposed strategy was further substantiated experimentally using an eight-storey steel plane frame model subjected to shaker and impulse hammer excitations. Both numerical and experimental results have shown that the proposed substructural strategy gave reasonably accurate identification in terms of locating and quantifying structural damage.