General Harris regularity criterion for non-linear Markov branching processes

We extend the Harris regularity condition for ordinary Markov branching process to a more general case of non-linear Markov branching process. A regularity criterion which is very easy to check is obtained. In particular, we prove that a super-linear Markov branching process is regular if and only if the per capita offspring mean is less than or equal to I while a sub-linear Markov branching process is regular if the per capita offspring mean is finite. The Harris regularity condition then becomes a special case of our criterion.

An application of quasi-Newton methods for the numerical solution of interface problems

Quasi-Newton methods are applied to solve interface problems which arise from domain decomposition methods. These interface problems are usually sparse systems of linear or nonlinear equations. We are interested in applying these methods to systems of linear equations where we are not able or willing to calculate the Jacobian matrices as well as to systems of nonlinear equations resulting from nonlinear elliptic problems in the context of domain decomposition. Suitability for parallel implementation of these algorithms on coarse-grained parallel computers is discussed.

Fitting host-parasitoid models with CV² > 1 using hierarchical generalized linear models

The powerful general Pacala-Hassell host-parasitoid model for a patchy environment, which allows host density–dependent heterogeneity (HDD) to be distinguished from between-patch, host density–independent heterogeneity (HDI), is reformulated within the class of the generalized linear model (GLM) family. This improves accessibility through the provision of general software within well–known statistical systems, and allows a rich variety of models to be formulated. Covariates such as age class, host density and abiotic factors may be included easily. For the case where there is no HDI, the formulation is a simple GLM. When there is HDI in addition to HDD, the formulation is a hierarchical generalized linear model. Two forms of HDI model are considered, both with between-patch variability: one has binomial variation within patches and one has extra-binomial, overdispersed variation within patches. Examples are given demonstrating parameter estimation with standard errors, and hypothesis testing. For one example given, the extra-binomial component of the HDI heterogeneity in parasitism is itself shown to be strongly density dependent.

Numerical experiments of some Krylov subspace methods for black oil model

This paper discusses preconditioned Krylov subspace methods for solving large scale linear systems that originate from oil reservoir numerical simulations. Two types of preconditioners, one being based on an incomplete LU decomposition and the other being based on iterative algorithms, are used together in a combination strategy in order to achieve an adaptive and efficient preconditioner. Numerical tests show that different Krylov subspace methods combining with appropriate preconditioners are able to achieve optimal performance.

A finite volume approach to geometrically non-linear stress analysis

In this past decade finite volume (FV) methods have increasingly been used for the solution of solid mechanics problems. This contribution describes a cell vertex finite volume discretisation approach to the solution of geometrically nonlinear (GNL) problems. These problems, which may well have linear material properties, are subject to large deformation. This requires a distinct formulation, which is described in this paper together with the solution strategy for GNL problem. The competitive performance for this procedure against the conventional finite element (FE) formulation is illustrated for a three dimensional axially loaded column.

A vertex-based finite volume method applied to non-linear material problems in computational solid mechanics

A vertex-based finite volume (FV) method is presented for the computational solution of quasi-static solid mechanics problems involving material non-linearity and infinitesimal strains. The problems are analysed numerically with fully unstructured meshes that consist of a variety of two- and threedimensional element types. A detailed comparison between the vertex-based FV and the standard Galerkin FE methods is provided with regard to discretization, solution accuracy and computational efficiency. For some problem classes a direct equivalence of the two methods is demonstrated, both theoretically and numerically. However, for other problems some interesting advantages and disadvantages of the FV formulation over the Galerkin FE method are highlighted.

Domain decomposition methods for welding problems

This paper, a 2-D non-linear electric arc-welding problem is considered. It is assumed that the moving arc generates an unknown quantity of energy which makes the problem an inverse problem with an unknown source. Robust algorithms to solve such problems e#ciently, and in certain circumstances in real-time, are of great technological and industrial interest. There are other types of inverse problems which involve inverse determination of heat conductivity or material properties [CDJ63][TE98], inverse problems in material cutting [ILPP98], and retrieval of parameters containing discontinuities [IK90]. As in the metal cutting problem, the temperature of a very hot surface is required and it relies on the use of thermocouples. Here, the solution scheme requires temperature measurements lied in the neighbourhood of the weld line in order to retrieve the unknown heat source. The size of this neighbourhood is not considered in this paper, but rather a domain decomposition concept is presented and an examination of the accuracy of the retrieved source are presented. This paper is organised as follows. The inverse problem is formulated and a method for the source retrieval is presented in the second section. The source retrieval method is based on an extension of the 1-D source retrieval method as proposed in [ILP].

Uniqueness criteria for continuous-time Markov chains with general transition structures

We derive necessary and sufficient conditions for the existence of bounded or summable solutions to systems of linear equations associated with Markov chains. This substantially extends a famous result of G. E. H. Reuter, which provides a convenient means of checking various uniqueness criteria for birth-death processes. Our result allows chains with much more general transition structures to be accommodated. One application is to give a new proof of an important result of M. F. Chen concerning upwardly skip-free processes. We then use our generalization of Reuter's lemma to prove new results for downwardly skip-free chains, such as the Markov branching process and several of its many generalizations. This permits us to establish uniqueness criteria for several models, including the general birth, death, and catastrophe process, extended branching processes, and asymptotic birth-death processes, the latter being neither upwardly skip-free nor downwardly skip-free.

Survival in 176 patients with inoperable non small cell lung cancer treated initially with chemotherapy in a district general hospital

In advanced non-small cell lung cancer (NSCLC) platinum based chemotherapy with second generation drugs improves median survival (MS) to 8 months and 29% and 10% at 1 and 2 years. Platinum with a third generation drug can improve survival further (BMJ 1995;311: 899) (Spiro et al. Thorax 2004;59:828 Big Lung Trial; N Engl J Med 2003;346:92 ECOG study). NICE now recommends chemotherapy with platinum and a third generation drug for inoperable NSCLC as the first treatment modality. Methods: We audited survival of 176/461 consecutive patients referred for at least 3 courses of platinum and either gemcitabine or vinorelbine from July 2001 to December 2005. Minimal follow up 17 months. Chemotherapy was given on site. Radical radiotherapy for stage IIIA, palliative radiotherapy and second line drugs were given as felt appropriate. Results: 64% were male. 30 (17%) were <55 years ; 66 (37.5%) age 55–65 years; 63 (35.8%) aged 66–75 and 16 (9.1%) >75 years. 5 (2.8%) were stage II; 46 (26%) stage IIIA; 68 (38%) stage IIIB and 55 (30.8%) stage IV. 68 (38%) had 0– 2 courses; 63 (36%) 3 courses and 44 (25%) had 4 or more.

Application of regularization methods to damage detection in large scale plane frame structures using incomplete noisy modal data

This paper presents an approach for detecting local damage in large scale frame structures by utilizing regularization methods for ill-posed problems. A direct relationship between the change in stiffness caused by local damage and the measured modal data for the damaged structure is developed, based on the perturbation method for structural dynamic systems. Thus, the measured incomplete modal data can be directly adopted in damage identification without requiring model reduction techniques, and common regularization methods could be effectively employed to solve the developed equations. Damage indicators are appropriately chosen to reflect both the location and severity of local damage in individual components of frame structures such as in brace members and at beam-column joints. The Truncated Singular Value Decomposition solution incorporating the Generalized Cross Validation method is introduced to evaluate the damage indicators for the cases when realistic errors exist in modal data measurements. Results for a 16-story building model structure show that structural damage can be correctly identified at detailed level using only limited information on the measured noisy modal data for the damaged structure.

An investigation into the rheological properties of different lead-free solder pastes for surface mount applications

Purpose – The purpose of this paper is to investigate the rheological behaviour of three different lead-free solder pastes used for surface mount applications in the electronic industry.Design/methodology/approach – This study concerns the rheological measurements of solder paste samples and is made up of three parts. The first part deals with the measurement of rhelogical properties with three different measuring geometries, the second part looks into the effect of frequencies on oscillatory stress sweep measurements and the final part reports on the characterisation and comparison of three different types of Pb-free solder pastes. Findings – Among the three geometries, the serrated parallel plate was found effective in minimising the wall-slip effect. From the oscillatory stresssweep data with different frequencies; it was observed that the linear visco-elastic region is independent of frequency for all the solder paste samples. To understand the shear thinning behaviour of solder paste, the well known Cross and Carreau models were fitted to the viscosity data. Moreover,creep-recovery and dynamic frequency-sweep tests were also carried out without destroying the sample’s structure and have yielded useful information on the pastes behaviour.Research limitations/implications – More extensive research is needed to fully characterise the wall-slip behaviour during the rheological measurements of solder pastes. Practical implications – The rheological test results presented in this paper will be of important value for research and development, quality control and facilitation of the manufacturing of solder pastes and flux mediums. Originality/value – This paper shows how wall-slip effects can be effectively avoided during rheological measurements of solder pastes. The paper also outlines how different rheological test methods can be used to characterise solder paste behaviours