An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Modeling of gas adsorption equilibrium over a wide range of pressure: A thermodynamic approach based on equation of state

A thermodynamic approach based on the Bender equation of state is suggested for the analysis of supercritical gas adsorption on activated carbons at high pressure. The approach accounts for the equality of the chemical potential in the adsorbed phase and that in the corresponding bulk phase and the distribution of elements of the adsorption volume (EAV) over the potential energy for gas-solid interaction. This scheme is extended to subcritical fluid adsorption and takes into account the phase transition in EAV The method is adapted to gravimetric measurements of mass excess adsorption and has been applied to the adsorption of argon, nitrogen, methane, ethane, carbon dioxide, and helium on activated carbon Norit R I in the temperature range from 25 to 70 C. The distribution function of adsorption volume elements over potentials exhibits overlapping peaks and is consistently reproduced for different gases. It was found that the distribution function changes weakly with temperature, which was confirmed by its comparison with the distribution function obtained by the same method using nitrogen adsorption isotherm at 77 K. It was shown that parameters such as pore volume and skeleton density can be determined directly from adsorption measurements, while the conventional approach of helium expansion at room temperature can lead to erroneous results due to the adsorption of helium in small pores of activated carbon. The approach is a convenient tool for analysis and correlation of excess adsorption isotherms over a wide range of pressure and temperature. This approach can be readily extended to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

Oxide ionic conductivity and microstructures of Sm- or La-doped CeO2-based systems

The concept of crystallographic index termed the effective index is suggested and applied to the design of ceria (CeO2)-based electrolytes to maximize oxide ionic conductivity. The suggested index considers the fluorite structure, and combines the expected oxygen vacancy level with the ionic radius mismatch between host and dopant cations. Using this approach, oxide ionic conductivity of Sm- or La-doped CeO2-based system has been optimized and tested under operating conditions of a solid oxide fuel cell. In the observation of microstructure in atomic scale, both Sm-doped CeO2 and La-doped CeO2 electrolytes had large micro-domains over 10 nm in the lattice. On the other hand, Sm or La and alkaline earth co-doped CeO2-based electrolytes with high effective index had small micro-domains around 1-3 nm in the microstructure. The large micro-domain would prevent oxide ion from passing through the lattice. Therefore, it is concluded that the improvement of ionic conductivity is reflected in changes of microstructure in atomic scale. (C) 2002 Elsevier Science B.V. All rights reserved.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Use of confocal microscopy to follow the development of penetrative hyphae during growth of Rhizopus oligosporus in an artificial solid-state fermentation system

Two methods were compared for determining the concentration of penetrative biomass during growth of Rhizopus oligosporus on an artificial solid substrate consisting of an inert gel and starch as the sole source of carbon and energy. The first method was based on the use of a hand microtome to make sections of approximately 0.2- to 0.4-mm thickness parallel to the substrate surface and the determination of the glucosamine content in each slice. Use of glucosamine measurements to estimate biomass concentrations was shown to be problematic due to the large variations in glucosamine content with mycelial age. The second method was a novel method based on the use of confocal scanning laser microscopy to estimate the fractional volume occupied by the biomass. Although it is not simple to translate fractional volumes into dry weights of hyphae due to the lack of experimentally determined conversion factors, measurement of the fractional volumes in themselves is useful for characterizing fungal penetration into the substrate. Growth of penetrative biomass in the artificial model substrate showed two forms of growth with an indistinct mass in the region close to the substrate surface and a few hyphae penetrating perpendicularly to the surface in regions further away from the substrate surface. The biomass profiles against depth obtained from the confocal microscopy showed two linear regions on log-linear plots, which are possibly related to different oxygen availability at different depths within the substrate. Confocal microscopy has the potential to be a powerful tool in the investigation of fungal growth mechanisms in solid-state fermentation. (C) 2003 Wiley Periodicals, Inc.

Determination of viscosity of some Australian honeys based on composition

Rheological properties of four unprocessed unifloral Australian honeys (heath, tea tree, yapunya, and yellow box) and an artificial honey were analysed at 20degreesC. A model previously used to describe viscosity data of various sugar and sugar mixtures was used to describe the concentration dependence of the viscosity of honey samples with varying moisture contents. The model successfully described the sugar concentration dependence of the unadulterated and medium moisture (70-85% solids) range honey samples.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.

Aproximação numérica – analítica para a modelagem da conversão termoquímica de combustíveis sólidos

Um algoritmo numérico foi criado para apresentar a solução da conversão termoquímica de um combustível sólido. O mesmo foi criado de forma a ser flexível e dependente do mecanismo de reação a ser representado. Para tanto, um sistema das equações características desse tipo de problema foi resolvido através de um método iterativo unido a matemática simbólica. Em função de não linearidades nas equações e por se tratar de pequenas partículas, será aplicado o método de Newton para reduzir o sistema de equações diferenciais parciais (EDP’s) para um sistema de equações diferenciais ordinárias (EDO’s). Tal processo redução é baseado na união desse método iterativo à diferenciação numérica, pois consegue incorporar nas EDO’s resultantes funções analíticas. O modelo reduzido será solucionado numericamente usando-se a técnica do gradiente bi-conjugado (BCG). Tal modelo promete ter taxa de convergência alta, se utilizando de um número baixo de iterações, além de apresentar alta velocidade na apresentação das soluções do novo sistema linear gerado. Além disso, o algoritmo se mostra independente do tamanho da malha constituidora. Para a validação, a massa normalizada será calculada e comparada com valores experimentais de termogravimetria encontrados na literatura, , e um teste com um mecanismo simplificado de reação será realizado.

Effect of fiber alignment in gelled electrospun membranes based on poly(vinylidene fluoride) for battery applications

Battery separators based on electrospun membranes of poly(vinylidene fluoride) (PVDF) have been prepared in order to study the effect of fiber alignment on the performance and characteristics of the membrane. The prepared membranes show an average fiber diameter of 272 nm and a degree of porosity of 87 %. The gel polymer electrolytes are prepared by soaking the membranes in the electrolyte solution. The alignment of the fibers improves the mechanical properties for the electrospun membranes. Further, the microstructure of the membrane also plays an important role in the ionic conductivity, being higher for the random electrospun membrane due to the lower tortuosity value. Independently of the microstructure, both membranes show good electrochemical stability up to 5.0 V versus Li/Li+. These results show that electrospun membranes based on PVDF are appropriate for battery separators in lithium-ion battery applications, the random membranes showing a better overall performance.

A mobile health care rule-based system

The relation between patient and physician in most modern Health Care Sys- tems is sparse, limited in time and very in exible. On the other hand, and in contradiction with several recent studies, most physicians do not rely their patient diagnostics evaluations on intertwined psychological and social nature factors. Facing these problems and trying to improve the patient/physician relation we present a mobile health care solution to im- prove the interaction between the physician and his patients. The solution serves not only as a privileged mean of communication between physicians and patients but also as an evolutionary intelligent platform delivering a mobile rule based system.

A creative solution for rapid and inexpensive model-making in product design processes

Rapid prototyping (RP) is an approach for automatically building a physical object through solid freeform fabrication. Nowadays, RP has become a vital aspect of most product development processes, due to the significant competitive advantages it offers compared to traditional manual model making. Even in academic environments, it is important to be able to quickly create accurate physical representations of concept solutions. Some of these can be used for simple visual validation, while others can be employed for ergonomic assessment by potential users or even for physical testing. However, the cost of traditional RP methods prevents their use in most academic environments on a regular basis, and even for very preliminary prototypes in many small companies. That results in delaying the first physical prototypes to later stages, or creating very rough mock-ups which are not as useful as they could be. In this paper we propose an approach for rapid and inexpensive model-making, which was developed in an academic context, and which can be employed for a variety of objects.

Estudo da lixiviação de um concentrado de zinco

Neste trabalho estudou-se a lixiviação em meio sulfúrico do zinco e outros metais de valor de um concentrado zinco, tendo-se realizado ensaios de lixiviação à pressão atmosférica e em autoclave. Nos estudos de lixiviação utilizou-se o ião férrico (sulfato férrico) como agente oxidante e avaliaram-se os efeitos de diversas variáveis como a razão sólido/líquido, concentração do ião Fe (III), temperatura, a pressão de oxigénio e a presença de enxofre elementar na eficiência da lixiviação. Os ensaios de lixiviação em autoclave sob pressão de oxigénio foram realizados para verificar o efeito da manutenção da quantidade de Fe (III) na lixívia, por oxidação do Fe(II) com oxigénio. Os resultados obtidos mostraram que à pressão atmosférica para uma razão sólido/líquido de 5% foi possível lixiviar no máximo 59% de zinco e 22% de cobre com solução de 0,25 M de Fe2(SO4)3 e 0,50 M de H2SO4 em 2 horas a 60ºC e com uma razão sólido/líquido de 5% foi possível lixiviar no máximo 65% de zinco e 23% de cobre com uma solução de 0,5 M de Fe2(SO4)3 e 0,25 M de H2SO4 em 2 horas a 80ºC. Efectuar a lixiviação do concentrado de zinco sobre pressão de oxigénio permitiu aumentar a cinética da reacção de lixiviação, tendo sido possível lixiviar 97% de zinco e 48% do cobre em 2 horas de lixiviação com uma solução de 0,25 M Fe2(SO4)3 e 0,5 M H2SO4 a 95 ºC e a 6 bar de pressão de oxigénio (à entrada do reactor) com uma razão sólido/líquido de 5%. Utilizando razão sólido/líquido de 10 % foi possível lixiviar 93% de zinco e 54% do cobre com uma solução de 0,50 M Fe2(SO4)3 e 1,25 M H2SO4 a 95 ºC e a 6 bar de pressão de oxigénio, e para uma razão sólido/líquido de 20 % foi possível lixiviar 84% de zinco e 39% do cobre com uma solução de 0,11 M Fe2(SO4)3 e 2,00 M H2SO4 a 95 ºC e a 10 bar de pressão. As análises de difracção de Raios X efectuados aos resíduos de lixiviação revelaram que o enxofre era maioritariamente oxidado a enxofre elementar. Assim, para um dos ensaios de lixiviação em autoclave, verificou-se que a remoção com tetracloreto de carbono do enxofre elementar formado num primeiro andar de lixiviação (s/l=20%, 0,11 M Fe2(SO4)3 e 2,00 M H2SO4 a 95 ºC e a 10 bar de pressão) permitia aumentar a percentagem de zinco no segundo andar de 42 para 68%. Por último, o estudo do efeito da temperatura permitiu calcular como base nas velocidades iniciais do zinco a energia de activação para a lixiviação do zinco que foi de 39 ± 1.40 kJ/mol para a lixiviação em autoclave e de 38 ± 1.40 kJ/mol para a lixiviação à pressão atmosférica, o que é indicativo do controlo reaccional.