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.

Modeling diffraction in free-space optical interconnects by the mode expansion method

Free-space optical interconnects (FSOIs), made up of dense arrays of vertical-cavity surface-emitting lasers, photodetectors and microlenses can be used for implementing high-speed and high-density communication links, and hence replace the inferior electrical interconnects. A major concern in the design of FSOIs is minimization of the optical channel cross talk arising from laser beam diffraction. In this article we introduce modifications to the mode expansion method of Tanaka et al. [IEEE Trans. Microwave Theory Tech. MTT-20, 749 (1972)] to make it an efficient tool for modelling and design of FSOIs in the presence of diffraction. We demonstrate that our modified mode expansion method has accuracy similar to the exact solution of the Huygens-Kirchhoff diffraction integral in cases of both weak and strong beam clipping, and that it is much more accurate than the existing approximations. The strength of the method is twofold: first, it is applicable in the region of pronounced diffraction (strong beam clipping) where all other approximations fail and, second, unlike the exact-solution method, it can be efficiently used for modelling diffraction on multiple apertures. These features make the mode expansion method useful for design and optimization of free-space architectures containing multiple optical elements inclusive of optical interconnects and optical clock distribution systems. (C) 2003 Optical Society of America.

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.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Seedling responses of four Australian tree species to toxic concentrations of manganese in solution culture

Little is known about the responses of Australian plants to excess metal, including Mn. It is important to remedy this lack of information so that knowledgeable decisions can be made about managing Mn contaminated sites where inhabited by Australian vegetation. Acacia holosericea, Melaleuca leucadendra, Eucalyptus crebra and Eucalyptus camaldulensis were grown in dilute solution culture for 10 weeks. The seedlings ( 42 days old) were exposed to six Mn treatments viz., 1, 8, 32, 128, 512 and 2048 muM. The order of tolerance to toxic concentrations of Mn was A. holosericea congruent to = E. crebra < M. leucadendra < E. camaldulensis, the critical external concentrations being approximately 5.1, 5.0, 21 and 330 muM, respectively. The critical tissue Mn concentrations for the youngest fully expanded leaf and total shoots were, respectively, 265 and 215 mug g(-1) DM for A. holosericea, 445 and 495 mug g(-1) DM for M. leucadendra, 495 and 710 mug g(-1) DM for E. crebra and 7230 and 6510 mug g(-1) DM for E. camaldulensis. The high tolerance of E. camaldulensis ( as opposed to the sensitivity of E. crebra) to excess Mn raises concern about fauna feeding on the plant and is consistent with hypotheses suggesting the Eucalyptus subgenus Symphomyrtus is particularly tolerant of stress, including excess Mn. The results from this paper provide the first comprehensive combination of growth responses, critical external concentrations, critical tissue concentrations and plant toxicity symptoms for three important Australian genera, viz., Eucalyptus, Acacia and Melaleuca, for use in the management of Mn toxic sites.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.

Seedling responses of three Australian tree species to toxic concentrations of zinc in solution culture

A frequently desired outcome when rehabilitating Zn toxic sites in Australia is to establish a self-sustaining native ecosystem. Hence, it is important to understand the tolerance of Australian native plants to high concentrations of Zn. Very little is known about the responses of Australian native plants, and trees in particular, to toxic concentrations of Zn. Acacia holosericea, Eucalyptus camaldulensis and Melaleuca leucadendra plants were grown in dilute solution culture for 10 weeks. The seedlings (42 days old) were exposed to six Zn treatments viz., 0.5, 5, 10, 25, 50 and 100 muM. The order of tolerance to toxic concentrations of Zn was E. camaldulensis > A. holosericea > M. leucadendra, the critical external concentrations being approximately 20, 12 and 1.5 muM, respectively. Tissue Zn concentrations increased as solution Zn increased for all species. Root tissue concentrations were higher than shoot tissue concentrations at all solution Zn concentrations. The critical tissue Zn concentrations were approximately 85 and 110 mug g(-1) DM for M. leucadendra, 115 and 155 mug g(-1) DM for A. holosericea and 415 and 370 mug g(-1) DM for E. camaldulensis for the youngest fully expanded leaf and total shoots, respectively. The results from this paper provide the first comprehensive combination of growth responses, critical external concentrations, critical tissue concentrations and plant toxicity symptoms for three important Australian genera, viz., Eucalyptus, Acacia and Melaleuca, for use in the rehabilitation of potentially Zn toxic sites.

Responses of four Australian tree species to toxic concentrations of copper in solution culture

In Australia, metal-contaminated sites, including those with elevated levels of copper (Cu), are frequently revegetated with endemic plants. Little is known about the responses of Australian plants to excess Cu. Acacia holosericea, Eucalyptus crebra, Eucalyptus camaldulensis, and Melaleuca leucadendra were grown in solution culture with six Cu treatments (0.1 to 40 mu M). While A. holosericea was the most tolerant to excess Cu, all of the species tested were sensitive to excess Cu when compared with exotic tree and agricultural species. The critical external concentrations for toxicity were < 0.7 mu M for all species tested. There was little differentiation between shoot-tissue Cu concentrations in normal versus treated plants, thus, the derivation of critical shoot concentrations was possible only for the most tolerant species, A. holosericea. Critical root Cu concentrations were approximately 210 mu g g(-1) (A. holosericea), 150 mu g g(-1) (E. crebra), 25 mu g g(-1) (E. camaldulensis), and 165 mu g g(-1) (M. leucadendra). These results provide the first comprehensive combination of growth responses, critical concentrations, and toxicity symptoms for three important Australian genera for use in the management of Cu-contaminated sites.

BM+-tree: A hyperplane-based index method for high-dimensional metric spaces

In this paper, we propose a novel high-dimensional index method, the BM+-tree, to support efficient processing of similarity search queries in high-dimensional spaces. The main idea of the proposed index is to improve data partitioning efficiency in a high-dimensional space by using a rotary binary hyperplane, which further partitions a subspace and can also take advantage of the twin node concept used in the M+-tree. Compared with the key dimension concept in the M+-tree, the binary hyperplane is more effective in data filtering. High space utilization is achieved by dynamically performing data reallocation between twin nodes. In addition, a post processing step is used after index building to ensure effective filtration. Experimental results using two types of real data sets illustrate a significantly improved filtering efficiency.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.