Diagnostic goal driven modelling and simulation of multiscale process systems

A systematic goal-driven top-down modelling methodology is proposed that is capable of developing a multiscale model of a process system for given diagnostic purposes. The diagnostic goal-set and the symptoms are extracted from HAZOP analysis results, where the possible actions to be performed in a fault situation are also described. The multiscale dynamic model is realized in the form of a hierarchical coloured Petri net by using a novel substitution place-transition pair. Multiscale simulation that focuses automatically on the fault areas is used to predict the effect of the proposed preventive actions. The notions and procedures are illustrated on some simple case studies including a heat exchanger network and a more complex wet granulation process.

A reinvestigation of phase equilibria in the system Al2O3-SiO2-ZnO

The phase equilibria and liquidus temperatures in the binary SiO2-ZnO system and in the ternary Al2O3-SiO2-ZnO system at low Al2O3 concentrations have been experimentally determined using the equilibration and quenching technique followed by electron probe X-ray microanalysis. In the SiO2-ZnO system, two binary eutectics involving the congruently melting willemite (Zn2SiO4) were found at 1448 +/- 5 degrees C and 0.52 +/- 0.01 mole fraction ZnO and at 1502 +/- 5 degrees C and 0.71 +/- 0.01 mole fraction ZnO, respectively. The two ternary eutectics involving willemite previously reported in the Al2O3SiO2-ZnO system were found to be at 1315 +/- 5 degrees C and 1425 +/- 25 T, respectively. The compositions of the eutectics are 0.07, 0.52, and 0.41 and 0.05, 0.28, and 0.67 mole fraction Al2O3, SiO2, and ZnO, respectively. The results of the present investigation are significantly different from the results of previous studies.

Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface

We present a new version of non-local density functional theory (NL-DFT) adapted to description of vapor adsorption isotherms on amorphous materials like non-porous silica. The novel feature of this approach is that it accounts for the roughness of adsorbent surface. The solid–fluid interaction is described in the same framework as in the case of fluid–fluid interactions, using the Weeks–Chandler–Andersen (WCA) scheme and the Carnahan–Starling (CS) equation for attractive and repulsive parts of the Helmholtz free energy, respectively. Application to nitrogen and argon adsorption isotherms on non-porous silica LiChrospher Si-1000 at their boiling points, recently published by Jaroniec and co-workers, has shown an excellent correlative ability of our approach over the complete range of pressures, which suggests that the surface roughness is mostly the reason for the observed behavior of adsorption isotherms. From the analysis of these data, we found that in the case of nitrogen adsorption short-range interactions between oxygen atoms on the silica surface and quadrupole of nitrogen molecules play an important role. The approach presented in this paper may be further used in quantitative analysis of adsorption and desorption isotherms in cylindrical pores such as MCM-41 and carbon nanotubes.

Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory

Adsorption of argon and nitrogen at their respective boiling points in cylindrical pores of MCM-41 type silica-like adsorbents is studied by means of a non-local density functional theory (NLDFT), which is modified to deal with amorphous solids. By matching the theoretical results of the pore filling pressure versus pore diameter against the experimental data, we arrive at a conclusion that the adsorption branch (rather than desorption) corresponds to the true thermodynamic equilibrium. If this is accepted, we derive the optimal values for the solid–fluid molecular parameters for the system amorphous silica–Ar and amorphous silica–N2, and at the same time we could derive reliably the specific surface area of non-porous and mesoporous silica-like adsorbents, without a recourse to the BET method. This method is then logically extended to describe the local adsorption isotherms of argon and nitrogen in silica-like pores, which are then used as the bases (kernel) to determine the pore size distribution. We test this with a number of adsorption isotherms on the MCM-41 samples, and the results are quite realistic and in excellent agreement with the XRD results, justifying the approach adopted in this paper.

Impact of potential models on adsorption of linear molecules on carbon black

In this paper, we investigate the effects of potential models on the description of equilibria of linear molecules (ethylene and ethane) adsorption on graphitized thermal carbon black. GCMC simulation is used as a tool to give adsorption isotherms, isosteric heat of adsorption and the microscopic configurations of these molecules. At the heart of the GCMC are the potential models, describing fluid-fluid interaction and solid-fluid interaction. Here we studied the two potential models recently proposed in the literature, the UA-TraPPE and AUA4. Their impact in the description of adsorption behavior of pure components will be discussed. Mixtures of these components with nitrogen and argon are also studied. Nitrogen is modeled a two-site plus discrete charges while argon as a spherical particle. GCMC simulation is also used for generating simulation mixture isotherms. It is found that co-operation between species occurs when the surface is fractionally covered while competition is important when surface is fully loaded.

Capillary phenomena in the framework of the two-dimensional density functional theory

We present results of application of the density functional theory (DFT) to adsorption and desorption in finite and infinite cylindrical pores accounting for the density distribution in radial and axial directions. Capillary condensation via formation of bridges is considered using canonical and grand canonical versions of the 2D DFT. The potential barrier of nucleation is determined as a function of the bulk pressure and the pore diameter. In the framework of the conventional assumptions on intermolecular interactions both 1D and 2D DFT versions lead to the same results and confirm the classical scenario of condensation and evaporation: the condensation occurs at the vapor-like spinodal point, and the evaporation corresponds to the equilibrium transition pressure. The analysis of experimental data on argon and nitrogen adsorption on MCM-41 samples seems to not completely corroborate this scenario, with adsorption branch being better described by the equilibrium pressure - diameter dependence. This points to the necessity of the further development of basic representations on the hysteresis phenomena.

Adsorption of N-2, Ar, acetone, chloroform and acetone-chloroform mixture on carbon black in the framework of molecular layer structure theory (MLST)

Adsorption of pure nitrogen, argon, acetone, chloroform and acetone-chloroform mixture on graphitized thermal carbon black is considered at sub-critical conditions by means of molecular layer structure theory (MLST). In the present version of the MLST an adsorbed fluid is considered as a sequence of 2D molecular layers, whose Helmholtz free energies are obtained directly from the analysis of experimental adsorption isotherm of pure components. The interaction of the nearest layers is accounted for in the framework of mean field approximation. This approach allows quantitative correlating of experimental nitrogen and argon adsorption isotherm both in the monolayer region and in the range of multi-layer coverage up to 10 molecular layers. In the case of acetone and chloroform the approach also leads to excellent quantitative correlation of adsorption isotherms, while molecular approaches such as the non-local density functional theory (NLDFT) fail to describe those isotherms. We extend our new method to calculate the Helmholtz free energy of an adsorbed mixture using a simple mixing rule, and this allows us to predict mixture adsorption isotherms from pure component adsorption isotherms. The approach, which accounts for the difference in composition in different molecular layers, is tested against the experimental data of acetone-chloroform mixture (non-ideal mixture) adsorption on graphitized thermal carbon black at 50 degrees C. (C) 2005 Elsevier Ltd. All rights reserved.

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

Evaluation of 1-site and 5-site models of methane on its adsorption on graphite and in graphitic slit pores

In this paper, we evaluate the performance of the 1- and 5-site models of methane on the description of adsorption on graphite surfaces and in graphitic slit pores. These models have been known to perform well in the description of the fluid-phase behavior and vapor-liquid equilibria. Their performance in adsorption is evaluated in this work for nonporous graphitized thermal carbon black, and simulation results are compared with the experimental data of Avgul and Kiselev (Chemistry and Physics of Carbon; Dekker: New York, 1970; Vol. 6, p 1). On this nonporous surface, it is found that these models perform as well on isotherms at various temperatures as they do on the experimental isosteric heat for adsorption on a graphite surface. They are then tested for their performance in predicting the adsorption isotherms in graphitic slit pores, in which we would like to explore the effect of confinement on the molecule packing. Pore widths of 10 and 20 angstrom are chosen in this investigation, and we also study the effects of temperature by choosing 90.7, 113, and 273 K. The first two are for subcritical conditions, with 90.7 K being the triple point of methane and 113 K being its boiling point. The last temperature is chosen to represent the supercritical condition so that we can investigate the performance of these models at extremely high pressures. We have found that for the case of slit pores investigated in this paper, although the two models yield comparable pore densities (provided the accessible pore width is used in the calculation of pore density), the number of particles predicted by the I-site model is always greater than that predicted by the 5-site model, regardless of whether temperature is subcritical or supercritical. This is due to the packing effect in the confined space such that a methane molecule modeled as a spherical particle in the I-site model would pack better than the fused five-sphere model in the case of the 5-site model. Because the 5-site model better describes the liquid- and solid-phase behavior, we would argue that the packing density in small pores is better described with a more detailed 5-site model, and care should be exercised when using the 1-site model to study adsorption in small pores.

Batchelor's theory extended to elongated cylindrical or ellipsoidal particles

Extensions to Batchelor's theory have been derived to take into account different shaped particles while relating extensional viscosity enhancement to three parameters - shape, volume fraction and particle aspect ratio. The extended theory now allows calculation of the extensional viscosity enhancement, at a given volume fraction of particles, for either ellipsoidal or cylindrical particles. The formula improves the predictive capability of Batchelor's theory when compared with measurements found in the literature for different rod-like polymer solutions. (c) 2005 Elsevier B.V. All rights reserved.

Effects of potential models in the vapor-liquid equilibria and adsorption of simple gases on graphitized thermal carbon black

In this paper, we investigate the effects of various potential models in the description of vapor–liquid equilibria (VLE) and adsorption of simple gases on highly graphitized thermal carbon black. It is found that some potential models proposed in the literature are not suitable for the description of VLE (saturated gas and liquid densities and the vapor pressure with temperature). Simple gases, such as neon, argon, krypton, xenon, nitrogen, and methane are studied in this paper. To describe the isotherms on graphitized thermal carbon black correctly, the surface mediation damping factor introduced in our recent publication should be used to calculate correctly the fluid–fluid interaction energy between particles close to the surface. It is found that the damping constant for the noble gases family is linearly dependent on the polarizability, suggesting that the electric field of the graphite surface has a direct induction effect on the induced dipole of these molecules. As a result of this polarization by the graphite surface, the fluid–fluid interaction energy is reduced whenever two particles are near the surface. In the case of methane, we found that the damping constant is less than that of a noble gas having the similar polarizability, while in the case of nitrogen the damping factor is much greater and this could most likely be due to the quadrupolar nature of nitrogen.

Modeling nitrogen adsorption in spherical pores of siliceous materials by density functional theory

Adsorption of nitrogen in spherical pores of FDU-1 silica at 77 K is considered by means of a nonlocal density functional theory (NLDFT) accounting for a disordered structure of pore walls. Pore size distribution analysis of various FDU-1 samples subject to different temperatures of calcination revealed three distinct groups of pores. The principal group of pores is identified as ordered spherical mesopores connected with each other by smaller interconnecting pores and irregular micropores present in the mesopore walls. To account for the entrances (connecting pores) into spherical mesopores, a concept of solid mass distribution with respect to the apparent density was introduced. It is shown that the introduction of the aforementioned distribution was sufficient to quantitatively describe experimental adsorption isotherms over the entire range of relative pressures spanning six decades.

Polyethylene multiwalled carbon nanotube composites

Polyethylene (PE) multiwalled carbon nanotubes (MWCNTs) with weight fractions ranging from 0.1 to 10 wt% were prepared by melt blending using a mini-twin screw extruder. The morphology and degree of dispersion of the MWCNTs in the PE matrix at different length scales was investigated using scanning electron microscopy (SEM), transmission electron microscopy (TEM), atomic force microscopy (AFM) and wide-angle X-ray diffraction (WAXD). Both individual and agglomerations of MWCNTs were evident. An up-shift of 17 cm(-1) for the G band and the evolution of a shoulder to this peak were obtained in the Raman spectra of the nanocomposites, probably due to compressive forces exerted on the MWCNTs by PE chains and indicating intercalation of PE into the MWCNT bundles. The electrical conductivity and linear viscoelastic behaviour of these nanocomposites were investigated. A percolation threshold of about 7.5 wt% was obtained and the electrical conductivity of PE was increased significantly, by 16 orders of magnitude, from 10(-20) to 10(-4) S/cm. The storage modulus (G') versus frequency curves approached a plateau above the percolation threshold with the formation of an interconnected nanotube structure, indicative of 'pseudo-solid-like' behaviour. The ultimate tensile strength and elongation at break of the nanocomposites decreased with addition of MWCNTs. The diminution of mechanical proper-ties of the nanocomposites, though concomitant with a significant increase in electrical conductivity, implies the mechanism for mechanical reinforcement for PE/MWCNT composites is filler-matrix interfacial interactions and not filler percolation. The temperature of crystallisation (T.) and fraction of PE that was crystalline (F-c) were modified by incorporating MWCNTs. The thermal decomposition temperature of PE was enhanced by 20 K on addition of 10 wt% MWCNT. (c) 2005 Elsevier Ltd. All rights reserved.