A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

The transform likelihood ratio method for rare event simulation with heavy tails

We present a novel method, called the transform likelihood ratio (TLR) method, for estimation of rare event probabilities with heavy-tailed distributions. Via a simple transformation ( change of variables) technique the TLR method reduces the original rare event probability estimation with heavy tail distributions to an equivalent one with light tail distributions. Once this transformation has been established we estimate the rare event probability via importance sampling, using the classical exponential change of measure or the standard likelihood ratio change of measure. In the latter case the importance sampling distribution is chosen from the same parametric family as the transformed distribution. We estimate the optimal parameter vector of the importance sampling distribution using the cross-entropy method. We prove the polynomial complexity of the TLR method for certain heavy-tailed models and demonstrate numerically its high efficiency for various heavy-tailed models previously thought to be intractable. We also show that the TLR method can be viewed as a universal tool in the sense that not only it provides a unified view for heavy-tailed simulation but also can be efficiently used in simulation with light-tailed distributions. We present extensive simulation results which support the efficiency of the TLR method.

Shift, scaling and derivative properties for the discrete cosine transform

A set of DCT domain properties for shifting and scaling by real amounts, and taking linear operations such as differentiation is described. The DCT coefficients of a sampled signal are subjected to a linear transform, which returns the DCT coefficients of the shifted, scaled and/or differentiated signal. The properties are derived by considering the inverse discrete transform as a cosine series expansion of the original continuous signal, assuming sampling in accordance with the Nyquist criterion. This approach can be applied in the signal domain, to give, for example, DCT based interpolation or derivatives. The same approach can be taken in decoding from the DCT to give, for example, derivatives in the signal domain. The techniques may prove useful in compressed domain processing applications, and are interesting because they allow operations from the continuous domain such as differentiation to be implemented in the discrete domain. An image matching algorithm illustrates the use of the properties, with improvements in computation time and matching quality.

Quantum noise in optical fibers. I. Stochastic equations

We analyze the quantum dynamics of radiation propagating in a single-mode optical fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We start from a fundamental Hamiltonian that includes the principal known nonlinear effects and quantum-noise sources, including linear gain and loss. Both Markovian and frequency-dependent, non-Markovian reservoirs are treated. This treatment allows quantum Langevin equations, which have a classical form except for additional quantum-noise terms, to be calculated. In practical calculations, it is more useful to transform to Wigner or 1P quasi-probability operator representations. These transformations result in stochastic equations that can be analyzed by use of perturbation theory or exact numerical techniques. The results have applications to fiber-optics communications, networking, and sensor technology.

Entanglement in a simple quantum phase transition

What entanglement is present in naturally occurring physical systems at thermal equilibrium? Most such systems are intractable and it is desirable to study simple but realistic systems that can be solved. An example of such a system is the one-dimensional infinite-lattice anisotropic XY model. This model is exactly solvable using the Jordan-Wigner transform, and it is possible to calculate the two-site reduced density matrix for all pairs of sites. Using the two-site density matrix, the entanglement of formation between any two sites is calculated for all parameter values and temperatures. We also study the entanglement in the transverse Ising model, a special case of the XY model, which exhibits a quantum phase transition. It is found that the next-nearest-neighbor entanglement (though not the nearest-neighbor entanglement) is a maximum at the critical point. Furthermore, we show that the critical point in the transverse Ising model corresponds to a transition in the behavior of the entanglement between a single site and the remainder of the lattice.

Multicomponent reversible reactions inside a porous catalyst pellet

This work deals with a solution method to handle multicomponents reversible reactions occurring inside a porous catalyst pellet. The complexity of this problem arises from the fact that the effective diffusivities and Biot number, which characterizes the external mass transfer, are different for each chemical species. In mathematical terms, this means that each chemical species has its own subspace and, therefore, when the technique of finite integral transform is applied to solve this multicomponent problem, each chemical species is associated with its own integral transform kernel. The analytical solutions obtained for this problem are compact and simple for any further manipulation. Application of this result to the catalytic reforming of C7 hydrocarbon system is shown in this paper.

Approximate analytical solutions for porous catalysts with nonuniform activity

Using a novel finite integral transform technique, the problem of diffusion and chemical reaction in a porous catalyst with general activity profile is investigated theoretically. Analytical expressions for the effectiveness factor are obtained for pth order and Michaelis-Menten kinetics. Perturbation methods are employed to provide useful asymptotic solutions for large or small values of Thiele modulus and Biot number.

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).

An improved novel method for solving nonlinear boundary value problems

A modified formula for the integral transform of a nonlinear function is proposed for a class of nonlinear boundary value problems. The technique presented in this paper results in analytical solutions. Iterations and initial guess, which are needed in other techniques, are not required in this novel technique. The analytical solutions are found to agree surprisingly well with the numerically exact solutions for two examples of power law reaction and Langmuir-Hinshelwood reaction in a catalyst pellet.

The myth of the nation: Historiography of New Zealand architecture since 1907

This paper considers the relationship between the recent historiography (of the last quarter century) of “New Zealand architecture” and the historical notion of “New Zealand-ness” invoked in contemporary architecture. It argues that a more recent programmatic uptake of post-War discussions on national identity and regional specificity has fed the tendencies of practicing architects to defer to history in rhetorical defences of their work: the beach-side mansion as a contemporary expression of the 1950s bach; a formal modernism divorced from the social discourse adherent to the historical moment that it “restates”; and so on. The paper will consider instances in the historiography of New Zealand architecture where historians have compounded, consciously or accidentally, a problem that is systemic to the uses made by architects of historical knowledge (in the most general examples), identifying the difficulties of relying upon the tentative conclusions of an under-studied field in developing principles of contemporary architectural practice under the banners of New Zealand-ness, regionalism, or localism, or with reference to icons of New Zealand architectural history. At the heart of this paper is a reflection on historiographical responsibility in presenting knowledge of a national past to an audience that is eager to transform that knowledge into principles of contemporary production. What, the paper asks, is the historical basis for speaking of a New Zealand architecture? Can we speak of a national history of architecture distinct from a regional history, or from an international history of architecture?

Bifurcation in growth patterns for arrays of parallel Griffith, edge and sliding cracks

This paper presents the recent finding by Muhlhaus et al [1] that bifurcation of crack growth patterns exists for arrays of two-dimensional cracks. This bifurcation is a result of the nonlinear effect due to crack interaction, which is, in the present analysis, approximated by the dipole asymptotic or pseudo-traction method. The nonlinear parameter for the problem is the crack length/ spacing ratio lambda = a/h. For parallel and edge crack arrays under far field tension, uniform crack growth patterns (all cracks having same size) yield to nonuniform crack growth patterns (i.e. bifurcation) if lambda is larger than a critical value lambda(cr) (note that such bifurcation is not found for collinear crack arrays). For parallel and edge crack arrays respectively, the value of lambda(cr) decreases monotonically from (2/9)(1/2) and (2/15.096)(1/2) for arrays of 2 cracks, to (2/3)(1/2)/pi and (2/5.032)(1/2)/pi for infinite arrays of cracks. The critical parameter lambda(cr) is calculated numerically for arrays of up to 100 cracks, whilst discrete Fourier transform is used to obtain the exact solution of lambda(cr) for infinite crack arrays. For geomaterials, bifurcation can also occurs when array of sliding cracks are under compression.

Carrageenans with complex substitution patterns from red algae of the genus Erythroclonium

Extracellular polysaccharides from three Erythroclonium spp. were shown, by a combination of compositional, linkage analyses, and Fourier transform infrared and C-13-nuclear magnetic resonance spectroscopy, to be highly substituted carrageenans with at least five types of repeating disaccharide units. These are the carrabiose 2,4'-disulfate of iota-carrageenan, carrabiose 2-sulfate of alpha-carrageenan, the 6'-O-methylated counterparts of each of these repeating units, and 4',6'-O-(1-carboxyethylidene)carrabiose 2-sulfate. The polysaccharides also contain significant amounts of unsubstituted, 4-linked galactopyranose and small amounts of 4-linked 3-O-methylgalactopyranose and terminal glycosyl residues. The carrageenan preparations of the three species are similar, differing only in the proportions of some components. (C) 1998 Elsevier Science Ltd.