The numerical simulation of enclosure fires using a Cfd fire field model coupled with a pyrolysis based solid fuel combustion submodel—a first approximation

In this paper, we present some early work concerned with the development of a simple solid fuel combustion model incorporated within a Computational Fluid Dynamics (CFD) framework. The model is intended for use in engineering applications of fire field modeling and represents an extension of this technique to situations involving the combustion of solid cellulosic fuels. A simple solid fuel combustion model consisting of a thermal pyrolysis model, a six flux radiation model and an eddy-dissipation model for gaseous combustion have been developed and implemented within the CFD code CFDS-FLOW3D. The model is briefly described and demonstrated through two applications involving fire spread in a compartment with a plywood lined ceiling. The two scenarios considered involve a fire in an open and closed compartment. The model is shown to be able to qualitatively predict behaviors similar to "flashover"—in the case of the open room—and "backdraft"— in the case of the initially closed room.

Simulation of transverse and longitudinal magnetic ripple structures induced by surface anisotropy

Micromagnetic ripple structures on the surfaces of thick specimens of ultra-soft magnetic material having strong surface anisotropy Ks favouring out-of-surface magnetization have been calculated. These ripples have wavelengths of the order of 0.1 μm and extend to a depth ∼ √A/Ms, where A is the exchange constant and Ms is the saturation magnetization. The wave-vectors of the ripple structures are either transverse or parallel to the bulk magnetization. Both structures have lower energy than the one-dimensional structure discussed by O'Handley and Woods, and they exhibit stronger normal magnetization. The transverse structure requires a surface anisotropy Ks ≥ 0.80K0, where is that required for the one-dimensional structure. The threshold for longitudinal ripples is 0.84K0. It is suggested that the transverse structure probably constitutes the ground state. The magnitudes of Ks and A should be obtainable from measurements of the ripple wavelength and amplitude, and Ms.

Coupled 3-D finite difference time domain and finite volume methods for solving microwave heating in porous media

Computational results for the microwave heating of a porous material are presented in this paper. Combined finite difference time domain and finite volume methods were used to solve equations that describe the electromagnetic field and heat and mass transfer in porous media. The coupling between the two schemes is through a change in dielectric properties which were assumed to be dependent both on temperature and moisture content. The model was able to reflect the evolution of temperature and moisture fields as the moisture in the porous medium evaporates. Moisture movement results from internal pressure gradients produced by the internal heating and phase change.

Modelling the fluid dynamics and coupled phenomena in arc weld pools

A 3D model of melt pool created by a moving arc type heat sources has been developed. The model solves the equations of turbulent fluid flow, heat transfer and electromagnetic field to demonstrate the flow behaviour phase-change in the pool. The coupled effects of buoyancy, capillary (Marangoni) and electromagnetic (Lorentz) forces are included within an unstructured finite volume mesh environment. The movement of the welding arc along the workpiece is accomplished via a moving co-ordinator system. Additionally a method enabling movement of the weld pool surface by fluid convection is presented whereby the mesh in the liquid region is allowed to move through a free surface. The surface grid lines move to restore equilibrium at the end of each computational time step and interior grid points then adjust following the solution of a Laplace equation.

A generic time and space accurate numerical approach to closely coupled fluid-structure interaction problems

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, flow in elastic pipes and blood vessels and extrusion of metals through dies. However a comprehensive computational model of these multi-physics phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply even to the extent in metal forming, for example, that the deformation of the die is totally ignored. More recently, strategies for solving the full coupling between the fluid and soild mechanics behaviour have developed. Conventionally, the computational modelling of fluid structure interaction is problematical since computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. In the past the concurrent, but rather disparate, development paths for the finite element and finite volume methods have resulted in numerical software tools for CFD and CSM that are different in almost every respect. Hence, progress is frustrated in modelling the emerging multi-physics problem of fluid structure interaction in a consistent manner. Unless the fluid-structure coupling is either one way, very weak or both, transferring and filtering data from one mesh and solution procedure to another may lead to significant problems in computational convergence. Using a novel three phase technique the full interaction between the fluid and the dynamic structural response are represented. The procedure is demonstrated on some challenging applications in complex three dimensional geometries involving aircraft flutter, metal forming and blood flow in arteries.

A computational model of coupled heat and moisture transfer with phase change in granular sugar during varying environmental conditions

As part of a comprehensive effort to predict the development of caking in granular materials, a mathematical model is introduced to model simultaneous heat and moisture transfer with phase change in porous media when undergoing temperature oscillations/cycling. The resulting model partial differential equations were solved using finite-volume procedures in the context of the PHYSICA framework and then applied to the analysis of sugar in storage. The influence of temperature on absorption/desorption and diffusion coefficients is coupled into the transport equations. The temperature profile, the depth of penetration of the temperature oscillation into the bulk solid, and the solids moisture content distribution were first calculated, and these proved to be in good agreement with experimental data. Then, the influence of temperature oscillation on absolute humidity, moisture concentration, and moisture migration for different parameters and boundary conditions was examined. As expected, the results show that moisture near boundary regions responds faster than farther away from them with surface temperature changes. The moisture absorption and desorption in materials occurs mainly near boundary regions (where interactions with the environment are more pronounced). Small amounts of solids moisture content, driven by both temperature and vapour concentration gradients, migrate between boundary and center with oscillating temperature.

Use of a coupled mechanical-acoustic computational model to identify failure mechanisms in paper production

In this paper, a couple mechanical-acoustic system of equations is solved to determine the relationship between emitted sound and damage mechanisims in paper under controlled stress conditions. The simple classical expression describing the frequency of a plucked string to its material properties is used to generate a numberical representation of the microscopic structue of the paper, and the resulting numerical model is then used to simulate the vibration of a range of simple fibre structures when undergoing two distinct types of damange mechanisms: (a)fibre/fibre bond failure, (b) fibre failure. The numercial results are analysed to determine whether there is any detectable systematic difference between the resulting acoustic emissions of the two damage processes. Fourier techniques are then used to compare th computeed results against experimental measurements. Distinct frequency components identifying each type of damage are shown to exist, and in this respect theory and experiments show good correspondece. Hence, it is shown, that althrough the mathematical model represents a grossly-simplified view of the complex structure of the paper, it nevertheless provides a good understanding of the underlying micro-mechanisms characterising its proeperties as a stress-resisting structure. Use of the model and acoompanying software will enable operators to identify approaching failure conditions in the continuous production of paper from emitted sound signals and take preventative action.

Predicting HCl concentrations in fire enclosures using an HCl decay model coupled to a CFD-based fire field model

The amount of atmospheric hydrogen chloride (HCl) within fire enclosures produced from the combustion of chloride-based materials tends to decay as the fire effluent is transported through the enclosure due to mixing with fresh air and absorption by solids. This paper describes an HCl decay model, typically used in zone models, which has been modified and applied to a computational fluid dynamics (CFD)-based fire field model. While the modified model still makes use of some empirical formulations to represent the deposition mechanisms, these have been reduced from the original three to two through the use of the CFD framework. Furthermore, the effect of HCl flow to the wall surfaces on the time to reach equilibrium between HCl in the boundary layer and on wall surfaces is addressed by the modified model. Simulation results using the modified HCl decay model are compared with data from three experiments. The model is found to be able to reproduce the experimental trends and the predicted HCl levels are in good agreement with measured values

Coupled FDTD-FVM simulation of microwave heating of polymer materials for micro-systems packaging applications

Heating in an idealised polymer load in a novel open-ended variable frequency microwave oven is numerically simulated using a couple solver approach. The frequency-agile microwave oven bonding system (FAMOBS)is developed to meet rapid polymer curing requirements in microelectronics and optoelectronics manufacturing. The heating of and idealised polymer load has been investigated through numerical modelling. Assessment of the system comprises of simulation of electromagnetic fields and of temperature distribution within the load. Initial simulation results are presented and contrasted with experimental analysis of field distribution

Numerical simulation of microwave thawing using a coupled solver approach

Thawing of a frozen food product in a domestic microwave oven is numerically simulated using a coupled solver approach. The approach consists of a dedicated electromagnetic FDTD solver and a closely coupled UFVM multi-physics package. Two overlapping numerical meshes are defined; the food material and container were meshed for heat transfer and phase change solution, whilst the microwave oven cavity and waveguide were meshed for the microwave irradiation. The two solution domains were linked using a cross-mapping routine. This approach allowed the rotation of the food load to be captured. Power densities obtained on the structured FDTD mesh were interpolated onto the UFVM mesh for each timestep/turntable position. The UFVM solver utilised the power density data to advance the temperature and phase distribution solution. The temperature-dependant dielectric and thermo-physical properties of the food load were updated prior to revising the electromagnetic solution. Changes in thermal/electric properties associated with the phase transition were fully accounted for as well as heat losses from product to cavity. Two scenarios were investigated: a centric and eccentric placement on the turntable. Developing temperature fields predicted by the numerical solution are validated against experimentally obtained data. Presented results indicate the feasibility of fully coupled simulations of the microwave heating of a frozen product. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Multiphysics simulation of microwave processing using a multidomain coupled solver approach

Microwave processing of materials is numerically simulated using a coupled solver approach. Microwave heating is a complex coupled process due to the variation in dielectric properties during heating. The effects of heating an object in a electromagnetic field directly influence the manner in which it interacts with the field. Simplifying assumptions and empirical solutions do not capture the fundamental physics involved and, in general, do not provide usefully accurate solutions in a number of practical problems. In order to capture the underlying processes involved in microwave heating, the problem must be looked at in a holistic manner rather than a number of discrete processes. This contribution outlines a coupled-solver multiphysics analysis approach to the solution of practical microwave heating problems.

Role of bonding time and temperature on the physical properties of coupled anisotropic conductive–nonconductive adhesive film for flip chip on glass technology

Using thermosetting epoxy based conductive adhesive films for the flip chip interconnect possess a great deal of attractions to the electronics manufacturing industries due to the ever increasing demands for miniaturized electronic products. Adhesive manufacturers have taken many attempts over the last decade to produce a number of types of adhesives and the coupled anisotropic conductive-nonconductive adhesive film is one of them. The successful formation of the flip chip interconnection using this particular type of adhesive depends on, among factors, how the physical properties of the adhesive changes during the bonding process. Experimental measurements of the temperature in the adhesive have revealed that the temperature becomes very close to the required maximum bonding temperature within the first 1s of the bonding time. The higher the bonding temperature the faster the ramp up of temperature is. A dynamic mechanical analysis (DMA) has been carried out to investigate the nature of the changes of the physical properties of the coupled anisotropic conductive-nonconductive adhesive film for a range of bonding parameters. Adhesive samples that are pre-cured at 170, 190 and 210°C for 3, 5 and 10s have been analyzed using a DMA instrument. The results have revealed that the glass transition temperature of this type of adhesive increases with the increase in the bonding time for the bonding temperatures that have been used in this work. For the curing time of 3 and 5s, the maximum glass transition temperature increases with the increase in the bonding temperature, but for the curing time of 10s the maximum glass transition temperature has been observed in the sample which is cured at 190°C. Based on these results it has been concluded that the optimal bonding temperature and time for this kind of adhesive are 190°C and 10s, respectively.

Coupled fire/evacuation analysis of the Station Nightclub fire

In this paper, coupled fire and evacuation simulation tools are used to simulate the Station Nightclub fire. This study differs from the analysis conducted by NIST in three key areas; (1)an enhanced flame spread model and (2)a toxicity generation model are used, (3)the evacuation is coupled to the fire simulation. Predicted early burning locations in the full-scale fire simulation are in line with photographic evidence and the predicted onset of flashover is similar to that produced by NIST. However, it is suggested that both predictions of the flashover time are approximately 15 sec earlier than actually occurred. Three evacuation scenarios are then considered, two of which are coupled with the fire simulation. The coupled fire and evacuation simulation suggests that 180 fatalities result from a building population of 460. With a 15 sec delay in the fire timeline, the evacuation simulation produces 84 fatalities which are in good agreement with actual number of fatalities. An important observation resulting from this work is that traditional fire engineering ASET/RSET calculations which do not couple the fire and evacuation simulations have the potential to be considerably over optimistic in terms of the level of safety achieved by building designs.

Elemental analysis of soil samples for forensic purposes by inductively coupled plasma spectrometry - precision considerations

Major and trace elemental composition provides a powerful basis for forensic comparison of soils, sediments and rocks. However, it is important that the potential 'errors' associated with the procedures are fully understood and quantified, and that standard protocols are applied for sample preparation and analysis. This paper describes such a standard procedure and reports results both for instrumental measurement precision (repeatability) and overall 'method' precision (reproducibility). Results obtained both for certified reference materials and example soils show that the instrumental measurement precision (defined by the coefficient of variation, CV) for most elements is better than 2-3%. When different solutions were prepared from the same sample powder, and from different sub-sample powders prepared from the same parent sample, the CV increased to c. 5-6% for many elements. The largest variation was found in results for certified reference materials generated from 23 instrument runs over an 18 month period (mean CV=c. 11%). Some elements were more variable than others. W was found to be the most variable and the elements V, Cr, Co, Cu, Ni and Pb also showed higher than average variability. SiO2, CaO, Al2O3 and Fe2O3, Rb, Sr, La, Ce, Nd and Sm generally showed lower than average variability, and therefore provided the most reliable basis for inter-sample comparison. It is recommended that, whenever possible, samples relating to the same investigation should be analysed in the same sample run, or at least sequential runs.