Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Diffusion analogy and shock capturing for solving aerodynamic equations

Improving the resolution of the shock is one of the most important subjects in computational aerodynamics. In this paper the behaviour of the solutions near the shock is discussed and the reason of the oscillation production is investigated heuristically. According to the differential approximation of the difference scheme the so-called diffusion analogy equation and the diffusion analogy coefficient are defined. Four methods for improving the resolution of the shock are presented using the concept of diffusion analogy.

Neutron and proton diffusion in fusion reactions for the synthesis of superheavy nuclei

The restriction of the one dimensional (1D) master equation (ME) with the mass number of the projectile-like fragment as a variable is studied, and a two-dimensional (2D) master equation with the neutron and proton numbers as independent variables is set up, and solved numerically. Our study showed that the 2D ME can describe the fusion process well in all projectile-target combinations. Therefore the possible channels to synthesize super-heavy nuclei can be studied correctly in wider possibilities. The available condition for employing 1D ME is pointed out.

Diffusion-limited kinetics modeling of one-step polyimide formation

A diffusion-limited kinetic model was developed to describe the imidization of one-step polythioetherimide formation based on an endgroup diffusion model. The changes of conversion and viscosity during the imidization were monitored with thermogravimetric analysis and dynamic stress rheometry, respectively. It was observed that the imidization rate began to decelerate after a fast early stage, whereas the viscosity in the system increased dramatically after a period of low value. Amic acid and imide formations concurrently take place in the one-step polyimide formation, but the formation of amic acid is much slower than that of imide and is the rate-limiting step of imidization. When a second-order kinetic model was used to describe the imidization, the effect of viscosity on the diffusion resistance of reactive groups needed to be included. In order to predict the change of viscosity during the imidization, the Lipshitz-Macosko model was modified and introduced into the diffusion-limited kinetic model by the Stokes-Einstein equation. The comparison of the modeled results with experimental data indicated that the diffusion-limited kinetic model and the modified Lipshitz-Macosko model were able to efficiently predict the changes of conversion and viscosity with temperature and time during the one-step polythioetherimide formation. (C) 2001 John Wiley & Sons, Inc.

Physical diffusion and electron-transfer dynamics of electroactive solutes in polymer electrolytes .2. Effect of the ionic size of supporting electrolytes

The diffusion coefficients(D-app) and the heterogeneous electron-transfer rate constants(k(s)) for ferrocene in MPEG/salt electrolytes were determined by using steady-state voltammetry. The temperature dependence of the two parameters obeys the Arrhenius equation. The effect of the ionic size of six supporting electrolytes on diffusion and electron transfer dynamics of ferrocene was discussed.

A study of the dependence of ferrocene diffusion on chain length in polymer electrolytes

The diffusion rates of ferrocene have been estimated in five kinds of poly(ethylene glycol) solution, containing the electrolyte LiClO4, by using non-steady-state chronoamperometry. The D-app of ferrocene increases with increasing temperature, and the dependency of D-app on temperature obeys the Williams-Landel-Ferry equation. The D-app of ferrocene decreases with increasing polymer chain length. Both the chain length and temperature dependence conform to a simple free volume model. A relation between current and polymer chain length is suggested at room and high temperatures.

Study on ferrocene derivatives diffusion dynamics in polymer electrolyte by solid-state voltammetry

The diffusion rates of seven ferrocene derivatives have been estimated in polyelectrolyte PEG . LiClO4 by using non-steady-state chronoamperometry. The D-app of ferrocene derivatives increases with temperature, and the dependency of D-app on temperature obeys the Arrhenius equation. The D-app of ferrocene derivatives decreases with increasing size of electroactive species. The Delta D-app values of D-T>Tm and D-Tdiffusion species. The dependency of D-app on the size of ferrocene at T < T-m is larger than that at T > T-m in the polyelectrolyte. On the other hand, the diffusion behaviour of ferrocene derivatives is qualitatively analyzed by using cyclic voltammetry. Copyright (C) 1996 Elsevier Science Ltd

STUDY ON MEASUREMENT OF DIFFUSION-COEFFICIENT OF LI ATOM IN LIQUID ALUMINUM

The tess potential-concentration curve was first applied to measure the concentration of an alloy. Attempt to use the V-c curve of Al-Li alloy in measuring the diffusion coeffictent of Li atom in liquid aluminium with anode chronopotentiometry at 720℃, was made and D_(Li/Al=4.94×10~(-5)cm~2·s~(-1) was obtained. The value is well consistent wish the theoretical value, D_(Li/Al)=4.85×10~(-5)cm~2·s~(-1) in terms of Stokes-Einstein equation.

Asymptotic periodicity of the Volterra equation with infinite delay

For some species, hereditary factors have great effects on their population evolution, which can be described by the well-known Volterra model. A model developed is investigated in this article, considering the seasonal variation of the environment, where the diffusive effect of the population is also considered. The main approaches employed here are the upper-lower solution method and the monotone iteration technique. The results show that whether the species dies out or not depends on the relations among the birth rate, the death rate, the competition rate, the diffusivity and the hereditary effects. The evolution of the population may show asymptotic periodicity, provided a certain condition is satisfied for the above factors. (c) 2006 Elsevier Ltd. All rights reserved.

The effect of hydrogen bonding on the diffusion of water in n-alkanes and n-alcohols measured with a novel single microdroplet method.

While the Stokes-Einstein (SE) equation predicts that the diffusion coefficient of a solute will be inversely proportional to the viscosity of the solvent, this relation is commonly known to fail for solutes, which are the same size or smaller than the solvent. Multiple researchers have reported that for small solutes, the diffusion coefficient is inversely proportional to the viscosity to a fractional power, and that solutes actually diffuse faster than SE predicts. For other solvent systems, attractive solute-solvent interactions, such as hydrogen bonding, are known to retard the diffusion of a solute. Some researchers have interpreted the slower diffusion due to hydrogen bonding as resulting from the effective diffusion of a larger complex of a solute and solvent molecules. We have developed and used a novel micropipette technique, which can form and hold a single microdroplet of water while it dissolves in a diffusion controlled environment into the solvent. This method has been used to examine the diffusion of water in both n-alkanes and n-alcohols. It was found that the polar solute water, diffusing in a solvent with which it cannot hydrogen bond, closely resembles small nonpolar solutes such as xenon and krypton diffusing in n-alkanes, with diffusion coefficients ranging from 12.5x10(-5) cm(2)/s for water in n-pentane to 1.15x10(-5) cm(2)/s for water in hexadecane. Diffusion coefficients were found to be inversely proportional to viscosity to a fractional power, and diffusion coefficients were faster than SE predicts. For water diffusing in a solvent (n-alcohols) with which it can hydrogen bond, diffusion coefficient values ranged from 1.75x10(-5) cm(2)/s in n-methanol to 0.364x10(-5) cm(2)/s in n-octanol, and diffusion was slower than an alkane of corresponding viscosity. We find no evidence for solute-solvent complex diffusion. Rather, it is possible that the small solute water may be retarded by relatively longer residence times (compared to non-H-bonding solvents) as it moves through the liquid.

Local kinetic interpretation of entropy production through reversed diffusion.

The time reversal of stochastic diffusion processes is revisited with emphasis on the physical meaning of the time-reversed drift and the noise prescription in the case of multiplicative noise. The local kinematics and mechanics of free diffusion are linked to the hydrodynamic description. These properties also provide an interpretation of the Pope-Ching formula for the steady-state probability density function along with a geometric interpretation of the fluctuation-dissipation relation. Finally, the statistics of the local entropy production rate of diffusion are discussed in the light of local diffusion properties, and a stochastic differential equation for entropy production is obtained using the Girsanov theorem for reversed diffusion. The results are illustrated for the Ornstein-Uhlenbeck process.

Modelling meso-scale diffusion processes in stochastic fluid bio-membranes

The space–time dynamics of rigid inhomogeneities (inclusions) free to move in a randomly fluctuating fluid bio-membrane is derived and numerically simulated as a function of the membrane shape changes. Both vertically placed (embedded) inclusions and horizontally placed (surface) inclusions are considered. The energetics of the membrane, as a two-dimensional (2D) meso-scale continuum sheet, is described by the Canham–Helfrich Hamiltonian, with the membrane height function treated as a stochastic process. The diffusion parameter of this process acts as the link coupling the membrane shape fluctuations to the kinematics of the inclusions. The latter is described via Ito stochastic differential equation. In addition to stochastic forces, the inclusions also experience membrane-induced deterministic forces. Our aim is to simulate the diffusion-driven aggregation of inclusions and show how the external inclusions arrive at the sites of the embedded inclusions. The model has potential use in such emerging fields as designing a targeted drug delivery system.

Fokker-Planck equation for a test-particle weakly coupled to a magnetized plasma

We consider the derivation of a kinetic equation for a charged test particle weakly interacting with an electrostatic plasma in thermal equilibrium, subject to a uniform external magnetic field. The Liouville equation leads to a generalized master equation to second order in the `weak' interaction; a Fokker-Planck-type equation then follows as a `Markovian' approximation. It is shown that such an equation does not preserve the positivity of the distribution function f(x,v;t). By applying techniques developed in the theory of open systems, a correct Fokker-Planck equation is derived. Explicit expressions for the diffusion and drift coefficients, depending on the magnetic field, are obtained.

Transport Properties Investigation of Aqueous Protic Ionic Liquid Solutions through Conductivity, Viscosity, and NMR Self-Diffusion Measurements

We present a study on the transport properties through conductivity (s), viscosity (?), and self-diffusion coefficient (D) measurements of two pure protic ionic liquids—pyrrolidinium hydrogen sulfate, [Pyrr][HSO4], and pyrrolidinium trifluoroacetate, [Pyrr][CF3COO]—and their mixtures with water over the whole composition range at 298.15 K and atmospheric pressure. Based on these experimental results, transport mobilities of ions have been then investigated in each case through the Stokes–Einstein equation. From this, the proton conduction in these PILs follows a combination of Grotthuss and vehicle-type mechanisms, which depends also on the water composition in solution. In each case, the displacement of the NMR peak attributed to the labile proton on the pyrrolidinium cation with the PILs concentration in aqueous solution indicates that this proton is located between the cation and the anion for a water weight fraction lower than 8%. In other words, for such compositions, it appears that this labile proton is not solvated by water molecules. However, for higher water content, the labile protons are in solution as H3O+. This water weight fraction appears to be the solvation limit of the H+ ions by water molecules in these two PILs solutions. However, [Pyrr][HSO4] and [Pyrr][CF3COO] PILs present opposed comportment in aqueous solution. In the case of [Pyrr][CF3COO], ?, s, D, and the attractive potential, Epot, between ions indicate clearly that the diffusion of each ion is similar. In other words, these ions are tightly bound together as ion pairs, reflecting in fact the importance of the hydrophobicity of the trifluoroacetate anion, whereas, in the case of the [Pyrr][HSO4], the strong H-bond between the HSO4– anion and water promotes a drastic change in the viscosity of the aqueous solution, as well as on the conductivity which is up to 187 mS·cm–1 for water weight fraction close to 60% at 298 K.

Comparative study on transport properties for LiFAP and LiPF6 in alkyl-carbonates as electrolytes through conductivity, viscosity and NMR self-diffusion measurements

We present in this work a comparative study on density and transport properties, such as the conductivity (sigma), viscosity (eta) and self-diffusion coefficients (D), for electrolytes based on the lithium hexafluorophosphate, LiPF6; or on the lithium tris(pentafluoroethane)-trifluorophosphate, LiFAP dissolved in a binary mixture of ethylene carbonate (EC) and dimethylcarbonate (DMC) (50:50 wt%). For each electrolyte, the temperature dependence on transport properties over a temperature range from 10 to 80 degrees C and 20 to 70 degrees C for viscosity and conductivity, respectively, exhibits a non-Arrhenius behavior. However, this dependence is correctly correlated by using the Vogel-Tamman-Fulcher (VTF) type fitting equation. In each case, the best-fit parameters, such as the pseudo activation energy and ideal glass transition temperature were then extracted. The self-diffusion coefficients (D) of the Li+ cation and PF6- or FAP(-) anions species, in each studied electrolyte, were then independently determined by observing Li-3, F-19 and P-31 nuclei with the pulsed-gradient spin-echo (PGSE) NMR technique over the same temperature range from 20 to 80 degrees C. Results show that even if the diffusion of the lithium cation is quite similar in both electrolytes, the anions diffusion differs notably. In the case of the LiPF6-based electrolyte, for example at T approximate to 75 degrees C (high temperature), the self-diffusion coefficients of Li+ cations in solution (D (Li+)approximate to 5 x 10(-19) m(2) s(-1)) is 1.6 times smaller than that of PF6- anions (D (PF6-) = 8.5 x 10(-19) m(2) s(-1)), whereas in the case of the LiFAP-based electrolyte, FAP(-) anions diffuse at same rate as the Li+ cations (D (FAP(-)) = 5 x 10(-1) m(2) s(-1)). Based on these experimental results, the transport mobility of ions were then investigated through Stokes-Einstein and Nernst-Einstein equations to determine the transport number of lithium t(Li)(+), effective radius of solvated Li+ and of PF6- and FAP(-) anions, and the degree of dissociation of these lithium salts in the selected EC/DMC (50:50 wt%) mixture over a the temperature range from 20 to 80 degrees C. This study demonstrates the conflicting nature of the requirements and the advantage of the well-balanced properties as ionic mobility and dissociation constant of the selected electrolytes. (C) 2013 Elsevier Ltd. All rights reserved.