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.

The dependence of power consumption of pneumatic conveying systems upon bulk material properties, pipeline bore and routing, and mode of flow

In this paper the dependence of the power consumption of pneumatic conveyors upon conveyed materials, pipeline route and bore, and mode of flow has been examined. The findings are that, with different materials and modes of flow, not only is the amount of power consumed very different but it varies in different ways with pipe bore and routing. Additionally it has been found that, for any given conveying system, the choice of air mover also has a strong influence on the power requirement.

Optimizing the operation of straight-grate iron-ore pellet induration systems using process models

Mathematical models of straight-grate pellet induration processes have been developed and carefully validated by a number of workers over the past two decades. However, the subsequent exploitation of these models in process optimization is less clear, but obviously requires a sound understanding of how the key factors control the operation. In this article, we show how a thermokinetic model of pellet induration, validated against operating data from one of the Iron Ore Company of Canada (IOCC) lines in Canada, can be exploited in process optimization from the perspective of fuel efficiency, production rate, and product quality. Most existing processes are restricted in the options available for process optimization. Here, we review the role of each of the drying (D), preheating (PH), firing (F), after-firing (AF), and cooling (C) phases of the induration process. We then use the induration process model to evaluate whether the first drying zone is best to use on the up- or down-draft gas-flow stream, and we optimize the on-gas temperature profile in the hood of the PH, F, and AF zones, to reduce the burner fuel by at least 10 pct over the long term. Finally, we consider how efficient and flexible the process could be if some of the structural constraints were removed (i.e., addressed at the design stage). The analysis suggests it should be possible to reduce the burner fuel lead by 35 pct, easily increase production by 5+ pct, and improve pellet quality.

Computational modelling of multi-physics processes on high performance parallel computer systems

The demands of the process of engineering design, particularly for structural integrity, have exploited computational modelling techniques and software tools for decades. Frequently, the shape of structural components or assemblies is determined to optimise the flow distribution or heat transfer characteristics, and to ensure that the structural performance in service is adequate. From the perspective of computational modelling these activities are typically separated into: • fluid flow and the associated heat transfer analysis (possibly with chemical reactions), based upon Computational Fluid Dynamics (CFD) technology • structural analysis again possibly with heat transfer, based upon finite element analysis (FEA) techniques.

Nanoscopic modelling of the adhesion, indentation and fracture characteristics of metallic systems via molecular dynamics simulations

Large-scale molecular dynamics simulations have been performed on canonical ensembles to model the adhesion and indentation characteristics of 3-D metallic nano-scale junctions in tip-substrate geometries, and the crack propagation in 2-D metallic lattices. It is shown that irreversible flows in nano-volumes of materials control the behaviour of the 3-D nano-contacts, and that local diffusional flow constitutes the atomistic mechanism underlying these plastic flows. These simulations show that the force of adhesion in metallic nano-contacts is reduced when adsorbate monolayers are present at the metal—metal junctions. Our results are in agreement with the conclusions of very accurate point-contact experiments carried out in this field. Our fracture simulations reveal that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures, the nucleation of dislocations is shown to cause a brittle-to-ductile transition. Limiting crack propagation velocities are computed for different strain rates and a dynamic instability is shown to control the crack movement beyond this limiting velocity, in line with the recent experimental results.

Turbulence modelling and its effects on electronic systems

This paper presents simulated computational fluid dynamics (CFD) results for comparison against experimental data. The performance of four turbulence models has been assessed for electronic application areas considering both fluid flow and heat transfer phenomenon. CFD is vast becoming a powerful and almost essential tool for design, development and optimization in engineering problems. However turbulence models remain to be the key problem issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the performance of the turbulence model employed together with the wall functions implemented. To be able to resolve the abrupt changes in the turbulent energy and other parameters near the wall a particularly fine mesh is necessary which unfortunately increases the computer storage capacity requirements. The objective of turbulence modelling is to enhance computational procdures of sufficient acccuracy and generality for engineers to anticipate the Reynolds stresses and the scalar transport terms.

An integrated approach to flow, thermal and mechanical modeling of electronics devices

The future success of many electronics companies will depend to a large extent on their ability to initiate techniques that bring schedules, performance, tests, support, production, life-cycle-costs, reliability prediction and quality control into the earliest stages of the product creation process. Earlier papers have discussed the benefits of an integrated analysis environment for system-level thermal, stress and EMC prediction. This paper focuses on developments made to the stress analysis module and presents results obtained for an SMT resistor. Lifetime predictions are made using the Coffin-Manson equation. Comparison with the creep strain energy based models of Darveaux (1997) shows the shear strain based method to underestimate the solder joint life. Conclusions are also made about the capabilities of both approaches to predict the qualitative and quantitative impact of design changes.

Turbulence modelling for thermal management of electronic systems

The electronics industry is developing rapidly together with the increasingly complex problem of microelectronic equipment cooling. It has now become necessary for thermal design engineers to consider the problem of equipment cooling at some level. The use of Computational Fluid Dynamics (CFD) for such investigations is fast becoming a powerful and almost essential tool for the design, development and optimisation of engineering applications. However turbulence models remain a key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt fluctuations experienced by the turbulent energy and other parameters located at near wall regions and shear layers a particularly fine computational mesh is necessary which inevitably increases the computer storage and run-time requirements. This paper will discuss results from an investigation into the accuract of currently used turbulence models. Also a newly formulated transitional hybrid turbulence model will be introduced with comparisonsaagainst experimental data.

No-flow underfill flip chip assembly–an experimental and modeling analysis

In the flip-chip assembly process, no-flow underfill materials have a particular advantage over traditional underfill: the application and curing of the former can be undertaken before and during the reflow process. This advantage can be exploited to increase the flip-chip manufacturing throughput. However, adopting a no-flow underfill process may introduce reliability issues such as underfill entrapment, delamination at interfaces between underfill and other materials, and lower solder joint fatigue life. This paper presents an analysis on the assembly and the reliability of flip-chips with no-flow underfill. The methodology adopted in the work is a combination of experimental and computer-modeling methods. Two types of no-flow underfill materials have been used for the flip chips. The samples have been inspected with X-ray and scanning acoustic microscope inspection systems to find voids and other defects. Eleven samples for each type of underfill material have been subjected to thermal shock test and the number of cycles to failure for these flip chips have been found. In the computer modeling part of the work, a comprehensive parametric study has provided details on the relationship between the material properties and reliability, and on how underfill entrapment may affect the thermal–mechanical fatigue life of flip chips with no-flow underfill.

Accuracy of turbulence models and CFD for thermal characterisation of electronic systems

A major percentage of the heat emitted from electronic packages can be extracted by air cooling whether by means of natural or forced convection. This flow of air throughout an electronic system and the heat extracted is highly dependable on the nature of turbulence present in the flow field. This paper will discuss results from an investigation into the accuracy of turbulence models to predict air cooling for electronic packages and systems.

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.

Effect of fluid dynamics and device mechanism on biofluid behaviour in microchannel systems: modelling biofluids in a microchannel biochip separator

Biofluid behaviour in microchannel systems is investigated in this paper through the modelling of a microfluidic biochip developed for the separation of blood plasma. Based on particular assumptions, the effects of some mechanical features of the microchannels on behaviour of the biofluid are explored. These include microchannel, constriction, bending channel, bifurcation as well as channel length ratio between the main and side channels. The key characteristics and effects of the microfluidic dynamics are discussed in terms of separation efficiency of the red blood cells with respect to the rest of the medium. The effects include the Fahraeus and Fahraeus-Lindqvist effects, the Zweifach-Fung bifurcation law, the cell-free layer phenomenon. The characteristics of the microfluid dynamics include the properties of the laminar flow as well as particle lateral or spinning trajectories. In this paper the fluid is modelled as a single-phase flow assuming either Newtonian or Non-Newtonian behaviours to investigate the effect of the viscosity on flow and separation efficiency. It is found that, for a flow rate controlled Newtonian flow system, viscosity and outlet pressure have little effect on velocity distribution. When the fluid is assumed to be Non-Newtonian more fluid is separated than observed in the Newtonian case, leading to reduction of the flow rate ratio between the main and side channels as well as the system pressure as a whole.

Sediment transfer and accumulation in two contrasting salt marsh/mudflat systems: the Seine estuary (France) and the Medway estuary (UK)

Understanding the dynamics of fine sediment transport across the upper intertidal zone is critical in managing the erosion and accretion of intertidal areas, and in managed realignment/estuarine habitat recreation strategies. This paper examines the transfer of sediments between salt marsh and mudflat environments in two contrasting macrotidal estuaries: the Seine (France) and the Medway (UK), using data collected during two joint field seasons undertaken by the Anglo-French RIMEW project (Rives-Manche Estuary Watch). High-resolution ADCP, Altimeter, OBS and ASM measurements from mudflat and marsh surface environments have been combined with sediment trap data to examine short-term sediment transport processes under spring tide and storm flow conditions. In addition, the longer-term accumulation of sediment in each salt marsh system has been examined via radiometric dating of sediment cores. In the Seine, rapid sediment accumulation and expansion of salt marsh areas, and subsequent loss of open intertidal mudflats, is a major problem, and the data collected here indicate a distinct net landward flux of sediments into the marsh interior. Suspended sediment fluxes are much higher than in the Medway estuary (averaging 0.09 g/m(3)/s), and vertical accumulation rates at the salt marsh/mudflat boundary exceed 3 cm/y. Suspended sediment data collected during storm surge conditions indicate that significant in-wash of fine sediments into the marsh interior can occur during (and following) these high-magnitude events. In contrast to the Seine, the Medway is undergoing erosion and general loss of salt marsh areas. Suspended sediment fluxes are of the order of 0.03 g/m(3)/s, and the marsh system here has much lower rates of vertical accretion (sediment accumulation rates are ca. 4 mm/y). Current velocity data for the Medway site indicate higher velocities on the ebb tide than occur on the flood tide, which may be sufficient to remobilise sediments deposited on the previous tide and so force net removal of material from the marsh.

A semi-Lagrangian finite volume method for solving viscoelastic flow problems

Recently, there has been considerable interest in solving viscoelastic problems in 3D particularly with the improvement in modern computing power. In many applications the emphasis has been on economical algorithms which can cope with the extra complexity that the third dimension brings. Storage and computer time are of the essence. The advantage of the finite volume formulation is that a large amount of memory space is not required. Iterative methods rather than direct methods can be used to solve the resulting linear systems efficiently.