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.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Distortional buckling of I-beams by finite element method

Distortional buckling, unlike the usual lateral-torsional buckling in which the cross-section remains rigid in its own plane, involves distortion of web in the cross-section. This type of buckling typically occurs in beams with slender web and stocky flanges. Most of the published studies assume the web to deform with a cubic shape function. As this assumption may limit the accuracy of the results, a fifth order polynomial is chosen here for the web displacements. The general line-type finite element model used here has two nodes and a maximum of twelve degrees of freedom per node. The model not only can predict the correct coupled mode but also is capable of handling the local buckling of the web.

Detecting finite bandwidth periodic signals in stationary noise using the signal coherence spectrum

All signals that appear to be periodic have some sort of variability from period to period regardless of how stable they appear to be in a data plot. A true sinusoidal time series is a deterministic function of time that never changes and thus has zero bandwidth around the sinusoid's frequency. A zero bandwidth is impossible in nature since all signals have some intrinsic variability over time. Deterministic sinusoids are used to model cycles as a mathematical convenience. Hinich [IEEE J. Oceanic Eng. 25 (2) (2000) 256-261] introduced a parametric statistical model, called the randomly modulated periodicity (RMP) that allows one to capture the intrinsic variability of a cycle. As with a deterministic periodic signal the RMP can have a number of harmonics. The likelihood ratio test for this model when the amplitudes and phases are known is given in [M.J. Hinich, Signal Processing 83 (2003) 1349-13521. A method for detecting a RMP whose amplitudes and phases are unknown random process plus a stationary noise process is addressed in this paper. The only assumption on the additive noise is that it has finite dependence and finite moments. Using simulations based on a simple RMP model we show a case where the new method can detect the signal when the signal is not detectable in a standard waterfall spectrograrn display. (c) 2005 Elsevier B.V. All rights reserved.

Central Cardiovascular Anatomy and Function in Crocodilia

Abstract: Among the vertebrates, crocodilians have the most complex anatomy of the heart and outflow channels. Their cardiovascular anatomy may also be the most func­tionally sophisticated, combining as it does the best features of both reptilian and mammalian (and avian) systems. The puzzlingly complex "plumbing" of crocodilians has fascinated ana­tomists and physiologists for a very long time, the first paper being that by Panizza (1833). Gradually, with the application of successive techniques of investigation as they became available, its functional significance has become reasonably clear, and the complexity is now revealed as a cardiovascular system of considerable elegance. In this paper I will review the main anatomical features of the heart and outflow channels, discuss what is known about the way they work, and speculate about the probable functional significance.

Simulations of Bose fields at finite temperature

We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.

Pair correlations in a finite-temperature 1D Bose gas

We calculate the two-particle local correlation for an interacting 1D Bose gas at finite temperature and classify various physical regimes. We present the exact numerical solution by using the Yang-Yang equations and Hellmann-Feynman theorem and develop analytical approaches. Our results draw prospects for identifying the regimes of coherent output of an atom laser, and of finite-temperature “fermionization” through the measurement of the rates of two-body inelastic processes, such as photoassociation.

Wolbachia infections are distributed throughout insect somatic and germ line tissues

Wolbachia are intracellular microorganisms that form maternally-inherited infections within numerous arthropod species. These bacteria have drawn much attention, due in part to the reproductive alterations that they induce in their hosts including cytoplasmic incompatibility (CI), feminization and parthenogenesis. Although Wolbachia's presence within insect reproductive tissues has been well described, relatively few studies have examined the extent to which Wolbachia infects other tissues. We have examined Wolbachia tissue tropism in a number of representative insect hosts by western blot, dot blot hybridization and diagnostic PCR. Results from these studies indicate that Wolbachia are much more widely distributed in host tissues than previously appreciated. Furthermore, the distribution of Wolbachia in somatic tissues varied between different Wolbachia/host associations. Some associations showed Wolbachia disseminated throughout most tissues while others appeared to be much more restricted, being predominantly limited to the reproductive tissues. We discuss the relevance of these infection patterns to the evolution of Wolbachia/host symbioses and to potential applied uses of Wolbachia.

Measurement Function Design for Visual Tracking Applications

Extracting human postural information from video sequences has proved a difficult research question. The most successful approaches to date have been based on particle filtering, whereby the underlying probability distribution is approximated by a set of particles. The shape of the underlying observational probability distribution plays a significant role in determining the success, both accuracy and efficiency, of any visual tracker. In this paper we compare approaches used by other authors and present a cost path approach which is commonly used in image segmentation problems, however is currently not widely used in tracking applications.

A unified and complete construction of all finite dimensional irreducible representations of gl(2|2)

Representations of the non-semisimple superalgebra gl(2/2) in the standard basis are investigated by means of the vector coherent state method and boson-fermion realization. All finite-dimensional irreducible typical and atypical representations and lowest weight (indecomposable) Kac modules of gl(2/2) are constructed explicity through the explicit construction of all gl(2) circle plus gl(2) particle states (multiplets) in terms of boson and fermion creation operators in the super-Fock space. This gives a unified and complete treatment of finite-dimensional representations of gl(2/2) in explicit form, essential for the construction of primary fields of the corresponding current superalgebra at arbitrary level.

A demonstration of artificial dissipation effects in some finite volume solution methods for the Navier-Stokes equations

The artificial dissipation effects in some solutions obtained with a Navier-Stokes flow solver are demonstrated. The solvers were used to calculate the flow of an artificially dissipative fluid, which is a fluid having dissipative properties which arise entirely from the solution method itself. This was done by setting the viscosity and heat conduction coefficients in the Navier-Stokes solvers to zero everywhere inside the flow, while at the same time applying the usual no-slip and thermal conducting boundary conditions at solid boundaries. An artificially dissipative flow solution is found where the dissipation depends entirely on the solver itself. If the difference between the solutions obtained with the viscosity and thermal conductivity set to zero and their correct values is small, it is clear that the artificial dissipation is dominating and the solutions are unreliable.