The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

On an indirect integral equation approach for stationary heat transfer in semi-infinite layered domains in R3 with cavities

Regions containing internal boundaries such as composite materials arise in many applications.We consider a situation of a layered domain in IR3 containing a nite number of bounded cavities. The model is stationary heat transfer given by the Laplace equation with piecewise constant conductivity. The heat ux (a Neumann condition) is imposed on the bottom of the layered region and various boundary conditions are imposed on the cavities. The usual transmission (interface) conditions are satised at the interface layer, that is continuity of the solution and its normal derivative. To eciently calculate the stationary temperature eld in the semi-innite region, we employ a Green's matrix technique and reduce the problem to boundary integral equations (weakly singular) over the bounded surfaces of the cavities. For the numerical solution of these integral equations, we use Wienert's approach [20]. Assuming that each cavity is homeomorphic with the unit sphere, a fully discrete projection method with super-algebraic convergence order is proposed. A proof of an error estimate for the approximation is given as well. Numerical examples are presented that further highlights the eciency and accuracy of the proposed method.

Travelling waves in a complex Ginzburg-Landau equation with time-delay feedback

One of the simplest ways to create nonlinear oscillations is the Hopf bifurcation. The spatiotemporal dynamics observed in an extended medium with diffusion (e.g., a chemical reaction) undergoing this bifurcation is governed by the complex Ginzburg-Landau equation, one of the best-studied generic models for pattern formation, where besides uniform oscillations, spiral waves, coherent structures and turbulence are found. The presence of time delay terms in this equation changes the pattern formation scenario, and different kind of travelling waves have been reported. In particular, we study the complex Ginzburg-Landau equation that contains local and global time-delay feedback terms. We focus our attention on plane wave solutions in this model. The first novel result is the derivation of the plane wave solution in the presence of time-delay feedback with global and local contributions. The second and more important result of this study consists of a linear stability analysis of plane waves in that model. Evaluation of the eigenvalue equation does not show stabilisation of plane waves for the parameters studied. We discuss these results and compare to results of other models.

Integration of the Manakov-PMD Equation with Precomputed M(w) Matrices for PMD Simulation

A novel direct integration technique of the Manakov-PMD equation for the simulation of polarisation mode dispersion (PMD) in optical communication systems is demonstrated and shown to be numerically as efficient as the commonly used coarse-step method. The main advantage of using a direct integration of the Manakov-PMD equation over the coarse-step method is a higher accuracy of the PMD model. The new algorithm uses precomputed M(w) matrices to increase the computational speed compared to a full integration without loss of accuracy. The simulation results for the probability distribution function (PDF) of the differential group delay (DGD) and the autocorrelation function (ACF) of the polarisation dispersion vector for varying numbers of precomputed M(w) matrices are compared to analytical models and results from the coarse-step method. It is shown that the coarse-step method achieves a significantly inferior reproduction of the statistical properties of PMD in optical fibres compared to a direct integration of the Manakov-PMD equation.

Application of the manakov-PMD equation to a computational investigation of low-PMD fibres

The phenomenon of low-PMD fibres is examined through numerical simulations. Instead of the coarse-step method we are using an algorithm developed through the Manakov-PMD equation. With the integration of the Manakov-PMD equation we have access to the fibre spin which relates to the orientation of the birefringence. The simulation results produced correspond to the behaviour of a low-PMD spun fibre. Furthermore we provide an analytical approximation compared to the numerical data. © 2005 Optical Society of America.

Numerical implementation of the Manakov-PMD equation with precomputed M(Ω) matrices

The Manakov-PMD equation can be integrated with the same numerical efficiency as the coarse-step method by using precomputed M(Ω) matrices, which entirely avoids the somewhat ad-hoc rescaling of coefficients necessary in the coarse-step method.

Equation of state for water in the small compressibility region

The equation of state for dense fluids has been derived within the framework of the Sutherland and Katz potential models. The equation quantitatively agrees with experimental data on the isothermal compression of water under extrapolation into the high pressure region. It establishes an explicit relationship between the thermodynamic experimental data and the effective parameters of the molecular potential.

The creep behaviour of pressure diecast zinc-aluminium based alloys

The creep behaviour of three pressure diecast commercial zinc-aluminium based alloys: Mazak 3, corresponding to BS 1004A, and the new alloys ZA.8 and ZA.27 with a series of alloys with compositions ranging from 0% to 30% aluminium was investigated. The total creep elongation of commercial alloys was shown to be well correlated using an empirical equation. Based on this a parametrical relationship was derived which allowed the total creep extension to be related to the applied stress, the temperature and the time of test, so that a quantitative assessment of creep of the alloys could be made under different conditions. Deviation from the normal creep kinetics occurred in alloys ZA.8 and ZA.27 at very low stresses, 150°C, due to structural coarsening combined with partial transformation of ε -phase into T' phase. The extent of primary creep was found to increase with aluminium content, but secondary creep rates decreased in the order Mazak 3, ZA.8 and ZA.27. Thus, based on the above equation, ZA.8 was found to have a substantially better total creep resistance than ZA.27, which in turn was marginally better than Mazak 3 for strains higher than 0.5%, but inferior for smaller strains, due to its higher primary creep extension. The superior creep resistance of ZA.8 was found to be due to the presence of strictly-orientated, thin plate-like precipitates of ε(CuZn4) phase in the zinc matrix of the eutectic and the lamellarly decomposed β phase, in which the precipitation morphology and orientation of ε in the zinc matrix was determined. Over broad ranges of temperature and stresses, the stress exponents and activation energies for creep were found to be consistent with some proposed creep rate mechanisms; i.e. viscous glide for Mazak 3, dislocation climb over second phase particles for ZA.8 and dislocation climb for ZA.27, controlled by diffusion in the zinc-rich phase. The morphology of aluminium and copper-rich precipitates formed from the solid solution of zinc was clearly revealed. The former were found to further increase the creep rate of inherently low creep resistant zinc, but the latter contributed significantly to the creep resistance. Excess copper in the composition, however, was not beneficial in improving the creep resistance. Decomposition of β in copper-containing alloys was found to be through a metastable Zn-Al phase which is strongly stabilised by copper, and the final products of the decomposition had a profound effect on the creep strength of the alloys. The poor creep resistance of alloy ZA.27 was due to the presence of particulate products derived from decomposed β-phase and a large volume of fine, equiaxed products of continuously decomposed α-dendrites.

Osmotic dehydration in plant tissues

The primary aim of the thesis is to provide a comprehensive investigation of the osmotic dehydration processes in plant tissue. Effort has been concentrated on the modelling for simulating the processes. Two mathematical models for simulating the mass transfer during osmotic dehydration processes in plant tissues are developed and verified using existing experimental data. Both models are based on the mechanism of diffusion and convection of any mobile material that can transport in plant tissues. The mass balance equation for the transport of each constituent is established separately for intracellular and extra-cellular volumes with taking into account the mass transfer across the cell membrane the intracellular and extra-cellular volumes and the shrinkage of the whole tissue. The contribution from turgor pressure is considered in both models. Model two uses Darcy’s law to build the relation between shrinkage velocity and hydrostatic pressure in each volume because the plant tissue can be considered as the porous medium. Moreover, it has been extended to solve the multi-dimensional problems. A lot of efforts have been made to the parameter study and the sensitivity analyses. The parameters investigated including the concentration of the osmotic solution, diffusion coefficient, permeability of the cell membrane, elastic modulus of the cell wall, critical cell volume etc. The models allow us to quantitatively simulate the time evolution of intracellular and extra-cellular volumes as well as the time evolution of concentrations in each cross-section.

The creep properties of a series of zinc-rich zinc-aluminium alloys

The compressive creep behaviour of six sand cast zinc-rich alloys: No3 and No5, corresponding to BS 1004A and BS 1004B, respectively, alloy No2, ILZRO,.16 and two newer alloys ACuZinc5 and ACuZinc10 was investigated. The total creep contraction of the alloys was found to be well correlated using an empirical equation. On the basis of this equation, a parametrical relationship was derived which allowed the total creep contraction to be related to the applied stress, the temperature and the time of test, so that a quantitative assessment of compressive creep of the alloys could be made under different testing conditions. The primary creep and secondary creep rates were found for the alloys at different temperatures and stresses. Generally, the primary creep contraction was found to increase with copper content, whereas secondary creep rates decreased in the order No3, ACuZinc10, ACuZinc5 and No2. ILZRO.16 was tested only at the highest stress and two higher temperatures. The results showed that ILZRO.16 had higher creep resistance than all the other alloys. Thus, based on the above empirical equation, alloy No2 was found to have a substantially better total creep resistance than alloys No3 and No5, and slightly better than ACuZinc5 and ACuZinc10 for strains up to 1%. Both ACuZinc alloys had higher creep strength than commercial alloys No3 and No5. Alloy No5 had much higher creep resistance than alloy No3 under all conditions. The superior creep resistance of alloy No2 was considered to be due to the presence of small precipitates of -phase in the zinc matrix and a regular eutectic morphology. The stress exponents and activation energies for creep under different testing conditions were found to be consistent with some established creep-controlling mechanisms; i.e. dislocation climb for alloy No3, dislocation climb over second phase particles for alloys No5, No2, ACuZinc10, controlled by lattice diffusion in the zinc-rich phase. The lower creep resistance of alloy No3 was mainly due to the lower creep strength of copper-free primary particles having greater volume than eutectic in the microstructure. Alloys No5, ACuZinc5 and ACuZinc10 showed much better creep resistance than alloy No3, based on the precipitation-hardening due to the presence of small -phase precipitates. The primary dendrites in both ACuZinc alloys however were not of much benefit in improving the creep resistance of the alloys.

Random input problem for the nonlinear Schrödinger equation

We consider the random input problem for a nonlinear system modeled by the integrable one-dimensional self-focusing nonlinear Schrödinger equation (NLSE). We concentrate on the properties obtained from the direct scattering problem associated with the NLSE. We discuss some general issues regarding soliton creation from random input. We also study the averaged spectral density of random quasilinear waves generated in the NLSE channel for two models of the disordered input field profile. The first model is symmetric complex Gaussian white noise and the second one is a real dichotomous (telegraph) process. For the former model, the closed-form expression for the averaged spectral density is obtained, while for the dichotomous real input we present the small noise perturbative expansion for the same quantity. In the case of the dichotomous input, we also obtain the distribution of minimal pulse width required for a soliton generation. The obtained results can be applied to a multitude of problems including random nonlinear Fraunhoffer diffraction, transmission properties of randomly apodized long period Fiber Bragg gratings, and the propagation of incoherent pulses in optical fibers.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.

Matrix-assisted diffusion-ordered spectroscopy:application of surfactant solutions to the resolution of isomer spectra

The component spectra of a mixture of isomers with nearly identical diffusion coefficients cannot normally be distinguished in a standard diffusion-ordered spectroscopy (DOSY) experiment but can often be easily resolved using matrix-assisted DOSY, in which diffusion behaviour is manipulated by the addition of a co-solute such as a surfactant. Relatively little is currently known about the conditions required for such a separation, for example, how the choice between normal and reverse micelles affects separation or how the isomer structures themselves affect the resolution. The aim of this study was to explore the application of sodium dodecyl sulfate (SDS) normal micelles in aqueous solution and sodium 1,4-bis(2-ethylhexyl)sulfosuccinate (AOT) aggregates in chloroform, at a range of concentrations, to the diffusion resolution of some simple model sets of isomers such as monomethoxyphenols and short chain alcohols. It is shown that SDS micelles offer better resolution where these isomers differ in the position of a hydroxyl group, whereas AOT aggregates are more effective for isomers differing in the position of a methyl group. For both the normal SDS micelles and the less well-defined AOT aggregates, differences in the resolution of the isomers can in part be rationalised in terms of differing degrees of hydrophobicity, amphiphilicity and steric effects. Copyright © 2012 John Wiley & Sons, Ltd.

Matrix-assisted diffusion-ordered spectroscopy:mixture resolution by NMR using SDSmicelles

Diffusion-ordered spectroscopy (DOSY) is a powerful technique for mixture analysis, but in its basic form it cannot separate the component spectra for species with very similar diffusion coefficients. It has been recently demonstrated that the component spectra of a mixture of isomers with nearly identical diffusion coefficients (the three dihydroxybenzenes) can be resolved using matrix-assisted DOSY (MAD), in which diffusion is perturbed by the addition of a co-solute such as a surfactant [R. Evans, S. Haiber, M. Nilsson, G. A. Morris, Anal. Chem. 2009, 81, 4548-4550]. However, little is known about the conditions required for such a separation, for example, the concentrations and concentration ratios of surfactant and solutes. The aim of this study was to explore the concentration range over whichmatrix-assisted DOSY using the surfactant SDS can achieve diffusion resolution of a simple model set of isomers, the monomethoxyphenols. The results show that the separation is remarkably robust with respect to both the concentrations and the concentration ratios of surfactant and solutes, supporting the idea that MAD may become a valuable tool formixture analysis. © 2010 John Wiley & Sons, Ltd.

On the numerical solution of a Cauchy problem for the Laplace equation via a direct integral equation approach

We investigate the problem of determining the stationary temperature field on an inclusion from given Cauchy data on an accessible exterior boundary. On this accessible part the temperature (or the heat flux) is known, and, additionally, on a portion of this exterior boundary the heat flux (or temperature) is also given. We propose a direct boundary integral approach in combination with Tikhonov regularization for the stable determination of the temperature and flux on the inclusion. To determine these quantities on the inclusion, boundary integral equations are derived using Green’s functions, and properties of these equations are shown in an L2-setting. An effective way of discretizing these boundary integral equations based on the Nystr¨om method and trigonometric approximations, is outlined. Numerical examples are included, both with exact and noisy data, showing that accurate approximations can be obtained with small computational effort, and the accuracy is increasing with the length of the portion of the boundary where the additionally data is given.