Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.

On the solute coupling at the moving solid/liquid interface during equiaxed solidification

Integral mass conservation was widely accepted for the solute coupling to solve solute redistribution during equiaxed solidification so far. The present study revealed that the integral form was invalid for moving boundary problems as it could not represent the mass balance at the moving interface. Accordingly, differential mass conservation at the solid/liquid interface was used to solve solute diffusion for spherical geometry. The model was applied for hydrogen diffusion in solidification to validate that the hydrogen enrichment was significant and depended on the growth rate. (c) 2006 American Institute of Physics.

Fragmented coherence: Solid waste management in Colombo

Notwithstanding the increasingly fragmented organizational relationships within Colombo's urban governance system, the cooperative nature of stakeholder relationships lends a high level of coherence to the overall system. Since 1995, Colombo's solid waste management system has been characterized by the increased role of the private sector, community-based organizations and NGOs. Whilst the increasingly fragmented nature of this system exhibits some deeply ingrained problems, there are also a number of positives associated with the increased role of civil society actors and, in particular, the informal sector. Reforming regulatory frameworks so as to integrate some of the social norms that are integral to the lives of the majority of urban residents will contribute to regulatory frameworks being considerably more enforceable than is currently the case. Such reform requires that institutional and regulatory frameworks need to be flexible enough to adapt to the changing social, political and economic context. In the Colombo case, effective cooperation between public sector and civil society stakeholders illustrates that adaptive institutional arrangements grounded in pragmatism are feasible. The challenge that arises is to translate these institutional arrangements into adaptive regulatory frameworks - something that would require a significant mind shift on the part of planners and urban managers.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

A new relaxation method for roll forming problems

Finite element analysis (FEA) of nonlinear problems in solid mechanics is a time consuming process, but it can deal rigorously with the problems of both geometric, contact and material nonlinearity that occur in roll forming. The simulation time limits the application of nonlinear FEA to these problems in industrial practice, so that most applications of nonlinear FEA are in theoretical studies and engineering consulting or troubleshooting. Instead, quick methods based on a global assumption of the deformed shape have been used by the roll-forming industry. These approaches are of limited accuracy. This paper proposes a new form-finding method - a relaxation method to solve the nonlinear problem of predicting the deformed shape due to plastic deformation in roll forming. This method involves applying a small perturbation to each discrete node in order to update the local displacement field, while minimizing plastic work. This is iteratively applied to update the positions of all nodes. As the method assumes a local displacement field, the strain and stress components at each node are calculated explicitly. Continued perturbation of nodes leads to optimisation of the displacement field. Another important feature of this paper is a new approach to consideration of strain history. For a stable and continuous process such as rolling and roll forming, the strain history of a point is represented spatially by the states at a row of nodes leading in the direction of rolling to the current one. Therefore the increment of the strain components and the work-increment of a point can be found without moving the object forward. Using this method we can find the solution for rolling or roll forming in just one step. This method is expected to be faster than commercial finite element packages by eliminating repeated solution of large sets of simultaneous equations and the need to update boundary conditions that represent the rolls.

Ultramicroporous membranes for hydrogen separation

Fuel cell systems offer excellent efficiencies when compared to internal combustion engines, which result in reduced fuel consumption and greenhouse gas emissions. One of the areas requiring research for the success of fuel cell technology is the H2 fuel purification to reduce CO, which is a poison to fuel cells. Molecular sieve silica (MSS) membranes have a potential application in this area. In this work showed activated transport, a characteristic of ultramicroporous (dp

Thermal Cycling Stability of Silica Membranes for Gas Separation

Hydrogen is being seen as an alternative energy carrier to conventional hydrocarbons to reduce greenhouse gas emissions. High efficiency separation technologies to remove hydrogen from the greenhouse gas, carbon dioxide, are therefore in growing demand. Traditional thermodynamic separation systems utilise distillation, absorption and adsorption, but are limited in efficiency at compact scales. Molecular sieve silica (MSS) membranes can perform this separation as they have high permselectivity of hydrogen to carbon dioxide, but their stability under thermal cycling is not well reported. In this work we exposed a standard MSS membrane and a carbonised template MSS (CTMSS) membrane to thermal cycling from 100 to 450°C. The standard MSS and carbonised template CTMSS membranes both showed permselectivity of helium to nitrogen dropping from around 10 to 6 in the first set of cycles, remaining stable until the last test. The permselectivity drop was due to small micropore collapse, which occurred via structure movement during cycling. Simulating single stage membrane separation with a 50:50 molar feed of H2:CO2, H2 exiting the permeate stream would start at 79% and stabilise at 67%. Higher selectivity membranes showed less of a purity drop, indicating the margin at which to design a stable membrane separation unit for CO2 capture.

Molecular sieve silica membranes for H2/CO separation

Efficient separation of fuel gas (H2) from other gases in reformed gas mixtures is becoming increasingly important in the development of alternative energy systems. A highly efficient and new technology available for these separations is molecular sieve silica (MSS) membranes derived from tetraethyl-orthosilicate (TEOS). A permeation model is developed from an analogous electronic system and compared to transport theory to determine permeation, selectivity and apparent activation of energy based on experimental values. Experimental results for high quality membranes show single gas permselectivity peaking at 57 for H2/CO at 150°C with a H2 permeation of 5.14 x 10^-8 mol.m^-2.s^-1.Pa^-1. Higher permeance was also achieved, but at the expense of selectivity. This is the case for low quality membranes with peak H2 permeation at 1.78 x 10-7 mol.m-2.s-1.Pa-1 at 22°C and H2/CO permselectivity of 4.5. High quality membranes are characterised with positive apparent activation energy while the low quality membranes have negative values. The model had a good fit of r-squared of 0.99-1.00 using the experimental data.

Molecular Sieve Silica (MSS) Membranes for Gas Separation and Reaction Processes

Weakly branched silica films formed by the two-step sol-gel process allow for the formation of high selectivity membranes for gas separation. 29Si NMR and gas permeation showed that reduced crosslinking leads to He/CH4 selectivity improvement from 300 to 1000. Applied in membrane reactor for cyclohexane conversion to benzene, conversions were achieved at 14 fold higher than a conventional reactor at 250°C. Hydrothermal stability studies showed that carbon templating of silica is required for hydrothermally stable membranes. From our work it was shown that with correct application of chemistry, practical membrane systems can be built to suit gas separation (e. g. hydrogen fuel) and reactor systems.

Hydrothermal stable templated molecular sieve silica (TMSS) membranes for gas separation

In this work we compare the hydrothermal stability performance of a Templated Molecular Sieve Silica (TMSS) membrane against a standard, non-templated Molecular Sieve Silica (MSS) membrane. The tests were carried under dry and wet (steam) conditions for single gas (He, H2, CO and CO2) at 1-2 atm membrane pressure drop at 200oC. Single gas TMSS membrane H2, permeance and H2/CO permselectivity was found to be 2.05 x 10-8 mols.m-2.s-1.Pa-1 and 15, respectively. The MSS membrane showed similar selectivity, but increased overall flux. He permeance through membranes decayed at a rate of 4-5 x 10-10 mols.m-2.s-1.Pa-1 per day regardless of membrane ambience (dry or wet). Although H2/CO permselectivity of the TMSS membrane slightly improved from 15 to 18 after steam testing, the MSS membrane resulted in significant reduction from 16 to 8.3. In addition, membrane regeneration after more than 50 days resulted in the TMSS membrane reverting to its original permeation levels while no significant improvements were observed for the MSS membra ne. Results showed that the TMSS membrane had enhanced hydrothermal stability and regeneration ability.

Hydrothermal Stability of Modified Silica Membranes for Gas Separation

A new class of hybrid molecular sieve silica (MSS) membranes is developed and tested against standard and organic templated membranes. The hybrid membrane is synthesized by the standard sol-gel process, integrating a template (methyltriethoxysilane - MTES) and a C6 surfactant (triethylhexylammonium bromide) into the silica film matrix. After hydro treatment under a relative humidity of 96% for 50h, the hybrid membrane shows no changes in its gas separation capabilities or energy of mobility. The structural characteristics and integrity of the hybrid membrane are retained due to a high concentration of organophilic functional groups and alkoxides observed using 29 Si NMR. In contrast, the structural integrity of the membranes prepared with non-templated films deteriorated during the hydro treatment due to a large percentage of silanol groups (Si-OH) which react with water. The hybrid membranes underwent a decrease in the H2/CO2 selectivity of only 1% whereas for the non-templated membrane a 21% decrease was observed. The transport mechanism of the hybrid membranes is activated as permeation increased with temperature. The activation energy for the permeation of H2 is positive while negative for CO2. The H2 permeation obtained was 3x 10 -8 mol.m -2 .s -1 .Pa -1 and permselectivities for H2/CO2 and H2/N2 varied between 1-7 and 31-34, respectively.