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.

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.

Pseudo-spectral solutions for fluid flow and heat transfer in electro-metallurgical applications

The pseudo-spectral solution method offers a flexible and fast alternative to the more usual finite element and volume methods, particularly when the long-time transient behaviour of a system is of interest. The exact solution is obtained at grid collocation points leading to superior accuracy on modest grids. Furthermore, the grid can be freely adapted in time and space to particular flow conditions or geometric variations, especially useful where strongly coupled, time-dependent, multi-physics solutions are investigated. Examples include metallurgical applications involving the interaction of electromagnetic fields and conducting liquids with a free surface. The electromagnetic field determines the instantaneous liquid volume shape, which then affects the electromagnetic field. A general methodology of the pseudo-spectral approach is presented, with several instructive example applications: the aluminium electrolysis MHD problem, induction melting in a cold crucible and the dynamics of AC/DC magnetically levitated droplets. Finally, comparisons with available analytical solutions and to experimental measurements are discussed.

A computational study of the characteristics of aircraft post-crash fires

Full-scale furnished cabin fires have been studied experimentally for the purpose of characterising the post-crash cabin fire environment by the US Federal Aviation Administration for many years. In this paper the Computational Fluid Dynamics fire field model SMARTFIRE is used to simulate one of these fires conducted in the C-133 test facility in order to provide further validation of the computational approach and the SMARTFIRE software. The experiment involves exposing the interior cabin materials to an external fuel fire, opening only one exit at the far end of the cabin (the same side as the rupture) for ventilation, and noting the subsequent spread of the external fire to the cabin interior and the onset of flashover at approximately 210 seconds. Through this analysis, the software is shown to be in good agreement with the experimental data, producing reasonable agreement with the fire dynamics prior to flashover and producing a reasonable prediction of the flashover time i.e. 225 seconds. The paper then proceeds to utilize the model to examine the impact on flashover time of the extent of cabin furnishings and cabin ventilation provided by available exits

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.

Run off Conditions on Rimming Power Law Fluids

The rimming ?ow of a power-law ?uid in the inner surface of a horizontal rotating cylinder is investigated. Exploiting the fact that the liquid layer is thin, the simplest lubrication theory is applied. The generalized run-off condition for the steady-state ?ow of the power-law liquid is derived. In the bounds implied by this condition, ?lm thickness admits a continuous solution. In the supercritical case when the mass of non-Newtonian liquid exceeds a certain value or the speed of rotation is less than an indicated limit, a discontinuous solution is possible and a hydraulic jump may occur in the steady-state regime. The location and height of the hydraulic jump for the power-law liquid is determined.

Design of a thick-walled screen for flow equalization in microstructured reactors

A systematic computational fluid dynamics (CFD) approach has been applied to design the geometry of the channels of a three-dimensional (thick-walled) screen comprising upstream and downstream sets of elongated channels positioned at an angle of 90 degrees with respect to each other. Such a geometry of the thick-wall screen can effectively drop the ratio of the maximum flow velocity to mean flow velocity below 1.005 in a downstream microstructured reactor at low Reynolds numbers. In this approach the problem of flow equalization reduces to that of flow equalization in the first and second downstream channels of the thick-walled screen. In turn, this requires flow equalization in the corresponding cross-sections of the upstream channels. The validity of the proposed design method was assessed through a case study. The effect of different design parameters on the flow non-uniformity in the downstream channels has been established. The design equation is proposed to calculate the optimum values of the screen parameters. The CFD results on flow distribution were experimentally validated by Laser Doppler Anemometry measurements in the range of Reynolds numbers from 6 to 113. The measured flow non-uniformity in the separate reactor channels was below 2%.

The Characteristics of Wave Impacts on an Oscillating Wave Surge Converter

Wave impacts on an oscillating wave surge converter are examined using experimental and numerical methods. The mechanics of the impact event are identified experimentally with the use of images recorded with a high-speed camera. It is shown that it is the device that impacts the wave rather than a breaking wave impacting the device. Numerical simulations using two different approaches are used to further understand the issue. Good agreement is shown between numerical simulations and experimental measurements at 25th scale.

Development and validation of Fluid Structure Interaction Methods for an Oscillating Wave Surge Converter

The term fatigue loads on the Oyster Oscillating Wave Surge Converter (OWSC) is used to describe hydrostatic loads due to water surface elevation with quasi-static changes of state. Therefore a procedure to implement hydrostatic pressure distributions into finite element analysis of the structure is desired. Currently available experimental methods enable one to measure time variant water surface elevation at discrete locations either on or around the body of the scale model during tank tests. This paper discusses the development of a finite element analysis procedure to implement time variant, spatially distributed hydrostatic pressure derived from discretely measured water surface elevation. The developed method can process differently resolved (temporal and spatial) input data and approximate the elevation over the flap faces with user defined properties. The structural loads, namely the forces and moments on the body can then be investigated by post processing the numerical results. This method offers the possibility to process surface elevation or hydrostatic pressure data from computational fluid dynamics simulations and can thus be seen as a first step to a fluid-structure interaction model.

The Characteristics of Wave Impacts on an Oscillating Wave Surge Converter

Wave impacts on an Oscillating Wave Surge Converter are examined using experimental and numerical methods. The mechanics of the impact event are identified experimentally with the use of images recorded with a high speed camera. It is shown that it is the device which impacts the wave rather than a breaking wave impacting the device. Numerical simulations using two different approaches are used to further understand the issue. Good agreement is shown between numerical simulations and experimental measurements at 25th scale.

Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

Coletores de distribuição e sua otimização hidráulica, aplicados a uma instalação de produção de água arrefecida em sistemas de AVAC ( Fórum Picoas)

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.