Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Simulation and optimization of heated, inviscid flows in scramjet ducts

A space-marching code for the simulation and optimization of inviscid supersonic flow in three dimensions is described. The now in a scramjet module with a relatively complex three-dimensional geometry is examined and wall-pressure estimates are compared with experimental data. Given that viscous effects are not presently included, the comparison is reasonable. The thermodynamic compromise of adding heat in a diverging combustor is also examined. The code is then used to optimize the shape of a thrust surface for a simpler (box-section) scramjet module in the presence of uniform and nonuniform heat distributions. The optimum two-dimensional profiles for the thrust surface are obtained via a perturbation procedure that requires about 30-50 now solutions. It is found that the final shapes are fairly insensitive to the details of the heat distribution.

Surface oxide and the role of magnesium during the sintering of aluminum

The effect of trace additions of magnesium on the sintering of aluminum and its alloys is examined. Magnesium, especially at low concentrations, has a disproportionate effect on sintering because it disrupts the passivating Al2O3 layer through the formation of a spinel phase. Magnesium penetrates the sintering compact by solid-state diffusion, and the oxide is reduced at the metal-oxide interface. This facilitates solid-state sintering, as well as wetting of the underlying metal by sintering liquids, when these are present. The optimum magnesium concentration is approximately 0.1 to 1.0 wt pet, but this is dependent on the volume of oxide and, hence, the particle size, as well as the sintering conditions. Small particle-size fractions require proportionally more magnesium than large-size fractions do.

Experimental study of phase equilibria in the system PbO-ZnO-SiO2

An experimental study on the ternary system PbO-ZnO-SiO2, in air by high-temperature equilibration and quenching techniques followed by electron probe X-ray microanalysis was carried out as part of the wider research program on the six-component system PbO-ZnO-SiO2-CaO-FeO-Fe2O3, which combines experimental and thermodynamic computer modeling techniques to characterize zinc and lead industrial slags. Liquidus and solidus data were reported for all primary phase fields in the system PbO-ZnO-SiO2 in the temperature range 640 degrees C to 1400 degrees C (913 to 1673 K).

Influence of microstructure on the corrosion of diecast AZ91D

The corrosion of die cast AZ91D was studied and related to its microstructure. For comparison and to more fully understand the behaviour of die cast AZ91D, corrosion studies and microstructural examinations were also carried out using slowly solidified high purity AZ91, Mg-2%Al, Mg-9%Al, low purity magnesium and high purity magnesium. Corrosion was studied in 1N NaCl at pH 11 by (1) observing the corrosion morphology, (2) measuring electrochemical polarisation curves and (3) simultaneously measuring both the hydrogen evolution rate and the magnesium dissolution rate. The skin of die cast AZ91D showed better corrosion resistance than the interior. This is attributed to a combination of(1) a higher volume fraction of the beta phase, (2) a more continuous beta phase distribution around finer alpha grains, and (3) lower porosity in the skin layer than in the interior of the die casting. This study showed that the casting method can influence the corrosion performance by its influence on the alloy microstructure. (C) 1999 Elsevier Science Ltd. All rights reserved.

A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Tidal effects on sea water intrusion in unconfined aquifers

A variable-density groundwater model is used to analyse the effects of tidal fluctuations on sea-water intrusion in an unconfined aquifer. It is shown that the tidal activity forces the sea-water to intrude further inland and it also creates a thicker interface than would occur without tidal effects. Moreover, the configuration of the interface is radically changed when the tidal fluctuations are included. This is because of the dramatic changes in the flow pattern and velocity of the groundwater near the shoreline. For aquifer depths much larger than tidal amplitudes, the tidal fluctuation does not have much effect on how far the sea-water intrudes into the aquifer; nevertheless, a significant change in the configuration of concentration contours because of the effect of tidal fluctuations is observed. This change is more noticeable at the top of the aquifer, near the water table, than at the bottom of the aquifer, and is caused by the infiltration of salt water into the top of the aquifer at higher tidal levels. A flatter beach slope, therefore, intensifies this phenomenon. The interface configurations do not change noticeably over the course of a tidal cycle. Neglecting tidal fluctuation effects results in an inaccurate evaluation of the water table elevation at the land end of the aquifer, although no distinguishable difference is seen between the water tables near the shoreline. Where the landward boundary condition is a constant head, the effects of tidal fluctuations on sea-water intrusion are more pronounced than for cases where the landward boundary condition is a specified flux. Also it is shown that the effects of tidal fluctuations are more significant for a sloping beach than for a vertical shoreline and the salt water intrudes further inland for the sloping case. (C) 1999 Elsevier Science B.V. All rights reserved.

Modeling the fluid-flow-induced stress and collapse in a dendritic network

An analytical approach to the stress development in the coherent dendritic network during solidification is proposed. Under the assumption that stresses are developed in the network as a result of the friction resisting shrinkage-induced interdendritic fluid flow, the model predicts the stresses in the solid. The calculations reflect the expected effects of postponed dendrite coherency, slower solidification conditions, and variations of eutectic volume fraction and shrinkage. Comparing the calculated stresses to the measured shear strength of equiaxed mushy zones shows that it is possible for the stresses to exceed the strength, thereby resulting in reorientation or collapse of the dendritic network.

A wavelet visible difference predictor

In this paper, we describe a model of the human visual system (HVS) based on the wavelet transform. This model is largely based on a previously proposed model, but has a number of modifications that make it more amenable to potential integration into a wavelet based image compression scheme. These modifications include the use of a separable wavelet transform instead of the cortex transform, the application of a wavelet contrast sensitivity function (CSP), and a simplified definition of subband contrast that allows us to predict noise visibility directly from wavelet coefficients. Initially, we outline the luminance, frequency, and masking sensitivities of the HVS and discuss how these can be incorporated into the wavelet transform. We then outline a number of limitations of the wavelet transform as a model of the HVS, namely the lack of translational invariance and poor orientation sensitivity. In order to investigate the efficacy of this wavelet based model, a wavelet visible difference predictor (WVDP) is described. The WVDP is then used to predict visible differences between an original and compressed (or noisy) image. Results are presented to emphasize the limitations of commonly used measures of image quality and to demonstrate the performance of the WVDP, The paper concludes with suggestions on bow the WVDP can be used to determine a visually optimal quantization strategy for wavelet coefficients and produce a quantitative measure of image quality.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

H-1-NMR structural studies of a cystine-linked peptide containing residues 71-93 of transthyretin and effects of a Ser84 substitution implicated in familial amyloidotic polyneuropathy

The Ile-->Ser84 substitution in the thyroid hormone transport protein transthyretin is one of over 50 variations found to be associated with familial amyloid polyneuropathy, a hereditary type of lethal amyloidosis. Using a peptide analogue of the loop containing residue 84 in transthyretin, we have examined the putative local structural effects of this substitution using H-1-NMR spectroscopy. The peptide, containing residues 71-93 of transthyretin with its termini linked via a disulfide bond, was found to possess the same helix-turn motif as in the corresponding region of the crystallographically derived structure of transthyretin in 20% trifluoroethanol (TFE) solution. It therefore, represents a useful model with which to examine the effects of amyloidogenic substitutions. In a peptide analogue containing the Ile84-->Ser substitution it was found that the substitution does not greatly disrupt the overall three-dimensional structure, but leads to minor local differences at the turn in which residue 84 is involved. Coupling constant and NOE measurements indicate that the helix-turn motif is still present, but differences in chemical shifts and amide-exchange rates reflect a small distortion. This is in keeping with observations that several other mutant forms of transthyretin display similar subunit interactions and those that have been structurally analysed possess a near native structure. We propose that the Ser84 mutation induces only subtle perturbations to the transthyretin structure which predisposes the protein to amyloid formation.

A comparative study of the fluidized-bed coating of cylindrical metal surfaces with various thermoplastic polymer powders

This paper reports the results of an experimental investigation into the fluidized-bed coating of cylindrical metal specimens using two types of thermoplastic powders, Rilsan(R) PA11, a nylon-11 powder produced by Elf Atochem, France and Cotene(TM) 4612, a linear low density polyethylene powder produced by J.R Courtenay (New Zealand). The effects of dipping time, preheat temperature and particle size distribution on coating thickness and surface finish were investigated. Consistent trends in coating thickness growth with dipping time were obtained for both nylon-11 and polyethylene powders with increases in coating thickness with preheat temperature. For the same preheat temperature, the lower melting point of polyethylene results in thicker coatings compared to those of nylon-11. There is a negligible change in the coating thickness for sieved powders compared to that for unsieved powders. A pre-heat temperatures of between 240 degrees C and 300 degrees C is necessary to achieve an acceptable surface finish with both nylon-11 and polyethylene powders. To minimize errors in achieving the desired coating thickness, dipping times shorter than 2 s are not recommended. The use of graphs of coating thickness versus dipping time in combination with the coating surface roughness plots presented in this paper enable the optimal choice of pre-heat temperature and dipping time to achieve acceptable surface finish. (C) 1999 Elsevier Science S.A. All rights reserved.

Equiaxed solidification of Al-Si alloys

An experimental programme has been undertaken to determine which of the grain formation mechanisms of equiaxed crystals are dominant in the solidification of Al-Si foundry alloys. Small ingots were cast from alloys of varying silicon concentration with and without gauze barriers, using different types of mould materials and different mould preheats. The results show that two mechanisms of grain nucleation are operating. The first is a wall mechanism where crystals are nucleated either on or near the mould wall owing to thermal undercooling. The second is a constitutional supercooling mechanism where nucleants are activated in the constitutionally undercooled zone ahead of the advancing interface. As a consequence, the grain size decreases with increasing silicon content. However a transition in the growth mode occurs once a critical degree of constitutional undercooling is exceeded. This change in growth is accompanied by an increase in grain size. The transition point can be shifted with respect to solute content by changing the casting conditions, and a mechanism is proposed to explain this effect. MST/4109