Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Prediction of flame liftoff height of diffusion/partially premixed jet flames and modeling of mild combustion burners

In this article, a new flame extinction model based on the k/epsilon turbulence time scale concept is proposed to predict the flame liftoff heights over a wide range of coflow temperature and O-2 mass fraction of the coflow. The flame is assumed to be quenched, when the fluid time scale is less than the chemical time scale ( Da < 1). The chemical time scale is derived as a function of temperature, oxidizer mass fraction, fuel dilution, velocity of the jet and fuel type. The present extinction model has been tested for a variety of conditions: ( a) ambient coflow conditions ( 1 atm and 300 K) for propane, methane and hydrogen jet flames, ( b) highly preheated coflow, and ( c) high temperature and low oxidizer concentration coflow. Predicted flame liftoff heights of jet diffusion and partially premixed flames are in excellent agreement with the experimental data for all the simulated conditions and fuels. It is observed that flame stabilization occurs at a point near the stoichiometric mixture fraction surface, where the local flow velocity is equal to the local flame propagation speed. The present method is used to determine the chemical time scale for the conditions existing in the mild/ flameless combustion burners investigated by the authors earlier. This model has successfully predicted the initial premixing of the fuel with combustion products before the combustion reaction initiates. It has been inferred from these numerical simulations that fuel injection is followed by intense premixing with hot combustion products in the primary zone and combustion reaction follows further downstream. Reaction rate contours suggest that reaction takes place over a large volume and the magnitude of the combustion reaction is lower compared to the conventional combustion mode. The appearance of attached flames in the mild combustion burners at low thermal inputs is also predicted, which is due to lower average jet velocity and larger residence times in the near injection zone.

A semi-analytical solution for multilayer diffusion in a composite medium consisting of a large number of layers

Diffusion in a composite slab consisting of a large number of layers provides an ideal prototype problem for developing and analysing two-scale modelling approaches for heterogeneous media. Numerous analytical techniques have been proposed for solving the transient diffusion equation in a one-dimensional composite slab consisting of an arbitrary number of layers. Most of these approaches, however, require the solution of a complex transcendental equation arising from a matrix determinant for the eigenvalues that is difficult to solve numerically for a large number of layers. To overcome this issue, in this paper, we present a semi-analytical method based on the Laplace transform and an orthogonal eigenfunction expansion. The proposed approach uses eigenvalues local to each layer that can be obtained either explicitly, or by solving simple transcendental equations. The semi-analytical solution is applicable to both perfect and imperfect contact at the interfaces between adjacent layers and either Dirichlet, Neumann or Robin boundary conditions at the ends of the slab. The solution approach is verified for several test cases and is shown to work well for a large number of layers. The work is concluded with an application to macroscopic modelling where the solution of a fine-scale multilayered medium consisting of two hundred layers is compared against an “up-scaled” variant of the same problem involving only ten layers.

Theory of laminar flame propagation with non-normal diffusion.

A study has been made of the problem of steady, one-dimensional, laminar flame propagation in premixed gases, with the Lewis number differing from (and equal to) unity. Analytical solutions, using the method of matched asymptotic expansions, have been obtained for large activation energies. Numerical solutions have been obtained for a wide range of the reduced activation temperature parameter (n {geometrically equal to} E/RTb), and the Lewis number δ. The studies reveal that the flame speed eigenvalue is linear in Lewis number for first order and quadratic in Lewis number for second order reactions. For a quick determination of flame speeds, with reasonable accuracy, a simple rule, expressing the flame speed eigenvalue as a function of the Lewis number and the centroid of the reaction rate function, is proposed. Comparisons have been made with some of the earlier works, for both first and second order reactions.

The significance of brittle reaction layers in fusing of dental ceramics to titanium

This thesis comprises four intercomplementary parts that introduce new approaches to brittle reaction layers and mechanical compatibility of metalloceramic joints created when fusing dental ceramics to titanium. Several different methods including atomic layer deposition (ALD), sessile drop contact angle measurements, scanning acoustic microscopy (SAM), three-point bending (TPB, DIN 13 927 / ISO 9693), cross-section microscopy, scanning electron microscopy (SEM), and energy dispersive X-ray spectroscopy (EDS) were employed. The first part investigates the effects of TiO2 layer structure and thickness on the joint strength of the titanium-metalloceramic system. Samples with all tested TiO2 thicknesses displayed good ceramics adhesion to Ti, and uniform TPB results. The fracture mode was independent of oxide layer thickness and structure. Cracking occurred deeper inside titanium, in the oxygen-rich Ti[O]x solid solution surface layer. During dental ceramics firing TiO2 layers dissociate and joints become brittle with increased dissolution of oxygen into metallic Ti and consequent reduction in the metal plasticity. To accomplish an ideal metalloceramic joint this needs to be resolved. The second part introduces photoinduced superhydrophilicity of TiO2. Test samples with ALD deposited anatase TiO2 films were produced. Samples were irradiated with UV light to induce superhydrophilicity of the surfaces through a cascade leading to increased amount of surface hydroxyl groups. Superhydrophilicity (contact angle ~0˚) was achieved within 2 minutes of UV radiation. Partial recovery of the contact angle was observed during the first 10 minutes after UV exposure. Total recovery was not observed within 24h storage. Photoinduced ultrahydrophilicity can be used to enhance wettability of titanium surfaces, an important factor in dental ceramics veneering processes. The third part addresses interlayers designed to restrain oxygen dissolution into Ti during dental ceramics fusing. The main requirements for an ideal interlayer material are proposed. Based on these criteria and systematic exclusion of possible interlayer materials silver (Ag) interlayers were chosen. TPB results were significantly better in when 5 μm Ag interlayers were used compared to only Al2O3-blasted samples. In samples with these Ag interlayers multiple cracks occurred inside dental ceramics, none inside Ti structure. Ag interlayers of 5 μm on Al2O3-blasted samples can be efficiently used to retard formation of the brittle oxygen-rich Ti[O]x layer, thus enhancing metalloceramic joint integrity. The most brittle component in metalloceramic joints with 5 μm Ag interlayers was bulk dental ceramics instead of Ti[O]x. The fourth part investigates the importance of mechanical interlocking. According to the results, the significance of mechanical interlocking achieved by conventional surface treatments can be questioned as long as the formation of the brittle layers (mainly oxygen-rich Ti[O]x) cannot be sufficiently controlled. In summary in contrast to former impressions of thick titanium oxide layers this thesis clearly demonstrates diffusion of oxygen from sintering atmosphere and SiO2 to Ti structures during dental ceramics firing and the following formation of brittle Ti[O]x solid solution as the most important factors predisposing joints between Ti and SiO2-based dental ceramics to low strength. This among other predisposing factors such as residual stresses created by the coefficient of thermal expansion mismatch between dental ceramics and Ti frameworks can be avoided with Ag interlayers.

A study of the slow reaction in Al-5 wt.% Ag alloy by resistivity measurements

The slow reaction in an Al-5 wt.% Ag alloy has been investigated by resistivity measurements. The "slope change" method gave an activation energy of 1.25 eV for silver diffusion during the slow reaction. The existence of an excess concentration of vacancies in equilibrium with the dislocation loops seems to be responsible for the slow reaction. The presence of silver inhibits the nucleation of dislocation loops by holding up the quenched-in vacancies in solution. There is no indication of the presence of a third stage in the low-temperature ageing process of this alloy.

Unsteady mixed convection with double diffusion over a horizontal cylinder and a sphere within a porous medium

The combined effects of the permeability of the medium, magnetic field, buoyancy forces and dissipation on the unsteady mixed convection flow over a horizontal cylinder and a sphere embedded in a porous medium have been studied. The nonlinear coupled partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The skin friction, heat transfer and mass transfer increase with the permeability of the medium, magnetic field and buoyancy parameter. The heat and mass transfer continuously decrease with the stream-wise distance, whereas the skin friction increases from zero, attains a maximum and then decreases to zero. The skin friction, heat transfer and mass transfer are significantly affected by the free stream velocity distribution. The effect of dissipation parameter is found to be more pronounced on the heat transfer than on the skin friction and mass transfer

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

On Application Of The Method Of Kinetic Invariants To The Description Of Dissolution Accompanied By A Chemical Reaction

The dissolution, accompanied by chemical reaction, of monodisperse solid particles has been analysed. The resulting model, which accounts for the variation of mass transfer coefficient with the size of the dissolving particles, yields an approximate analytical form of a kinetic function. Rigorous numerical and approximate analytical solutions have been obtained for the governing system of nonlinear ordinary differential equations. The transient nature of the dissolution process as well as the accuracy of the analytical solution is brought out by the rigorous numerical solution. The analytical solution is fairly accurate for the major part of the range of operational times encountered in practice.

Study of the reaction between 1,2,3-trihydroxybenzene and 8-hydroxyquinoline in the solid state

1.2,3-Trihydroxybenzene (THB) reacts with 8-hydroxyquinoline (8HQ) in the solid state forming an orange-coloured charge transfer complex THB* (8HQ)(2). When the reaction was carried out in a petri dish, or when the vapours of 8HQ were allowed to react with solid THB (gravimetric study), the reaction product separated out as good quality, shiny single crystals. X-Ray diffraction studies on single crystals showed that they belong to the orthorhombic system with a = 15.408(1), b = 16.276(1), c = 7.825(1) Angstrom, Z = 4, D-x = 1.413 g cm(-3) and space group Pnaa. From the crystallographic evidence it has been found that the proton of the middle OH group of THB is transferred to the N atom of 8HQ. This accounts for the observed colour change. Kinetic studies on the solid state reaction showed that the 8HQ molecules diffuse towards THB, and the lateral diffusion occurs through surface migration, grain boundary diffusion and vapour phase diffusion. Gravimetric studies of the reaction between solid THB and 8HQ vapour showed that the diffusion of 8HQ molecules into the crystal lattice of THB has a higher energy of activation than that observed when the reactants are in contact. The nature of the crystal packing in the reaction product indicates diffusion of 8HQ molecules into the crystal lattice of THB along the c-axis, to occupy the cavities present between the THB molecules in the unit cell.

Shear-induced enhancement of self-diffusion in interacting colloidal suspensions

We set up the generalized Langevin equations describing coupled single-particle and collective motion in a suspension of interacting colloidal particles in a shear how and use these to show that the measured self-diffusion coefficients in these systems should be strongly dependent on shear rate epsilon. Three regimes are found: (i) an initial const+epsilon(.2), followed by (ii) a large regime of epsilon(.1/2) behavior, crossing over to an asymptotic power-law approach (iii) D-o - const x epsilon(.-1/2) to the Stokes-Einstein value D-o. The shear dependence is isotropic up to very large shear rates and increases with the interparticle interaction strength. Our results provide a straightforward explanation of recent experiments and simulations on sheared colloids.

Some Basic Aspects of Reaction Engineering of Precipitation Processes

Analysis of precipitation reactions is extremely important in the technology of production of fine particles from the liquid phase. The control of composition and particle size in precipitation processes requires careful analysis of the several reactions that comprise the precipitation system. Since precipitation systems involve several, rapid ionic dissociation reactions among other slower ones, the faster reactions may be assumed to be nearly at equilibrium. However, the elimination of species, and the consequent reduction of the system of equations, is an aspect of analysis fraught with the possibility of subtle errors related to the violation of conservation principles. This paper shows how such errors may be avoided systematically by relying on the methods of linear algebra. Applications are demonstrated by analyzing the reactions leading to the precipitation of calcium carbonate in a stirred tank reactor as well as in a single emulsion drop. Sample calculations show that supersaturation dynamics can assume forms that can lead to subsequent dissolution of particles that have once been precipitated.

Sandwich propellant combustion: Modeling and experimental comparison

This paper reports reacting fluid dynamics calculations for an ammonium percholrate binder sandwich and extracts experimentally observed features including surface profiles and maximum regression rates as a function of pressure and binder thickness. These studies have been carried out by solving the two-dimensional unsteady Navier-Stokes equations with energy and species conservation equations and a kinetic model of three reaction steps (ammonium perchlorate decomposition flame, primary diffusion flame, and final diffusion flame) in the gas phase. The unsteady two-dimensional conduction equation is solved in the condensed phase. The regressing surface is unsteady and two dimensional. Computations have been carried out for a binder thickness range of 25-125 mum and a pressure range of 1.4 to 6.9 MPa. Good comparisons at several levels of detail are used to demonstrate the need for condensed-phase two-dimensional unsteady conduction and three-step gas-phase reactions. The choice of kinetic and thermodynamic parameters is crucial to good comparison with experiments. The choice of activation energy parameters for ammonium percholrate combustion has been made with stability of combustion in addition to experimentally determined values reported in literature. The choice of gas-phase parameters for the diffusion flames are made considering that (a) primary diffusion flame affects the low-pressure behavior and (b) final diffusion flame affects high-pressure behavior. The predictions include the low-pressure deflagration limit of the sandwich apart from others noted above. Finally, this study demonstrates the possibility of making meaningful comparisons with experimental observations on sandwich propellant combustion.

A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Unsteady free convection flow over an infinite vertical porous plate due to the combined effects of thermal and mass diffusion, magnetic field and Hall currents

The unsteady free convection flow over an infinite vertical porous plate, which moves with time-dependent velocity in an ambient fluid, has been studied. The effects of the magnetic field and Hall current are included in the analysis. The buoyancy forces arise due to both the thermal and mass diffusion. The partial differential equations governing the flow have been solved numerically using both the implicit finite difference scheme and the difference-differential method. For the steady case, analytical solutions have also been obtained. The effect of time variation on the skin friction, heat transfer and mass transfer is very significant. Suction increases the skin friction coefficient in the primary flow, and also the Nusselt and Sherwood numbers, but the skin friction coefficient in the secondary flow is reduced. The effect of injection is opposite to that of suction. The buoyancy force, injection and the Hall parameter induce an overshoot in the velocity profiles in the primary flow which changes the velocity gradient from a negative to a positive value, but the magnetic field and suction reduce this velocity overshoot.