Toward an understanding of the turbidity measurement of fleterocoagulation rate constants of dispersions containing particles of different sizes

Our previous studies have shown that the determination of coagulation rate constants by turbidity measurement becomes impossible for a certain operating wavelength (that is, its blind point) because at this wavelength the change in the turbidity of a dispersion completely loses its response to the coagulation process. Therefore, performing the turbidity measurement in the wavelength range near the blind point should be avoided. In this article, we demonstrate that the turbidity measurement of the rate constant for coagulation of a binary dispersion containing particles of two different sizes (heterocoagulation) presents special difficulties because the blind point shifts with not only particle size but also with the component fraction. Some important aspects of the turbidity measurement for the heterocoagulation rate constant are discussed and experimentally tested. It is emphasized that the T-matrix method can be used to correctly evaluate extinction cross sections of doublets formed during the heterocoagulation process, which is the key data determining the rate constant from the turbidity measurement, and choosing the appropriate operating wavelength and component fraction are important to achieving a more accurate rate constant. Finally, a simple scheme in experimentally determining the sensitivity of the turbidity changes with coagulation over a wavelength range is proposed.

Sequential Monte Carlo computation of the score and observed information matrix in state-space models with application to parameter estimation

Sequential Monte Carlo (SMC) methods are popular computational tools for Bayesian inference in non-linear non-Gaussian state-space models. For this class of models, we propose SMC algorithms to compute the score vector and observed information matrix recursively in time. We propose two different SMC implementations, one with computational complexity $\mathcal{O}(N)$ and the other with complexity $\mathcal{O}(N^{2})$ where $N$ is the number of importance sampling draws. Although cheaper, the performance of the $\mathcal{O}(N)$ method degrades quickly in time as it inherently relies on the SMC approximation of a sequence of probability distributions whose dimension is increasing linearly with time. In particular, even under strong \textit{mixing} assumptions, the variance of the estimates computed with the $\mathcal{O}(N)$ method increases at least quadratically in time. The $\mathcal{O}(N^{2})$ is a non-standard SMC implementation that does not suffer from this rapid degrade. We then show how both methods can be used to perform batch and recursive parameter estimation.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

The annular crack in a nonhomogeneous matrix surrounding a fiber under torsional loading .2. The crack going through the bimaterial interface into the fiber

The axisymmetric problem of an elastic fiber perfectly bonded to a nonhomogeneous elastic matrix which contains an annular crack going through the interface into the fiber under axially symmetric shear stress is considered. The nature of the stress singularity is studied. It is shown that at the irregular point on the interface, whether the shear modulus is continuous or discontinuous the stresses are bounded. The problem is formulated in terms of a singular integral equation and can be solved by a regular method. The stress intensity factors and crack surface displacement are given.

Upwind compact scheme with dispersion control for aerodynamic equations

A compact upwind scheme with dispersion control is developed using a dissipation analogy of the dispersion term. The term is important in reducing the unphysical fluctuations in numerical solutions. The scheme depends on three free parameters that may be used to regulate the size of dissipation as well as the size and direction of dispersion. A coefficient to coordinate the dispersion is given. The scheme has high accuracy, the method is simple, and the amount of computation is small. It also has a good capability of capturing shock waves. Numerical experiments are carried out with two-dimensional shock wave reflections and the results are very satisfactory.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.