Classification and analysis of integrating frameworks in multiscale modelling

Chemical engineers are turning to multiscale modelling to extend traditional modelling approaches into new application areas and to achieve higher levels of detail and accuracy. There is, however, little advice available on the best strategy to use in constructing a multiscale model. This paper presents a starting point for the systematic analysis of multiscale models by defining several integrating frameworks for linking models at different scales. It briefly explores how the nature of the information flow between the models at the different scales is influenced by the choice of framework, and presents some restrictions on model-framework compatibility. The concepts are illustrated with reference to the modelling of a catalytic packed bed reactor. (C) 2004 Elsevier Ltd. All rights reserved.

Evaluation of the BEASY program using linear and piecewise linear approaches for the boundary conditions

The boundary element method (BEM) was used to study galvanic corrosion using linear and logarithmic boundary conditions. The linear boundary condition was implemented by using the linear approach and the piecewise linear approach. The logarithmic boundary condition was implemented by the piecewise linear approach. The calculated potential and current density distribution were compared with the prior analytical results. For the linear boundary condition, the BEASY program using the linear approach and the piecewise linear approach gave accurate predictions of the potential and the galvanic current density distributions for varied electrolyte conditions, various film thicknesses, various electrolyte conductivities and various area ratio of anode/cathode. The 50-point piecewise linear method could be used with both linear and logarithmic polarization curves.

A diffusion and T2 relaxation MRI study of the ovine lumbar intervertebral disc under compression in vitro

The ovine lumbar intervertebral disc is a useful model for the human lumbar disc. We present preliminary estimates of diffusion coefficients and T-2 relaxation times in a pilot MRI study of the ovine lumbar intervertebral disc during uniaxial compression in vitro, and identify factors that hamper the ability to accurately monitor the temporal evolution of the effective diffusion tensor at high spatial resolution.

Microstructure of MmM(5)/Mg multi-layer hydrogen storage films prepared by magnetron sputtering

Multi-layer hydrogen storage thin films with Mg and MmNi(3.5)(CoAlMn)(1.5) (here Mm denotes La-rich mischmetal) as alternative layers were prepared by direct current magnetron sputtering. Transmission electron microscopy investigation shows that the microstructure of the MmNi(3.5)(CoAlMn)(1.5) and Mg layers are significantly different although their deposition conditions are the same. The MmNi(3.5)(CoAlMn)(1.5) layer is composed of two regions: one is an amorphous region approximately 4 nm thick at the bottom of the layer and the other is a nanocrystalline region on top of the amorphous region. The Mg layer is also composed of two regions: one is a randomly orientated nanocrystalline region 50 nm thick at the bottom of the layer and the other is a columnar crystallite region on top of the nanocrystalline region. These Mg columnar crystallites have their [001] directions parallel to the growth direction and the average lateral size of these columnar crystallites is about 100 nm. A growth mechanism of the multi-layer thin films is discussed based on the experiment results. Wiley-Liss, Inc.

Superplasticity and superplastic forming ability of a Zr-Ti-Ni-Cu-Be bulk metallic glass in the supercooled liquid recgion

The superplastic deformation behavior and superplastic forming ability of the Zr41.25Ti13.75Ni10Cu12.5Be22.5 (at.%) bulk metallic glass (BMG) in the supercooled liquid region were investigated. The isothermal tensile results indicate (hat the BMG exhibits a Newtonian behavior at low strain rates but a non-Newtonian behavior at hiqh-strain rates in the initial deformation stage. The maximum elongation reaches as high as 1624% at 656 K. and nanocrystallization was found to occur during the deformation process. Based cm the analysis on tensile deformation. a gear-like micropart is successfully die-forged via a superplastic forgings process. demonstrating that the BMG has excellent workability in the supercooled liquid region. (C) 2004 Elsevier B.V. All rights reserved.

Effects of hydrogen-air non-equilibrium chemistry on the performance of a model scramjet thrust nozzle

Two aspects of hydrogen-air non-equilibrium chemistry related to scramjets are nozzle freezing and a process called 'kinetic afterburning' which involves continuation of combustion after expansion in the nozzle. These effects were investigated numerically and experimentally with a model scramjet combustion chamber and thrust nozzle combination. The overall model length was 0.5m, while precombustion Mach numbers of 3.1 +/- 0.3 and precombustion temperatures ranging from 740K to 1,400K were involved. Nozzle freezing was investigated at precombustion pressures of 190kPa and higher, and it was found that the nozzle thrusts were within 6% of values obtained from finite rate numerical calculations, which were within 7% of equilibrium calculations. When precombustion pressures of 70kPa or less were used, kinetic afterburning was found to be partly responsible for thrust production, in both the numerical calculations and the experiments. Kinetic afterburning offers a means of extending the operating Mach number range of a fixed geometry scramjet.

Numerical modelling of radiating superorbital flows

Experimental aerodynamic studies of the flows around new aerocapture spacecraft configurations are presently being done in the superorbital expansion tubes at The University of Queensland. Short duration flows at speeds of 10--13 km/s are produced in the expansion tube facility and are then applied to the model spacecraft. Although high-temperature effects, such as molecular dissociation, have long been a part of the computational modelling of the expansion tube flows for speeds below 10 km/s, radiation may now be a significant mechanism of energy transfer within the shock layer on the model. This paper will study the coupling of radiation energy transport for an optically thin gas to the flow dynamics in order to obtain accurate predictions of thermal loads on the spacecraft. The results show that the effect of radiation on the flowfields of subscale models for expansion tube experiments can be assessed by measurements of total heat transfer and radiative heat transfer.

Fuel cells, hydrogen and energy supply in Australia

Australia is unique in terms of its geography, population distribution, and energy sources. It has an abundance of fossil fuel in the form of coal, natural gas, coal seam methane (CSM), oil, and a variety renewable energy sources that are under development. Unfortunately, most of the natural gas is located so far away from the main centres of population that it is more economic to ship the energy as LNG to neighboring countries. Electricity generation is the largest consumer of energy in Australia and accounts for around 50% of greenhouse gas emissions as 84% of electricity is produced from coal. Unless these emissions are curbed, there is a risk of increasing temperatures throughout the country and associated climatic instability. To address this, research is underway to develop coal gasification and processes for the capture and sequestration Of CO2. Alternative transport fuels such as biodiesel are being introduced to help reduce emissions from vehicles. The future role of hydrogen is being addressed in a national study commissioned this year by the federal government. Work at the University of Queensland is also addressing full-cycle analysis of hydrogen production, transport, storage, and utilization for both stationary and transport applications. There is a modest but growing amount of university research in fuel cells in Australia, and an increasing interest from industry. Ceramic Fuel Cells Ltd. (CFCL) has a leading position in planar solid oxide fuel cells (SOFCs) technology, which is being developed for a variety of applications, and next year Perth in Western Australia is hosting a trial of buses powered by proton-exchange fuel cells. (C) 2004 Elsevier B.V. All rights reserved.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Numerical evaluation of three dimensional effects in planar flow birefringence

The influence of three dimensional effects on isochromatic birefringence is evaluated for planar flows by means of numerical simulation. Two fluid models are investigated in channel and abrupt contraction geometries. In practice, the flows are confined by viewing windows, which alter the stresses along the optical path. The observed optical properties differ therefore from their counterpart in an ideal two-dimensional flow. To investigate the influence of these effects, the stress optical rule and the differential propagation Mueller matrix are used. The material parameters are selected so that a retardation of multiple orders is achieved, as is typical for highly birefringent melts. Errors due to three dimensional effects are mainly found on the symmetry plane, and increase significantly with the flow rate. Increasing the geometric aspect ratio improve the accuracy provided that the error on the retardation is less than one order. (C) 2004 Elsevier B.V. All rights reserved.

Numerical simulation of abrupt contraction flows using the Double Convected Pom-Pom model

The Double Convected Pom-Pom model was recently introduced to circumvent some numerical and theological defects found in other formulations of the Pom-Pom concept. It is used here for the simulation of a benchmark problem: the flow in an abrupt planar contraction. The predictions are compared with birefringence measurements and show reasonable quantitative agreement with experimental data. A parametric study is also carried out with the aim of analysing the effect of the branching parameter on vortex dynamics and extrudate swell. The results show that the Double Convected Pom-Pom model (DCPP) model is able to discriminate between branched and linear macromolecular structures in accordance with experimental observations. In that respect, the role of the extensional properties in determining complex flow behaviour is stressed. Also, the ratio of the first normal stress difference to the shear stress appears to play a major role in die swell observation. For the time being, the role of the second normal stress difference appears to be less obvious to evaluate in this complex flow. (C) 2004 Elsevier B.V. All rights reserved.

Permeation of binary gas mixtures in ultramicroporous membranes

High-quality nanometer thick ultramicroporous membranes were prepared from silica sol-gel processes and tested for the permeation of binary gas mixtures of He, H-2, CO2, and CH4 across different temperature and partial pressure regimens. Pore size distribution by molecular probing showed that the majority of pore sizes had dimensions below 2.9 Angstrom. In 50:50 binary mixtures, the fluxes of gases increased as a function of temperature, indicating an activated transport mechanism. The ultramicroporous membranes showed high selectivities at 150 degreesC for He/CO2 (30), He/CH4 (93), H-2/CO2 (10), and H-2/CH4 (9) with lower selectivities for CO2/CH4 (5). High activation energies (E-a) were observed for the permeance of 50:50 binary mixtures containing He and H-2 of 22.1-27.5 and 17.6-23.1 kJ.mol(-1), respectively. The E-a for the permeance of the total mixture approached the E-a for the permeance of the molecule with the smaller kinetic diameter (He or H-2).