Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Modelling of species in hood fires by conditional moment closure

Carbon monoxide, the chief killer in fires, and other species are modelled for a series of enclosure fires. The conditions emulate building fires where CO is formed in the rich, turbulent, nonpremixed flame and is transported frozen to lean mixtures by the ceiling jet which is cooled by radiation and dilution. Conditional moment closure modelling is used and computational domain minimisation criteria are developed which reduce the computational cost of this method. The predictions give good agreement for CO and other species in the lean, quenched-gas stream, holding promise that this method may provide a practical means of modelling real, three-dimensional fire situations. (c) 2005 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Numerical simulation of dendrite growth during solidification

A comprehensive probabilistic model for simulating dendrite morphology and investigating dendritic growth kinetics during solidification has been developed, based on a modified Cellular Automaton (mCA) for microscopic modeling of nucleation, growth of crystals and solute diffusion. The mCA model numerically calculated solute redistribution both in the solid and liquid phases, the curvature of dendrite tips and the growth anisotropy. This modeling takes account of thermal, curvature and solute diffusion effects. Therefore, it can simulate microstructure formation both on the scale of the dendrite tip length. This model was then applied for simulating dendritic solidification of an Al-7%Si alloy. Both directional and equiaxed dendritic growth has been performed to investigate the growth anisotropy and cooling rate on dendrite morphology. Furthermore, the competitive growth and selection of dendritic crystals have also investigated.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

A Dynamical Modeling to Study the Adaptive Immune System and the Influence of Antibodies in the Immune Memory

Immunological systems have been an abundant inspiration to contemporary computer scientists. Problem solving strategies, stemming from known immune system phenomena, have been successfully applied to chall enging problems of modem computing. Simulation systems and mathematical modeling are also beginning use to answer more complex immunological questions as immune memory process and duration of vaccines, where the regulation mechanisms are not still known sufficiently (Lundegaard, Lund, Kesmir, Brunak, Nielsen, 2007). In this article we studied in machina a approach to simulate the process of antigenic mutation and its implications for the process of memory. Our results have suggested that the durability of the immune memory is affected by the process of antigenic mutation.and by populations of soluble antibodies in the blood. The results also strongly suggest that the decrease of the production of antibodies favors the global maintenance of immune memory.

Analysis of flexural waves excited by rectangular contact type transducers

In this paper, an attempt was made to investigate a fundamental problem related to the flexural waves excited by rectangular transducers. Due to the disadvantages of the Green's function approach for solving this problem, a direct and effective method is proposed using a multiple integral transform method and contour integration technique. The explicit frequency domain solutions obtained from this newly developed method are convenient for understanding transducer behavior and theoretical optimization and experimental calibration of rectangular transducers. The time domain solutions can then be easily obtained by using the fast Fourier transform technique. (C) 2001 Elsevier Science B.V. All rights reserved.

The measurement of Pitot pressure in high enthalpy expansion tubes

Expansion tubes operating at total flow enthalpies of 100 MJ kg(-1) or more have characteristical test times of 30-50 mus. Under these conditions, the response time of the Pitot pressure measuring device is critical when performing flow calibration studies. The conventional technique of using a commercial pressure transducer protected by shielding has not always proven to be effective, due to the relatively large (and variable) response time caused by the shielding. A device called the stress wave bar gauge has been designed and calibrated and shown to be an effective way to measure the Pitot pressure with a response time of only 2-3 mus.

A comparison of thermal condition between pilot- and full-scale furnaces for studying slagging and fouling propensity in PF boilers

The Australian Coal Industry Research Laboratory (ACIRL) furnace is scaled to simulate slagging and fouling in operating boilers. This requires that the gas and target temperatures, the heat flux, and the flow pattern be the same as those in real boilers. The gas and target temperatures are maintained by insulating the wall and cooling the target respectively. The flow pattern of a small burner cannot be the same as a large furnace. However, this flow pattern is partially compensated for by placing the slagging panels in three vertical locations. The paper develops the models of radiant heat transfer from the flame to the deposits both in pilot-scale and full-scale furnaces. They are used to compare the effective radiant heat transfer of the pilot- and full-scale furnaces. The experimental data both from the pilot- and full-scale furnaces are used to verify the incident heat flux and temperature profiles in the pilot- and full-scale furnaces. The results showed that the thermal condition in the pilot-scale furnace meets the requirements for studying the slagging regarding the gas temperature and the incident heat flux, particularly for the panel #1. The gas temperature in the convective section also meets the requirement for studying the fouling.

A methodology for discipline-specific curriculum development

A methodology and framework for discipline-specific curriculum development in a local context is described. These activities, as part of the Thailand-Australia Science and Engineering Assistance Project, were in response to a needs analysis for curriculum assistance to a number of publicly-funded Thai universities in the engineering priority area of Materials Processing and Manufacturing. The paper outlines a strategy for the delivery of a centralised curriculum development workshop for academic staff follow-up visits and local curriculum activities with participating universities, and the presentation of technical short courses as guidance for such activity in other settings and/or discipline areas. This paper is part of a process of documentation so that others can apply the developed methodology and framework for curriculum development. While the paper is a report on curriculum activities in a particular setting, it is written in a manner that allows application of the methodology to other settings. The reader is advised that each curriculum activity needs to adopt a methodology and strategy to fit the particular circumstances being considered To assist in applying this approach elsewhere, a description of the various steps in the curriculum process, and typical responses to some of the more global issues, have been presented. Full details are available in the various TASEAP reports prepared by the authors. Specific detail has been omitted where this detail does not provide any information for generalized consumption.

A flexible method for rapid force computation in elliptical magnets

The design of open-access elliptical cross-section magnet systems has recently come under consideration. Obtaining values for the forces generated within these unusual magnets is important to progress the designs towards feasible instruments. This paper presents a novel and flexible method for the rapid computation of forces within elliptical magnets. The method is demonstrated by the analysis of a clinical magnetic resonance imaging magnet of elliptical cross-section and open design. The analysis reveals the non-symmetric nature of the generated Maxwell forces, which are an important consideration, particularly in the design of superconducting systems.

An inverse design method for elliptical MRI magnets

An inverse, current density mapping (CDM) method has been developed for the design of elliptical cross-section MRI magnets. The method provides a rapid prototyping system for unusual magnet designs, as it generates a 3D current density in response to a set of target field and geometric constraints. The emphasis of this work is on the investigation of new elliptical coil structures for clinical MRI magnets. The effect of the elliptical aspect ratio on magnet performance is investigated. Viable designs are generated for symmetric, asymmetric and open architecture elliptical magnets using the new method. Clinically relevant attributes such as reduced stray field and large homogeneous regions relative to total magnet length are included in the design process and investigated in detail. The preliminary magnet designs have several novel features.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Synthesis of polymer-montmorillonite nanocomposites by in situ intercalative polymerization

Two polymer-montmorillonite (MMT) nanocomposites have been synthesized by in situ intercalative polymerization. The styrene monomer is intercalated into the interlayer space of organically modified MMT, a layered clay mineral. Upon the intercalation, the complex is subsequently polymerized in the confinement environment of the interlayer space with a free radical initiator, 2,2-azobis isobutyronitrile. The aniline monomer is also intercalated and then polymerized within the interlayer space of sodium- and copper-MMT initiated by ammonium peroxodisulphate and interlayer copper cations respectively. X-ray diffraction indicates that the MMT layers are completely dispersed in the polystyrene matrix and an exfoliated structure has been obtained. The resulting polyaniline-MMT nanocomposites show a highly ordered structure of a single polyaniline layer stacked with the MMT layers. Fourier transform infrared spectra further confirm the intercalation and formation of both polymer-MMT nanocomposites.