Diffusion of Tracer in Altered Tonalite: Experiments and Simulations with Heterogeneous Distribution of Porosity

Numerical time-domain-diffusion simulations were used for studying the diffusion behavior of tracer molecules in rock matrix with homogeneous and heterogeneous porosity. As the heterogeneous sample in these simulations, a 3D tomographic image of altered tonalite was used, in which the mineral components and the pores resolved by X-ray microtomography were represented by their respective intragranular porosities determined previously by the ¹⁴C-PMMA method. The apparent diffusion coefficient of a tracer in altered tonalite was determined experimentally, and was then used in the simulations. In the altered tonalite analyzed, inclusion of heterogeneity in the porosity increased the diffusion coefficient by 16 %. Altered and pristine feldspar was the main mineral component in the sample (72 %), and it also provided the dominant contribution to tracer diffusion, explaining alone 52 % of the diffusion coefficient. The large pores resolved by microtomography (6 %) and altered and pristine mica (22 %) gave an equal contribution to the diffusion coefficient. The simulation method applied was also validated by comparing the results to both an analytical and a numerical solution to the diffusion equation in a homogenous medium. In addition, the method was compared to discrete-time random-walk simulations in the case of randomly placed overlapping spheres.     

Solutions to and Validation of Matrix-Diffusion Models

A model transport system is considered in which a pulse of tracer molecules is advected along a flow channel with porous walls. The advected tracer thus undergoes diffusion, matrix-diffusion, inside the walls, which affects the breakthrough curve of the tracer. Analytical solutions in the form of series expansions are derived for a number of situations which include a finite depth of the porous matrix, varying aperture of the flow channel, and longitudinal diffusion and Taylor dispersion of the tracer in the flow channel. Novel expansions for the Laplace transforms of the concentration in the channel facilitated closed-form expressions for the solutions. A rigorous result is also derived for the asymptotic form of the breakthrough curve for a finite depth of the porous matrix, which is very different from that for an infinite matrix. A specific experimental system was created for validation of matrix-diffusion modeling for a matrix of finite depth. A previously reported fracture-column experiment was also modeled. In both cases model solutions gave excellent fits to the measured breakthrough curves with very consistent values for the diffusion coefficients used as the fitting parameters. The matrix-diffusion modeling performed could thereby be validated.     

High‐molar‐mass polypropene with tunable elastic properties by hafnocene/borate catalysts

Elastic polypropene has gained growing industrial and academic interest as a thermoplastic elastomer. In this study, “rac”‐ and “meso”‐dimethylsilyl(3‐benzylindenyl)(2‐methylindenyl)hafnium dichloride complexes (Hfr and Hfm, respectively), activated with [NHMe₂Ph][B(C₆F₅)₄]/triisobutyl aluminum, were used in propene polymerization. Using these catalyst systems, we obtained polymers with high molar masses, up to 550 kg/mol, and moderate isotacticities between 34 and 52%. By varying the polymerization conditions, we could modify the polymer microstructure and molar mass. ¹³C NMR was used to calculate the polymer pentad sequence distributions. The crystalline parts of the polymers were analyzed with the differential scanning calorimetry successive self‐nucleation and annealing (SSA) technique. The SSA thermograms revealed that Hfr produced polypropene with a more uniform lamellar structure than Hfm. The mechanical properties were tested with dynamic mechanical analysis creep‐recovery tests. In the series, the polymers with the lowest isotacticities and therefore lowest crystallinities showed the best elastic properties. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 4743–4751, 2006     

Arbitrage without borrowing or short selling?

We show that a trader, who starts with no initial wealth and is not allowed to borrow money or short sell assets, is theoretically able to attain positive wealth by continuous trading, provided that she has perfect foresight of future asset prices, given by a continuous semimartingale. Such an arbitrage strategy can be constructed as a process of finite variation that satisfies a seemingly innocuous self-financing condition, formulated using a pathwise Riemann–Stieltjes integral. Our result exemplifies the potential intricacies of formulating economically meaningful self-financing conditions in continuous time, when one leaves the conventional arbitrage-free framework.     

A novel method to prepare water dispersible poly(benzimidazobenzophenanthroline) (BBL) by partial substitution of chain ends with poly(ethylene oxide)

Poly(benzimidazobenzophenanthroline) (BBL) was prepared according to literature method and modified with poly(ethylene oxide) in a one pot synthesis. After precipitation in aqueous sodium carbonate solution and subsequent purification, aqueous dispersions were prepared by ultrasonication. Particle sizes in the dispersions ranged from few tens of nanometers to several micrometers and most of the particles had sizes of 50–250 nm. Further studies indicated that the colloidal stability is a combined result of steric stabilization caused by excluded volume interactions of PEO chains on particle surface and electrostatic stabilization by the dissociated carboxylic acid groups on the particle surface. The product could be processed into uniform films 20–30 nm in thickness by spin coating onto gold-plated silicon substrates having aminethiol monolayer as the top most layer.     

Fluctuating lattice-Boltzmann model for complex fluids

We develop and test numerically a lattice-Boltzmann (LB) model for nonideal fluids that incorporates thermal fluctuations. The fluid model is a momentum-conserving thermostat, for which we demonstrate how the temperature can be made equal at all length scales present in the system by having noise both locally in the stress tensor and by shaking the whole system in accord with the local temperature. The validity of the model is extended to a broad range of sound velocities. Our model features a consistent coupling scheme between the fluid and solid molecular dynamics objects, allowing us to use the LB fluid as a heat bath for solutes evolving in time without external Langevin noise added to the solute. This property expands the applicability of LB models to dense, strongly correlated systems with thermal fluctuations and potentially nonideal equations of state. Tests on the fluid itself and on static and dynamic properties of a coarse-grained polymer chain under strong hydrodynamic interactions are used to benchmark the model. The model produces results for single-chain diffusion that are in quantitative agreement with theory.     

Structural principles of semiconducting Group 14 clathrate frameworks

We have performed a comprehensive theoretical investigation of the structural principles of semiconducting clathrate frameworks composed of the Group 14 elements carbon, silicon, germanium, and tin. We have investigated the basic clathrate frameworks, together with their polytypes, intergrowth clathrate frameworks, and extended frameworks based on larger icosahedral building blocks. Quantum chemical calculations with the PBE0 hybrid density functional method provided a clear overview of the structural trends and electronic properties among the various clathrate frameworks. In agreement with previous experimental and theoretical studies, the clathrate II framework proved to be the energetically most favorable, but novel hexagonal polytypes of clathrate II also proved to be energetically very favorable. In the case of silicon, several of the studied clathrate frameworks possess direct and wide band gaps. The band structure diagrams and simulated powder X-ray patterns of the studied frameworks are provided and systematic preliminary evaluation of guest-occupied frameworks is conducted to shed light on the characteristics of novel, experimentally feasible clathrate compositions.     

Stencils with isotropic discretization error for differential operators

We derive stencils, i.e., difference schemes, for differential operators for which the discretization error becomes isotropic in the lowest order. We treat the Laplacian, Bilaplacian (= biharmonic operator), and the gradient of the Laplacian both in two and three dimensions. For three dimensions ??(h²) results are given while for two dimensions both ??(h²) and ??(h⁴) results are presented. The results are also available in electronic form as a Mathematica file. It is shown that the extra computational cost of an isotropic stencil usually is less than 20%. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Nucleation, Growth, and Scaling in Slow Combustion

We study the nucleation and growth of flame fronts in slow combustion. This is modeled by a set of reaction-diffusion equations for the temperature field, coupled to a background of reactants and augmented by a term describing random temperature fluctuations for ignition. We establish connections between this model and the classical theories of nucleation and growth of droplets from a metastable phase. Our results are in good agreement with theoretical predictions.     

Asymptotic Lipschitz Regularity for Tug-of-War Games with Varying Probabilities

We prove an asymptotic Lipschitz estimate for value functions of tug-of-war games with varying probabilities defined in Ω ? ?ⁿ. The method of the proof is based on a game-theoretic idea to estimate the value of a related game defined in Ω ×Ω via couplings.