Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

Generalized optimal lattice covering of finite-dimensional Euclidean space

A generalization of the classical problem of optimal lattice covering of R-n is considered. Solutions to this generalized problem are found in two specific classes of lattices. The global optimal solution of the generalization is found for R-2. (C) 1998 Elsevier Science Inc. All rights reserved.

5-cycle systems of K-v-F with a hole

Recently the problem of the existence of a 5-cycle system of K-v with a hole of size u was completely solved. In this paper we prove necessary and sufficient conditions on v and u for the existence of a 5-cycle system of K-v - F, with a hole of size u.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

The conserved quantity theory of causation and chance raising

In this paper I offer an 'integrating account' of singular causation, where the term 'integrating' refers to the following program for analysing causation. There are two intuitions about causation, both of which face serious counterexamples when used as the basis for an analysis of causation. The 'process' intuition, which says that causes and effects are linked by concrete processes, runs into trouble with cases of misconnections', where an event which serves to prevent another fails to do so on a particular occasion and yet the two events are linked by causal processes. The chance raising intuition, according to which causes raise the chance of their effects, easily accounts for misconnections but faces the problem of chance lowering causes, a problem easily accounted for by the process approach. The integrating program attempts to provide an analysis of singular causation by synthesising the two insights, so as to solve both problems. In this paper I show that extant versions of the integrating program due to Eells, Lewis, and Menzies fail to account for the chance-lowering counterexample. I offer a new diagnosis of the chance lowering case, and use that as a basis for an integrating account of causation which does solve both cases. In doing so, I accept various assumptions of the integrating program, in particular that there are no other problems with these two approaches. As an example of the process account, I focus on the recent CQ theory of Wesley Salmon (1997).

Bounds on integrals of the Wigner function

The integral of the Wigner function over a subregion of the phase space of a quantum system may be less than zero or greater than one. It is shown that for systems with 1 degree of freedom, the problem of determining the best possible upper and lower bounds on such an integral, over an possible states, reduces to the problem of finding the greatest and least eigenvalues of a Hermitian operator corresponding to the subregion. The problem is solved exactly in the case of an arbitrary elliptical region. These bounds provide checks on experimentally measured quasiprobability distributions.

Public health implications of new guidelines for lead in drinking water: a case study in an area with historically high water lead levels

Concern about the neurotoxicity of lead, particularly in infants and young children, has led to a revision of blood lead levels which are considered to involve an acceptable level of human exposure. Drinking water guidelines have also been reviewed in order to reduce this source of population exposure to lead. In the last 20 years, guidelines have been reduced from 100 to 50 to 10 mu g/litre. Lead in tap water used to be a major public health problem in Glasgow because of the high prevalence of houses with lead service pipes, the low pH of the public water supply and the resulting high levels of lead in water used for public consumption. Following two separate programmes of water treatment, involving the addition of lime and, a decade later, lime supplemented with orthophosphate, it is considered that maximal measures have been taken to reduce lead exposure by chemical treatment of the water supply. Any residual problem of public exposure would require large scale replacement of lead service pipes. In anticipation of the more stringent limits for lead in drinking water, we set out to measure current lead exposure From tap water in the population of Glasgow served by the Loch Katrine water supply. to compare the current situation with 12 years previously and to assess the public health implications of different limits. The study was based on mothers of young children since maternal blood lead concentrations and the domestic water that mothers use to prepare bottle feeds are the principal sources of foetal and infant lead exposure. An estimated 17% of mothers lived in households with tap water lead concentrations of 10 mu g/litre (the WHO guideline) or above in 1993 compared with 49% in 1981. Mean maternal blood lead concentrations fell by 69% in 12 years. For a given water lead concentration, maternal blood lead concentrations were 67% lower. The mean maternal blood lead concentration was 3.7 mu g/litre in the population at large, compared with 3.3 mu g/litre in households with negligible or absent tap water lead. Nevertheless, between 63% and 76% of cases of mothers with blood lead concentrations of 10 mu g/dl or above were attributable to tap water lead. The study found that maternal blood lead concentrations were well within limits currently considered safe for human health. About 15% of infants may be exposed via bottle feeds to tap water lead concentrations that exceed the WHO guideline of 10 mu g/litre. In the context of the health and social problems which affect the well-being and development of infants and children in Glasgow, however, current levels of lend exposure are considered to present a relatively minor health problem. (C) 2000 Elsevier Science Ltd. All rights reserved.

The number of 4-cycles in 2-factorizations of K2n minus a 1-factor

In this paper we give a complete solution to problem of determining the number of 4-cycles in a 2-factorization of K-2n\ 1-factor. (C) 2000 Elsevier Science B.V. All rights reserved.

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.

Interpretation of negative values of the interaction parameter in the adsorption equation through the effects of surface layer heterogeneity

The problem of the negative values of the interaction parameter in the equation of Frumkin has been analyzed with respect to the adsorption of nonionic molecules on energetically homogeneous surface. For this purpose, the adsorption states of a homologue series of ethoxylated nonionic surfactants on air/water interface have been determined using four different models and literature data (surface tension isotherms). The results obtained with the Frumkin adsorption isotherm imply repulsion between the adsorbed species (corresponding to negative values of the interaction parameter), while the classical lattice theory for energetically homogeneous surface (e.g., water/air) admits attraction alone. It appears that this serious contradiction can be overcome by assuming heterogeneity in the adsorption layer, that is, effects of partial condensation (formation of aggregates) on the surface. Such a phenomenon is suggested in the Fainerman-Lucassen-Reynders-Miller (FLM) 'Aggregation model'. Despite the limitations of the latter model (e.g., monodispersity of the aggregates), we have been able to estimate the sign and the order of magnitude of Frumkin's interaction parameter and the range of the aggregation numbers of the surface species. (C) 2004 Elsevier B.V All rights reserved.

d-cube decompositions of K-n/K-m

Necessary conditions on n, m and d are given for the existence of an edge-disjoint decomposition of K-n\K-m into copies of the graph of a d-dimensional cube. Sufficiency is shown when d = 3 and, in some cases, when d = 2(t). We settle the problem of embedding 3-cube decompositions of K-m into 3-cube decompositions of K-n; where n greater than or equal to m.

Standard errors of fitted component means of normal mixtures

In this paper use consider the problem of providing standard errors of the component means in normal mixture models fitted to univariate or multivariate data by maximum likelihood via the EM algorithm. Two methods of estimation of the standard errors are considered: the standard information-based method and the computationally-intensive bootstrap method. They are compared empirically by their application to three real data sets and by a small-scale Monte Carlo experiment.

Embeddings of m-cycle systems and incomplete m-cycle systems: m<=14

In this paper we completely settle the embedding problem for m-cycle systems with m less than or equal to 14. We also solve the more general problem of finding m-cycle systems of K-v - K-u when m is an element of {4,6,7,8,10,12,14}.

Interconnected-tubes model of hepatic elimination

The distributed-tubes model of hepatic elimination is extended to include intermixing between sinusoids, resulting in the formulation of a new, interconnected-tubes model. The new model is analysed for the simple case of two interconnected tubes, where an exact solution is obtained. For the case of many strongly-interconnected tubes, it is shown that a zeroth-order approximation leads to the convection-dispersion model. As a consequence the dispersion number is expressed, for the first time, in terms of its main physiological determinants: heterogeneity of flow and density of interconnections between sinusoids. The analysis of multiple indicator dilution data from a perfused liver preparation using the simplest version of the model yields the estimate 10.3 for the average number of interconnections. The problem of boundary conditions for the dispersion model is considered from the viewpoint that the dispersion-convection equation is a zeroth-order approximation to the equations for the interconnected-tubes model. (C) 1997 Academic Press Limited.

An Economic Analysis of the Legal Concept of Trustees' Power of Investment

Authorised trustee status is a legal concept which has economic implications; one of the major implications is that it assists in the direction of investment funds into particular securities and areas of the economy. The concept of authorised trustee status, while attempting to achieve specific outcomes for the beneficiaries of trusts cannot be relied upon to secure these results. Economic analysis of the role of the trustee maintains that this role is one of portfolio manager; a role which is complex but which is explicable in terms of definable procedures and practices. The role of trustee as portfolio manager is one which requires greater financial knowledge than can be assumed is possessed by all trustees. The trustee as portfolio manager is required to maintain a review of decisions make under powers to invest trust assets. A solution to the problem of authorised trustee status is proposed. The solution takes two parts: the first is the adoption of the prudent person approach but with the codification of duties of the trustee and the explicit listing of the factors that a trustee should consider in using the investment powers. The second part of the proposed solution is to link the investment powers of trustees to the best practice features of securities advisers who are now licensed by a regulatory body, the Australian Securities and Investment Commission.