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

Unbiased Linear Property Estimation For Spheres From Sections Exhibiting Overprojection And Truncation

In the biological sciences, stereological techniques are frequently used to infer changes in structural parameters (volume fraction, for example) between samples from different populations or subject to differing treatment regimes. Non-homogeneity of these parameters is virtually guaranteed, both between experimental animals and within the organ under consideration. A two-stage strategy is then desirable, the first stage involving unbiased estimation of the required parameter, separately for each experimental unit, the latter being defined as a subset of the organ for which homogeneity can reasonably be assumed. In the second stage, these point estimates are used as data inputs to a hierarchical analysis of variance, to distinguish treatment effects from variability between animals, for example. Techniques are therefore required for unbiased estimation of parameters from potentially small numbers of sample profiles. This paper derives unbiased estimates of linear properties in one special case—the sampling of spherical particles by transmission microscopy, when the section thickness is not negligible and the resulting circular profiles are subject to lower truncation. The derivation uses the general integral equation formulation of Nicholson (1970); the resulting formulae are simplified, algebraically, and their efficient computation discussed. Bias arising from variability in slice thickness is shown to be negligible in typical cases. The strategy is illustrated for data examining the effects, on the secondary lysosomes in the digestive cells, of exposure of the common mussel to hydrocarbons. Prolonged exposure, at 30 μg 1−1 total oil-derived hydrocarbons, is seen to increase the average volume of a lysosome, and the volume fraction that lysosomes occupy, but to reduce their number.

Can diagnostic triage by general practitioners or rheumatology nurses improve the positive predictive value of referrals to early arthritis clinics?

Objectives: To determine whether diagnostic triage by general practitioners (GPs) or rheumatology nurses (RNs) can improve the positive predictive value of referrals to early arthritis clinics (EACs).

Methods: Four GPs and two RNs were trained in the assessment of early in?ammatory arthritis (IA) by four visits to an EAC supervised by hospital rheumatologists. Patients referred to one of three EACs were recruited for study and assessed independently by a GP, an RN and one of six rheumatologists. Each assessor was asked to record their clinical ?ndings and whether they considered the patient to have IA. Each was then asked to judge the appropriateness of the referral according to predetermined guidelines. The rheumatologists had been shown previously to have a satisfactory level of agreement in the assessment of IA.

Results: Ninety-six patients were approached and all consented to take part in the study. In 49 cases (51%), the rheumatologist judged that the patient had IA and that the referral was appropriate. The assessments of GPs and RNs were compared with those of the rheumatologists. Levels of agreement were measured using the kappa value, where 1.0 represents total unanimity. The kappa value was
0.77 for the GPs when compared with the rheumatologists and 0.79 for the RNs. Signi?cant stiffness in the morning or after rest and objective joint swelling were the most important clinical features enabling the GPs and RNs to discriminate between IA and non-IA conditions.

Conclusion: Diagnostic triage by GPs or RNs improved the positive predictive value of referrals to an EAC with a degree of accuracy approaching that of a group of experienced rheumatologists.

The SPHERE Study: Secondary prevention of heart disease in general practice: protocol of a randomised controlled trial of tailored practice and patient care plans with parallel qualitative, economic and policy analyses. [ISRCTN24081411]

Background: The aim of the SPHERE study is to design, implement and evaluate tailored practice and personal care plans to improve the process of care and objective clinical outcomes for patients with established coronary heart disease (CHD) in general practice across two different health systems on the island of Ireland.CHD is a common cause of death and a significant cause of morbidity in Ireland. Secondary prevention has been recommended as a key strategy for reducing levels of CHD mortality and general practice has been highlighted as an ideal setting for secondary prevention initiatives. Current indications suggest that there is considerable room for improvement in the provision of secondary prevention for patients with established heart disease on the island of Ireland. The review literature recommends structured programmes with continued support and follow-up of patients; the provision of training, tailored to practice needs of access to evidence of effectiveness of secondary prevention; structured recall programmes that also take account of individual practice needs; and patient-centred consultations accompanied by attention to disease management guidelines.

Methods: SPHERE is a cluster randomised controlled trial, with practice-level randomisation to intervention and control groups, recruiting 960 patients from 48 practices in three study centres (Belfast, Dublin and Galway). Primary outcomes are blood pressure, total cholesterol, physical and mental health status (SF-12) and hospital re-admissions. The intervention takes place over two years and data is collected at baseline, one-year and two-year follow-up. Data is obtained from medical charts, consultations with practitioners, and patient postal questionnaires. The SPHERE intervention involves the implementation of a structured systematic programme of care for patients with CHD attending general practice. It is a multi-faceted intervention that has been developed to respond to barriers and solutions to optimal secondary prevention identified in preliminary qualitative research with practitioners and patients. General practitioners and practice nurses attend training sessions in facilitating behaviour change and medication prescribing guidelines for secondary prevention of CHD. Patients are invited to attend regular four-monthly consultations over two years, during which targets and goals for secondary prevention are set and reviewed. The analysis will be strengthened by economic, policy and qualitative components.

Pullout test correlations and in-place strength assessment - the European Concrete Frame Building Project

A full-scale, seven-story, reinforced concrete building frame was constructed in-place at the Building Research Establishment's Cardington Laboratory, which encompassed a range of different concrete mixtures and advanced construction techniques. This provided an opportunity to assess in-place nondestructive test methods, namely the pullout test, and more specifically the Danish version, which has been known as the Lok test, on a systematic basis during the construction of the building. It was used in conjunction with both standard and temperature-matched cube specimens to assess its practicality and accuracy under site conditions. Strength correlations were determined using linear and power function regression analysis. Strength predictions from these were found to be in very good agreement with the compressive strengths of temperature-matched cube specimens. When a general correlation is used, however, estimates for compressive strength are likely to have 95% confidence limits of around '20% of the mean value of four results.

Decision making in childhood asthma exacerbation – a clinical judgement analysis

Background: Clinical decisions which impact directly on patient safety and quality of care are made during acute asthma attacks by individual doctors on the basis of their knowledge and experience. These include administration of systemic corticosteroids (CS), oral antibiotics, and admission to hospital. Clinical judgement analysis provides a methodology for comparing decisions between practitioners with different training and experience, and improving decision making. Methods: Stepwise linear regression was used to select clinical cues based on visual analogue scale assessments of the propensity of 62 clinicians to prescribe a short course of oral CS (decision 1), a course of antibiotics (decision 2), and/or admit to hospital (decision 3) for 60 â??paperâ?? patients. Results:When compared by specialty, paediatriciansâ?? models for decision 1 were more likely to include as a cue level of alertness (54% v. 16%); for decision 2 presence of crepitations (49% v. 16%), and less likely to include inhaled CS (8% v. 40%), respiratory rate (0% v. 24%), and air entry (70% v. 100%). When compared to other grades, the models derived for decision 3 by consultants/general practitioners were more likely to include wheeze severity as a cue (39% v. 6%). Conclusions: Clinicians differed in their use of individual cues and the number included in their models. Patient safety and quality of care will benefit from clarification of decision making strategies as general learning points during medical training, in the development of guidelines and care pathways, and by clinicians developing self-awareness of their own preferences.

DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Extraction of electrode kinetic parameters from microdisc voltammetric data measured under transport conditions intermediate between steady-state convergent and transient linear diffusion as typically applies to room temperature ionic liquids

The extraction of electrode kinetic parameters for electrochemical couples in room-temperature ionic liquids (RTILs) is currently an area of considerable interest. Electrochemists typically measure electrode kinetics in the limits of either transient planar or steady-state convergent diffusion for which the voltammetic response is well understood. In this paper we develop a general method allowing the extraction of this kinetic data in the region where the diffusion is intermediate between the planar and convergent limits, such as is often encountered in RTILs using microelectrode voltammetry. A general working surface is derived, allowing the inference of Butler-Volmer standard electrochemical rate constants for the peak-to-peak potential separation in a cyclic voltammogram as a function of voltage scan rate. The method is applied to the case of the ferrocene/ferrocenium couple in [C(2)mim][N(Tf)(2)] and [C(4)mim][N(Tf)(2)].