The implementation of an aeronautical CFD flow code onto distributed memory parallel systems

The parallelization of an industrially important in-house computational fluid dynamics (CFD) code for calculating the airflow over complex aircraft configurations using the Euler or Navier–Stokes equations is presented. The code discussed is the flow solver module of the SAUNA CFD suite. This suite uses a novel grid system that may include block-structured hexahedral or pyramidal grids, unstructured tetrahedral grids or a hybrid combination of both. To assist in the rapid convergence to a solution, a number of convergence acceleration techniques are employed including implicit residual smoothing and a multigrid full approximation storage scheme (FAS). Key features of the parallelization approach are the use of domain decomposition and encapsulated message passing to enable the execution in parallel using a single programme multiple data (SPMD) paradigm. In the case where a hybrid grid is used, a unified grid partitioning scheme is employed to define the decomposition of the mesh. The parallel code has been tested using both structured and hybrid grids on a number of different distributed memory parallel systems and is now routinely used to perform industrial scale aeronautical simulations. Copyright © 2000 John Wiley & Sons, Ltd.

Flo/Stress: an integrated software module to predict stress in electronic products

Software technology that predicts stress in electronic systems and packages, developed as part of TCS Programme, is described. The software is closely integrated within a thermal design tool providing the ability to simulate the coupled effects of airflow, temperature and stress on product performance. This integrated approach to analysis will help decrease the number of design cycles.

Simulating the Swissair flight 111 in-flight fire using the CFD fire simulation software SMARTFIRE

At 8.18pm on 2 September 1998, Swissair Flight 111 (SR 111), took off from New York’s JFK airport bound for Geneva, Switzerland. Tragically, the MD-11 aircraft never arrived. According to the crash investigation report, published on 27 March 2003, electrical arcing in the ceiling void cabling was the most likely cause of the fire that brought down the aircraft. No one on board was aware of the disaster unfolding in the ceiling of the aircraft and, when a strange odour entered the cockpit, the pilots thought it was a problem with the air-conditioning system. Twenty minutes later, Swissair Flight 111 plunged into the Atlantic Ocean five nautical miles southwest of Peggy’s Cove, Nova Scotia, with the loss of all 229 lives on board. In this paper, the Computational Fluid Dynamics (CFD) analysis of the in-flight fire that brought down SR 111 is described. Reconstruction of the wreckage disclosed that the fire pattern was extensive and complex in nature. The fire damage created significant challenges to identify the origin of the fire and to appropriately explain the heat damage observed. The SMARTFIRE CFD software was used to predict the “possible” behaviour of airflow as well as the spread of fire and smoke within SR 111. The main aims of the CFD analysis were to develop a better understanding of the possible effects, or lack thereof, of numerous variables relating to the in-flight fire. Possible fire and smoke spread scenarios were studied to see what the associated outcomes would be. This assisted investigators at Transportation Safety Board (TSB) of Canada, Fire & Explosion Group in assessing fire dynamics for cause and origin determination.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Heat is extracted away from an electronic package by convection, conduction, and/or radiation. The amount of heat extracted by forced convection using air is highly dependent on the characteristics of the airflow around the package which includes its velocity and direction. Turbulence in the air is also important and is required to be modeled accurately in thermal design codes that use computational fluid dynamics (CFD). During air cooling the flow can be classified as laminar, transitional, or turbulent. In electronics systems, the flow around the packages is usually in the transition region, which lies between laminar and turbulent flow. This requires a low-Reynolds number numerical model to fully capture the impact of turbulence on the fluid flow calculations. This paper provides comparisons between a number of turbulence models with experimental data. These models included the distance from the nearest wall and the local velocity (LVEL), Wolfshtein, Norris and Reynolds, k-ε, k-ω, shear-stress transport (SST), and kε/kl models. Results show that in terms of the fluid flow calculations most of the models capture the difficult wake recirculation region behind the package reasonably well, although for packages whose heights cause a high degree of recirculation behind the package the SST model appears to struggle. The paper also demonstrates the sensitivity of the models to changes in the mesh density; this study is aimed specifically at thermal design engineers as mesh independent simulations are rarely conducted in an industrial environment.

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).

Multi-physics modelling for microelectronics and microsystems - current capabilities and future challenges

At present the vast majority of Computer-Aided- Engineering (CAE) analysis calculations for microelectronic and microsystems technologies are undertaken using software tools that focus on single aspects of the physics taking place. For example, the design engineer may use one code to predict the airflow and thermal behavior of an electronic package, then another code to predict the stress in solder joints, and then yet another code to predict electromagnetic radiation throughout the system. The reason for this focus of mesh-based codes on separate parts of the governing physics is essentially due to the numerical technologies used to solve the partial differential equations, combined with the subsequent heritage structure in the software codes. Using different software tools, that each requires model build and meshing, leads to a large investment in time, and hence cost, to undertake each of the simulations. During the last ten years there has been significant developments in the modelling community around multi- physics analysis. These developments are being followed by many of the code vendors who are now providing multi-physics capabilities in their software tools. This paper illustrates current capabilities of multi-physics technology and highlights some of the future challenges

Modelling and optimization of a pneumatic microfeeder system

This paper presents modelling and design optimization of a microfeeder which, as part of a microassembly system, is used for contactless object delivery. The microfeeder consists of an array of microactuators which are controlled by electrostatic actuation and used for maneuvering outcoming air jet for object hovering and delibery. The airflow behaviour in the microactuator is analysed by means of fluid mechanics and Computational Fluid Dynamics (CFD) simulation from three aspects, theoretical analysis, initial design assessment, and design modifications. The focus is put on the basic types of the microfeeder structure and the effects of structural details to the systematic performance. The structural pattern of the microactuator for forming airflow nozzle is identified and two design plans are proposed as basic structure patterns of pneumatic microactuators. The optimized design numerically shows the ability of delivering objects. This paper analyses the flow distribution pattern in microactuators and points out a way for effective design of pneumatic microfeeder systems. The optimization strategy provided by the present paper has close relevance to the design and manufacture of pneumatic microfeeder systems.

Observations and analysis of a 63-year-old reinforced concrete promenade railing exposed to the North Sea

Despite being exposed to the harsh sea-spray environment of the North Sea at Arbroath, Scotland, for over 63 years, many of the reinforced concrete precast beam elements of the 1.5 km long promenade railing are still in very good condition and show little evidence of reinforcement corrosion. In contrast, railing replacements constructed in about 1968 and in 1993 are almost all badly cracked as a result of extensive corrosion of the longitudinal reinforcement. This is despite the newer concrete appearing to be of better quality than the 1943 concrete. Statistics for maximum crack width for each of the three populations, based on measurements made in 2004 and in 2006, are presented. In situ and laboratory measurements show that the 1943 concrete appears to have high permeability but it also shows high electrical resistivity. Chloride penetration measurements show the 1943 and 1993 concretes to have similar chloride profiles and similar chloride concentrations at the reinforcement bars. This is inconsistent with the 1943 beams showing much less reinforcement corrosion than their later replacements and casts doubt on the conventional practice for durability design focusing on reducing concrete permeability through denser concretes or greater cover.

Condition assessment of a 63-year old reinforced concrete structure exposed to the North Sea

This paper describes the condition of a reinforced concrete balustrade consisting of some 1000 individual beam elements all exposed similarly to the hostile marine environment of the North Sea at Arbroath, Scotland since 1943. A comparison is made of the condition of the original construction with the condition of repairs carried out in 1968 and in 1993. It is shown that the 1943 construction shows very little corrosion-induced cracking and little rust staining even though it does not appear to be of high construction quality. Only a very low percentage of the balustrade beams have been replaced. In contrast the beam installed in 1968 and later in 1993 show very considerable and large concrete cracks directly attributable to the corrosion of the longitudinal reinforcement, even though the concrete is of a higher quality and density. A detailed condition survey and statistics of crack sizes are presented in the paper. It is found that the higher corrosion resistance of the 1943 concrete is generally consistent with the concrete electrical resistivity measurements but the degree of corrosion of the reinforcing bars is inconsistent with chloride penetration measurements. The results are compared with the very few observations available in the literature for ageing concrete structures in marine environments. The results cast doubt on the conventional wisdom that chloride content at the reinforcement correlates well with reinforcement corrosion.

Holocene climate change and past Irish societal response

The extent to which North Atlantic Holocene climatic perturbations influenced past human societies is an area of considerable uncertainty and fierce debate. Ireland is ideally placed to help resolve this issue, being occupied for over 9000 yr and located on the eastern Atlantic seaboard, a region dominated by westerly airflow. Irish bog and lake tree populations provide unambiguous evidence of major shifts in surface moisture through the Holocene similar to cycles recorded in the marine realm of the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. To test for human response to these cycles we summed the probabilities of 465 radiocarbon ages obtained from Irish archaeological contexts and observe enhanced archaeological visibility during periods of sustained wet conditions. These results suggest either increasing density of human populations in key, often defensive locations, and/or the development of subsistence strategies to overcome changing conditions, the latter recently proposed as a significant factor in avoiding societal collapse. Regardless, we demonstrate environmental change is a significantly more important factor in influencing human activity in the landscape than has hitherto been acknowledged.

LPCVD of tungsten by selective deposition on silicon

Beta-phase W, selectively grown at 440C had resistivity 20 micro-ohm cm and maximum layer thickness 100nm. Hydrogen passivation proved essential in this process. Higher deposition temperatures resulted in increased layer thickness but incorporated WSi2 and alpha- phase W. Self limiting W grown on polycrystalline and heavily doped silicon yielded reduced thickness. Boron is involved in the WF6 reduction reaction but phosphorus is not and becomes incorporated in the W layer. The paper establishes an optimised and novel CVD process suited to IC contact technology. A funded technology transfer contract with National Semiconductor Greenock (M Fallon) resulted from this work.

Non-linear conductance of disordered quantum wires

The self-consistent electron potential in a current-carrying disordered quantum wire is spatially inhomogeneous due to the formation of resistivity dipoles across scattering centres. In this paper it is argued that these inhomogeneities in the potential result in a suppression of the differential conductance of such a wire at finite applied voltage. A semi-classical argument allows this suppression, quadratic in the voltage, to be related directly to the amount of intrinsic defect scattering in the wire. This result is then tested against numerical calculations.

Geophysical imaging of stimulated microbial biomineralization

Understanding how microorganisms influence the physical and chemical properties of the subsurface is hindered by our inability to observe microbial dynamics in real time and with high spatial resolution. Here, we investigate the use of noninvasive geophysical methods to monitor biomineralization at the laboratory scale during stimulated sulfate reduction under dynamic flow conditions. Alterations in sediment characteristics resulting from microbe-mediated sulfide mineral precipitation were concomitant with changes in complex resistivity and acoustic wave propagation signatures. The sequestration of zinc and iron in insoluble sulfides led to alterations in the ability of the pore fluid to conduct electrical charge and of the saturated sediments to dissipate acoustic energy. These changes resulted directly from the nucleation, growth, and development of nanoparticulate precipitates along grain surfaces and within the pore space. Scanning and transmission electron microscopy (SEM and TEM) confirmed the sulfides to be associated with cell surfaces, with precipitates ranging from aggregates of individual 3-5 nm nanocrystals to larger assemblages of up to 10-20 m in diameter. Anomalies in the geophysical data reflected the distribution of mineral precipitates and biomass over space and time, with temporal variations in the signals corresponding to changes in the aggregation state of the nanocrystalline sulfides. These results suggest the potential for using geophysical techniques to image certain subsurface biogeochemical processes, such as those accompanying the bioremediation of metal-contaminated aquifers.

Non-linear electronic transport in Pt nanowires deposited by focused ion beam

Electrical transport and structural properties of platinum nanowires, deposited using the focussed ion beam method have been investigated. Energy dispersive X-ray spectroscopy reveals metal-rich grains (atomic composition 31% Pt and 50% Ga) in a largely non-metallic matrix of C, O and Si. Resistivity measurements (15-300 K) reveal a negative temperature coefficient with the room-temperature resistivity 80-300 times higher than that of bulk Pt. Temperature dependent current-voltage characteristics exhibit non-linear behaviour in the entire range investigated. The conductance spectra indicate increasing non-linearity with decreasing temperature, reaching 4% at 15 K. The observed electrical behaviour is explained in terms of a model for inter-grain tunnelling in disordered media, a mechanism that is consistent with the strongly disordered nature of the nanowires observed in the structure and composition analysis.

An efficient method of radiation attenuation through aircraft cooling vents

Some critical avionic systems require cooling air via vents on the side of the aircraft, thus creating leakage points for high-intensity electromagnetic radiation. This paper presents a novel application of high-intensity radiated field (HIRF) shielding using a rectangular waveguide array, while maintaining cooling airflow requirements. Signal attenuation versus frequency and depth of the array has been calculated using closed-form equations. The simulation and measurement results are in good agreement with the calculated values. (C) 2004 Wiley Periodicals, Inc.