A finite volume unstructured mesh approach to dynamic fluid–structure interaction: an assessment of the challenge of predicting the onset of flutter

Computational modelling of dynamic fluid–structure interaction (DFSI) is a considerable challenge. Our approach to this class of problems involves the use of a single software framework for all the phenomena involved, employing finite volume methods on unstructured meshes in three dimensions. This method enables time and space accurate calculations in a consistent manner. One key application of DFSI simulation is the analysis of the onset of flutter in aircraft wings, where the work of Yates et al. [Measured and Calculated Subsonic and Transonic Flutter Characteristics of a 45° degree Sweptback Wing Planform in Air and Freon-12 in the Langley Transonic Dynamic Tunnel. NASA Technical Note D-1616, 1963] on the AGARD 445.6 wing planform still provides the most comprehensive benchmark data available. This paper presents the results of a significant effort to model the onset of flutter for the AGARD 445.6 wing planform geometry. A series of key issues needs to be addressed for this computational approach. • The advantage of using a single mesh, in order to eliminate numerical problems when applying boundary conditions at the fluid-structure interface, is counteracted by the challenge of generating a suitably high quality mesh in both the fluid and structural domains. • The computational effort for this DFSI procedure, in terms of run time and memory requirements, is very significant. Practical simulations require even finer meshes and shorter time steps, requiring parallel implementation for operation on large, high performance parallel systems. • The consistency and completeness of the AGARD data in the public domain is inadequate for use in the validation of DFSI codes when predicting the onset of flutter.

[Book review] Owen Gingerich, The Book Nobody Read: In Pursuit of the Revolutions of Nicolaus Copernicus, Heinemann, 2004

Tony Mann reviews: Owen Gingerich, The Book Nobody Read: In Pursuit of the Revolutions of Nicolaus Copernicus, Heinemann, 2004, 0-434-01315-3, £12.99.

CFD fire simulation of the Swissair flight 111 in-flight fire - part II: fire spread analysis

In 1998, Swissair Flight I I I (SR111) developed an in-flight fire shortly after take-off which resulted in the loss of the aircraft, a McDonnell Douglas MD-I 1, and all passengers and crew. The Transportation Safety Board (TSB) of Canada, Fire and Explosion Group launched a four year investigation into the incident in an attempt to understand the cause and subsequent mechanisms which lead to the rapid spread of the in-flight fire. As part of this investigation, the SMARTFIRE Computational Fluid Dynamics (CFD) software was used to predict the 'possible' development of the fire and associated smoke movement. In this paper the CFD fire simulations are presented and model predictions compared with key findings from the investigation. The model predictions are shown to be consistent with a number of the investigation findings associated with the early stages of the fire development. The analysis makes use of simulated pre-fire airflow conditions within the MD-11 cockpit and above ceiling region presented in an earlier publication (Part 1) which was published in The Aeronautical Journal in January 2006(4).

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Thermo-mechanical analysis of the packaging process for micro-electronic displays

Micro-electronic displays are sensitive devices and its performance is easily affected by external environmental factors. To enable the display to perform in extreme conditions, the device must be structurally strengthened, the effects of this packaging process was investigated. A thermo-mechanical finite element analysis was used to discover potential problems in the packaging process and to improve the overall design of the device. The main concern from the analysis predicted that displacement of the borosilicate glass and the Y stress of the adhesive are important. Using this information a design which reduced the variation of displacement and kept the stress to a minimum was suggested

An analysis of the passenger to cabin crew ratio and exit reliability based on past survivable aviation accidents

A hotly debated issue in the area of aviation safety is the number of cabin crew members required to evacuate an aircraft in the event of an emergency. Most countries regulate the minimum number required for the safe operation of an aircraft, but these rulings are based on little if any scientific evidence. Another issue of concern is the failure rate of exits and slides. This paper examines these issues using the latest version of Aircraft Accident Statistics and Knowledge database AASK V4.0, which contains information from 105 survivable crashes and more than 2,000 survivors, including accounts from 155 cabin crew members.

Computational modeling of mold filling and related free-surface flows in shape casting: an overview of the challenges involved

Accurate representation of the coupled effects between turbulent fluid flow with a free surface, heat transfer, solidification, and mold deformation has been shown to be necessary for the realistic prediction of several defects in castings and also for determining the final crystalline structure. A core component of the computational modeling of casting processes involves mold filling, which is the most computationally intensive aspect of casting simulation at the continuum level. Considering the complex geometries involved in shape casting, the evolution of the free surface, gas entrapment, and the entrainment of oxide layers into the casting make this a very challenging task in every respect. Despite well over 30 years of effort in developing algorithms, this is by no means a closed subject. In this article, we will review the full range of computational methods used, from unstructured finite-element (FE) and finite-volume (FV) methods through fully structured and block-structured approaches utilizing the cut-cell family of techniques to capture the geometric complexity inherent in shape casting. This discussion will include the challenges of generating rapid solutions on high-performance parallel cluster technology and how mold filling links in with the full spectrum of physics involved in shape casting. Finally, some indications as to novel techniques emerging now that can address genuinely arbitrarily complex geometries are briefly outlined and their advantages and disadvantages are discussed.

Reliability study of the electroless Ni–P layer against solder alloy

Copper (Cu) has been widely used in the under bump metallurgy of chip and substrate metallization for chip packaging. However, due to the rapid formation of Cu–Sn intermetallic compound (IMC) at the tin-based solder/Cu interface during solder reaction, the reliability of this type of solder joint is a serious concern. In this work, electroless nickel–phosphorous (Ni–P) layer was deposited on the Cu pad of the flexible substrate as a diffusion barrier between Cu and the solder materials. The deposition was carried out in a commercial acidic sodium hypophosphite bath at 85 °C for different pH values. It was found that for the same deposition time period, higher pH bath composition (mild acidic) yields thicker Ni–P layer with lower phosphorous content. Solder balls having composition 62%Sn–36%Pb–2%Ag were reflowed at 240 °C for 1 to 180 min on three types of electroless Ni–P layers deposited at the pH value of 4, 4.8 and 6, respectively. Thermal stability of the electroless Ni–P barrier layer against the Sn–36%Pb–2%Ag solder reflowed for different time periods was examined by scanning electron microscopy equipped with energy dispersed X-ray. Solder ball shear test was performed in order to find out the relationship between the mechanical strength of solder joints and the characteristics of the electroless Ni–P layer deposited. The layer deposited in the pH 4 acidic bath showed the weak barrier against reflow soldering whereas layer deposited in pH 6 acidic bath showed better barrier against reflow soldering. Mechanical strength of the joints were deteriorated quickly in the layer deposited at pH 4 acidic bath, which was found to be thin and has a high phosphorous content. From the cross-sectional studies and fracture surface analyses, it was found that the appearance of the dark crystalline phosphorous-rich Ni layer weakened the interface and hence lower solder ball shear strength. Ni–Sn IMC formed at the interfaces was found to be more stable at the low phosphorous content (∼14 at.%) layer. Electroless Ni–P deposited at mild acidic bath resulting phosphorous content of around 14 at.% is suggested as the best barrier layer for Sn–36%Pb–2%Ag solder.

Computer modelling analysis of the globtop's effects on aluminium wirebond reliability

In power electronics modules, heavy aluminium wires, i.e. wire diameters greater than 100 microns, are bonded to the active semiconductor devices and conductor metallization to form electric circuits of the power electronic module. Due to the high currents that may flow through these wires, a single connection usually contains several wires and thus, a large number of wires are used in a power electronics module. Under normal operation or test condition, a significant amount of stresses and strains induced in the wire and bonding interfaces, resulting in failure over time. In this paper, computer modelling techniques are used to analyse the effect of globtop design on the reliability of aluminium wirebonds under cyclic thermal-mechanical loading conditions. The results will show the sensitivity of the reliability of the wirebonds to the changes in the geometry and the material properties of the wirebond globtop.

Three-dimensional simulation of the influence of convection on dendritic solidification

Bulk and interdendritic flow during solidification alters the microstructure development, potentially leading to the formation of defects. In this paper, a 3D numerical model is presented for the simulation of dendritic growth in the presence of fluid flow in both liquid and semi-solid zones during solidification. The dendritic growth was solved by the combination of a stochastic nucleation approach with a finite difference solution of the solute diffusion equation and. a projection method solution of the Navier-Stokes equations. The technique was applied first to simulate the growth of a single dendrite in 2D and 3D in an isothermal environment with forced fluid flow. Significant differences were found in the evolution of dendritic morphology when comparing the 2D and 3D results. In 3D the upstream arm has a faster growth velocity due to easier flow around the perpendicular arms. This also promotes secondary arm formation on the upstream arm. The effect of fluid flow on columnar dendritic growth and micro-segregation in constrained solidification conditions is then simulated. For constrained growth, 2D simulations lead to even greater inaccuracies as compared to 3D.

Impact of the vertical potshell deformation on the MHD cell stability behavior of a 500 kA aluminium electrolysis cell

In previous publications [1,2], it was rationalized that a large vertical potshell deformation may have a negative impact on the operations of very high amperage cells. The MHD-Valdis non-linear Magneto-Hydro-Dynamic model was therefore extended to take into account the displacement of the potshell. The MHD cell stability behavior of a 500 kA cell with a 17.3 meters long potshell was then studied.

Approximating the evacuation of the World Trade Center North Tower using computer simulation

This article concerns an investigation of the full scale evacuation of a building with a configuration similar to that of the World Trade Center (WTC) North Tower using computer simulation. A range of evacuation scenarios is explored in order to better understand the evacuation of the WTC on 11 September 2001. The analysis makes use of response time data derived from a study of published WTC survivor accounts. Geometric details of the building are obtained from architects' plans while the total building population used in the scenarios is based on estimates produced by the National Institute of Standards and Technology formal investigation into the evacuation. This paper attempts to approximate the events of 11 September 2001 and pursue several `what if' questions concerning the evacuation. In particular, the study explores the likely outcome had a single staircase survived intact from top to bottom. More generally, this paper explores issues associated with the practical limits of building size that can be expected to be efficiently evacuated using stairs alone.

Time-dependent modelling and experimental validation of the metal/flux interface in a continuous casting mould

A finite volume computer model of the continuous casting process for steel flat products has been developed. In this first stage, the model concentrates on the hydrodynamic aspects of the process and in particular the dynamic behavior of the metal/slag interface. The model was validated against experimental measurements obtained in a water model apparatus.

Projecting the prison: the depiction of the uncanny in The Shawshank Redemption

The space of the prison is no longer on the margins in relation to societal `centres', but instead acts as an adjunct to the urban environment. With the disappearance of the Gothic prison from the archi-texture of contemporary cities, the meaning conveyed by its façade has lost much of its potency. It is now contemporary prison drama, as opposed to the physical façade, that represents the interface between the public and the prison. This article explores a dramatic representation of the prison (The Shawshank Redemption) through the lens of Freud's (1919/1955) notion of the uncanny and Bachelard's (1958/1994) poetics of domestic space. Incarceration, as depicted in film and television, reinforces the `place myths' of the prison (Shields, 1991). Contemporary prison drama portrays the prison as a marginal space in much the way that the Gothic façades of the 19th-century prison projected a particular message. The prison, as depicted on screen, is a simulacrum. It is a facsimile of an architectural idea that only ever existed as a façade - a façade that occluded as much as it projected.