Simulations of Two-Dimensional Foam under Couette Shear

null Sponsorship: Financial support is acknowledged from the University of Wales Aberystwyth Senate Fund, the Ulysses France-Ireland Exchange Scheme and EPSRC (EP/D014956/1, EP/D048397/1, EP/D071127/1).

Observations of spatial flow patterns at the coral colony scale on a shallow reef flat

Although small-scale spatial flow variability can affect both larger-scale circulation patterns and biological processes on coral reefs, there are few direct measurements of spatial flow patterns across horizontal scales <100 m. Here flow patterns on a shallow reef flat were measured at scales from a single colony to several adjacent colonies using an array of acoustic Doppler velocimeters on a diver-operated traverse. We observed recirculation zones immediately behind colonies, reduced currents and elevated dissipation rates in turbulent wakes up to 2 colony diameters downstream and enhanced Reynolds stresses in shear layers around wake peripheries. Flow acceleration zones were observed above and between colonies. Coherent flow structures varied with incident flow speeds; recirculation zones were stronger and wakes were more turbulent in faster flows. Low-frequency (<0.03 Hz) flow variations, for which water excursions were large compared with the colony diameters (Keulegan-Carpenter number, KC >1), had similarspatial patterns to wakes, while higher-frequency variations (0.05-0.1 Hz, KC<1) had no observable spatial structure. On the reef flat, both drag and inertial forces exerted by coral colonies could have significant effects on flow, but within different frequency ranges; drag dominates for low-frequency flow variations and inertial forces dominate for higher frequency variations, including the wave band. Our scaling analyses suggest that spatial flow patterns at colony and patch scales could have important implications or both physical and biological processes at larger reef scales through their effects on forces exerted on the flow, turbulent mixing, and dispersion. © 2013. American Geophysical Union. All Rights Reserved.

Liquid metal turbulent flow dynamics in a cylindrical container with free surface: experiment and numerical analysis

The liquid metal flow in induction crucible models is known to be unstable, turbulent and difficult to predict in the regime of medium frequencies when the electromagnetic skin-layer is of considerable extent. We present long term turbulent flow measurements by a permanent magnet incorporated potential difference velocity probe in a cylindrical container filled with eutectic melt In-Ga-Sn. The parallel numerical simulation of the long time scale development of the turbulent average flow is presented. The numerical flow model uses an implicit pseudo-spectral code and k-w turbulence model, which was recently developed for the transitional flow modelling. The results compare reasonably to the experiment and demonstrate the time development of the turbulent flow field and the turbulence energy.

Liquid metal turbulent flow dynamics in a cylindrical container with free surface: experiment and numerical analysis

The liquid metal flow in inducation crucible models is known to be higly unstable and turbutlen in the regim e of medium frequecies when the elctronmagnetic skin-layer is of considerable extent. We present long term turbulent flow measurements by a permanent magnet incorporated potential difference veolocity probe in a cylindirical container filled with eutecti mlt In-Ga-SN. The parallel numerical simulation of the long time scale development of the turbulen average flow is presented. The numerical lfow model uses a pseud-spectral code and k-w turbulence model, which was recently developed for the transitional flow modelling. The result compare reasonably to the experiment and demonstrate the time development of the turbulent flow field.

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 and its impact on CFD predictions for cooling of electronic components

This paper will discuss Computational Fluid Dynamics (CFD) results from an investigation into the accuracy of several turbulence models to predict air cooling for electronic packages and systems. Also new transitional turbulence models will be proposed with emphasis on hybrid techniques that use the k-ε model at an appropriate distance away from the wall and suitable models, with wall functions, near wall regions. A major proportion of heat emitted from electronic packages can be extracted by air cooling. This flow of air throughout an electronic system and the heat extracted is highly dependent on the nature of turbulence present in the flow. The use of CFD for such investigations is fast becoming a powerful and almost essential tool for the design, development and optimization 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. The PHYSICA Finite Volume code was used for this investigation. With the exception of the k-ε and k-ω models which are available as standard within PHYSICA, all other turbulence models mentioned were implemented via the source code by the authors. The LVEL, LVEL CAP, Wolfshtein, k-ε, k-ω, SST and kε/kl models are described and compared with experimental data.

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.

Experience of using a CFD code for estimating the noise generated by gusts along the sun-roof of a car

The problem to be examined here is the fluctuating pressure distribution along the open cavity of the sun-roof at the top of a car compartment due to gusts passing over the sun-roof. The aim of this test is to investigate the capability of a typical commercial CFD package, PHOENICS, in recognising pressure fluctuations occurring in an important automotive industrial problem. In particular to examine the accuracy of transporting pulsatory gusts traveling along the main flow through the use of finite volume methods with higher order schemes in the numercial solutins of the unsteady compressible Navier-Stokes equations. The Helmholtz equation is used to solve the sound distribution inside the car compartment, resulting from the externally induced fluctuations.

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.

Effect of machine dependence on wall failure property measurements, Jenike versus an annular shear tester

Annular, ring or torsional shear testers are commonly used in bulk solids handling research for the purpose of powder characterisation or equipment design. This paper reports from a DEFRA sponsored project which aims to develop an industrial powder flow-ability tester, (based on the annular shear tester) that is economic to buy and quick and easy to use in trained but unskilled hands. This paper compares the wall failure loci measured with an annular shear cell with measurements obtained using the accepted standard wall friction tester, the Jenike shear cell. These wall failure loci have been measured for several bulk solids which range from fine cohesive powders to free-flowing granular materials, on a stainless steel 304 2B wall surface.

Modelling of the time-dependent flow behaviour of lead-free solder pastes used for flip-chip assembly applications

The market for solder paste materials in the electronic manufacturing and assembly sector is very large and consists of material and equipment suppliers and end users. These materials are used to bond electronic components (such as flip-chip, CSP and BGA) to printed circuit boards (PCB's) across a range of dimensions where the solder interconnects can be in the order of 0.05mm to 5mm in size. The non-Newtonian flow properties exhibited by solder pastes during its manufacture and printing/deposition phases have been of practical concern to surface mount engineers and researchers for many years. The printing of paste materials through very small-sized stencil apertures is known to lead to increased stencil clogging and incomplete transfer of paste to the substrate pads. At these very narrow aperture sizes the paste rheology and particle-wall interactions become crucial for consistent paste withdrawal. These non-Newtonian effects must be understood so that the new paste formulations can be optimised for consistent printing. The focus of the study reported in this paper is the characterisation of the rheological properties of solder pastes and flux mediums, and the evaluation of the effect of these properties on the pastes' printing performance at the flip-chip assembly application level. Solder pastes are known to exhibit a thixotropic behaviour, which is recognised by the decrease in apparent viscosity of paste material with time when subjected to a constant shear rate. The proper characterisation of this time-dependent theological behaviour of solder pastes is crucial for establishing the relationships between the pastes' structure and flow behaviour; and for correlating the physical parameters with paste printing performance. In this paper, we present a number of methods which have been developed for characterising the time-dependent and non-Newtonian rheological behaviour of solder pastes and flux mediums as a function of shear rates. We also present results of the study of the rheology of the solder pastes and flux mediums using the structural kinetic modelling approach, which postulates that the network structure of solder pastes breaks down irreversibly under shear, leading to time and shear dependent changes in the flow properties. Our results show that for the solder pastes used in the study, the rate and extent of thixotropy was generally found to increase with increasing shear rate. The technique demonstrated in this study has wide utility for R&D personnel involved in new paste formulation, for implementing quality control procedures used in solder paste manufacture and packaging; and for qualifying new flip-chip assembly lines

Examination of the flow rheological and textural properties of polymer gels composed of poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone): rheological and mathematical interpretation of textural parameters.

The purpose of this study was to mathematically characterize the effects of defined experimental parameters (probe speed and the ratio of the probe diameter to the diameter of sample container) on the textural/mechanical properties of model gel systems. In addition, this study examined the applicability of dimensional analysis for the rheological interpretation of textural data in terms of shear stress and rate of shear. Aqueous gels (pH 7) were prepared containing 15% w/w poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone) (PVP) (0, 3, 6, or 9% w/w). Texture profile analysis (TPA) was performed using a Stable Micro Systems texture analyzer (model TA-XT 2; Surrey, UK) in which an analytical probe was twice compressed into each formulation to a defined depth (15 mm) and at defined rates (1, 3, 5, 8, and 10 mm s-1), allowing a delay period (15 s) between the end of the first and beginning of the second compressions. Flow rheograms were performed using a Carri-Med CSL2-100 rheometer (TA Instruments, Surrey, UK) with parallel plate geometry under controlled shearing stresses at 20.0°?±?0.1°C. All formulations exhibited pseudoplastic flow with no thixotropy. Increasing concentrations of PVP significantly increased formulation hardness, compressibility, adhesiveness, and consistency. Increased hardness, compressibility, and consistency were ascribed to enhanced polymeric entanglements, thereby increasing the resistance to deformation. Increasing probe speed increased formulation hardness in a linear manner, because of the effects of probe speed on probe displacement and surface area. The relationship between formulation hardness and probe displacement was linear and was dependent on probe speed. Furthermore, the proportionality constant (gel strength) increased as a function of PVP concentration. The relationship between formulation hardness and diameter ratio was biphasic and was statistically defined by two linear relationships relating to diameter ratios from 0 to 0.4 and from 0.4 to 0.563. The dramatically increased hardness, associated with diameter ratios in excess of 0.4, was accredited to boundary effects, that is, the effect of the container wall on product flow. Using dimensional analysis, the hardness and probe displacement in TPA were mathematically transformed into corresponding rheological parameters, namely shearing stress and rate of shear, thereby allowing the application of the power law (??=?k?n) to textural data. Importantly, the consistencies (k) of the formulations, calculated using transformed textural data, were statistically similar to those obtained using flow rheometry. In conclusion, this study has, firstly, characterized the relationships between textural data and two key instrumental parameters in TPA and, secondly, described a method by which rheological information may be derived using this technique. This will enable a greater application of TPA for the rheological characterization of pharmaceutical gels and, in addition, will enable efficient interpretation of textural data under different experimental parameters.

Fundamentals of steady state, non-Newtonian rimming flow

Rimming flow on the inner surface of a horizontal rotating cylinder is investigated. Using a scale analysis, a theoretical description is obtained for steady-state non-Newtonian flow. Simple lubrication theory is applied since the Reynolds number is small and the liquid film is thin. Since the Deborah number is very small the flow is viscometric. The shear-thinning number, which characterizes the shear-thinning effect, may be small or large. A general constitutive law for this kind of flow requires only a single function relating shear stress and shear rate that corresponds to a generalized Newtonian liquid. For this case the run-off condition for rimming flow is derived. Provided the run-off condition is satisfied, the existence of a continuous steady-state solution is proved. The rheological models, which show Newtonian behavior at low shear rates with transition to power-law shear thinning at moderate shear rates, are considered. Numerical results are carried out for the Carreau and Ellis models, which exhibit Newtonian behavior near the free surface and power-law behavior near the wall of the rotating cylinder.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

An Experimental and Predictive Evaluation of Unsteady Gas Flow through Automotive Catalyst Elements

The incorporation of one-dimensional simulation codes within engine modelling applications has proved to be a useful tool in evaluating unsteady gas flow through elements in the exhaust system. This paper reports on an experimental and theoretical investigation into the behaviour of unsteady gas flow through catalyst substrate elements. A one-dimensional (1-D) catalyst model has been incorporated into a 1-D simulation code to predict this behaviour.

Experimental data was acquired using a ‘single pulse’ test rig. Substrate samples were tested under ambient conditions in order to investigate a range of regimes experienced by the catalyst during operation. This allowed reflection and transmission characteristics to be quantified in relation to both geometric and physical properties of substrate elements. Correlation between measured and predicted results is demonstrably good and the model provides an effective analysis tool for evaluating unsteady gas flow through different catalytic converter designs.