An improved novel method for solving nonlinear boundary value problems

A modified formula for the integral transform of a nonlinear function is proposed for a class of nonlinear boundary value problems. The technique presented in this paper results in analytical solutions. Iterations and initial guess, which are needed in other techniques, are not required in this novel technique. The analytical solutions are found to agree surprisingly well with the numerically exact solutions for two examples of power law reaction and Langmuir-Hinshelwood reaction in a catalyst pellet.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

A hybrid planning method for transmission networks in a deregulated environment

The reconstruction of power industries has brought fundamental changes to both power system operation and planning. This paper presents a new planning method using multi-objective optimization (MOOP) technique, as well as human knowledge, to expand the transmission network in open access schemes. The method starts with a candidate pool of feasible expansion plans. Consequent selection of the best candidates is carried out through a MOOP approach, of which multiple objectives are tackled simultaneously, aiming at integrating the market operation and planning as one unified process in context of deregulated system. Human knowledge has been applied in both stages to ensure the selection with practical engineering and management concerns. The expansion plan from MOOP is assessed by reliability criteria before it is finalized. The proposed method has been tested with the IEEE 14-bus system and relevant analyses and discussions have been presented.

A Step-by-Step Approach to Transit Oriented Development Project Delivery

A major challenge in successfully implementing transit-oriented development (TOD) is having a robust process that ensures effective appraisal, initiation and delivery of multi-stakeholder TOD projects. A step-by step project development process can assist in the methodic design, evaluation, and initiation of TOD projects. Successful TOD requires attention to transit, mixed-use development and public space. Brisbane, Australia provides a case-study where recent planning policies and infrastructure documents have laid a foundation for TOD, but where barriers lie in precinct level planning and project implementation. In this context and perhaps in others, the research effort needs to shift toward identification of appropriate project processes and strategies. This paper presents the outcomes of research conducted to date. Drawing on the mainstream approach to project development and financial evaluation for property projects, key steps for potential use in successful delivery of TOD projects have been identified, including: establish the framework; location selection; precinct context review; preliminary precinct design; the initial financial viability study; the decision stage; establishment of project structure; land acquisition; development application; and project delivery. The appropriateness of this mainstream development and appraisal process will be tested through stakeholder research, and the proposed process will then be refined for adoption in TOD projects. It is suggested that the criteria for successful TOD should be broadened beyond financial concerns in order to deliver public sector support for project initiation.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Bifurcation in growth patterns for arrays of parallel Griffith, edge and sliding cracks

This paper presents the recent finding by Muhlhaus et al [1] that bifurcation of crack growth patterns exists for arrays of two-dimensional cracks. This bifurcation is a result of the nonlinear effect due to crack interaction, which is, in the present analysis, approximated by the dipole asymptotic or pseudo-traction method. The nonlinear parameter for the problem is the crack length/ spacing ratio lambda = a/h. For parallel and edge crack arrays under far field tension, uniform crack growth patterns (all cracks having same size) yield to nonuniform crack growth patterns (i.e. bifurcation) if lambda is larger than a critical value lambda(cr) (note that such bifurcation is not found for collinear crack arrays). For parallel and edge crack arrays respectively, the value of lambda(cr) decreases monotonically from (2/9)(1/2) and (2/15.096)(1/2) for arrays of 2 cracks, to (2/3)(1/2)/pi and (2/5.032)(1/2)/pi for infinite arrays of cracks. The critical parameter lambda(cr) is calculated numerically for arrays of up to 100 cracks, whilst discrete Fourier transform is used to obtain the exact solution of lambda(cr) for infinite crack arrays. For geomaterials, bifurcation can also occurs when array of sliding cracks are under compression.

A PCR-based method of screening individuals of all ages, from neonates to the elderly, for familial hyperaldosteronism type I

Aim: Unless specifically treated (glucocorticoids in low doses), Familial Hyperaldosteronism Type I(FH-I) may result in early death from stroke. We report the successful application of a rapid, polymerase chain reaction (PCR)-based method of detecting the 'hybrid' 11 beta-hydroxylase (11 beta-OHase)/aldosterone synthase (AS) gene as a screening test for FH-I. Methods: 'Long-PCR' was used to amplify, concurrently, a 4 kb fragment of AS gene (both primers AS-specific) and a 4 kb fragment of the hybrid gene (5' primer 11 beta-OHase-specific, 3'primer AS-specific) from DNA extracted from blood either collected locally or transported from elsewhere. Sample collection and transport were straightforward. This 4 kb fragment contains all the currently recognised hybrid gene 'crossover' points. Results: Within a single family, long-PCR identified all 21 individuals known to have FH-I. Hypertension was corrected in all 11 treated with glucocorticoids. Nine with normal blood pressure are being closely followed for development of hypertension. Long-PCR cord blood analysis excluded FH-I in three neonates born to affected individuals. Long-PCR newly identified two other affected families: (1) a female (60 years) with a personal and family history of stroke and her normotensive daughter (40 years), and (2) a female (51 years) previously treated for primary aldosteronism with amiloride, her two hypertensive sons (14 and 16 years) and her hypertensive mother (78 years). No false negative or false positive results have yet been encountered. At least seven other centres have successfully performed this test. Conclusion: Long-PCR is a reliable method of screening individuals of all ages for FH-I.

Electrical impedance particle size method (Coulter Counter) detects the large fat globules in poorly homogenised UHT processed milk

The large fat globules that can be present in UHT milk due to inadequate homogenisation cause a cream layer to form that limits the shelf life of UHT milk. Four different particle size measurement techniques were used to measure the size of fat globules in poorly homogenised UHT milk processed in a UHT pilot plant. The thickness of the cream layer that formed during storage was negatively correlated with homogenisation pressure. It was positively correlated with the mass mean diameter and the percentage volume of particles between 1.5 and 2 mu m diameter, as determined by laser light scattering using the Malvern Mastersizer. Also, the thickness of the cream layer was positively correlated with the volume mode diameter and the percentage volume of particles between 1.5 and 2 mu m diameter, as determined by electrical impedance using the Coulter Counter. The cream layer thickness did not correlate significantly with the Coulter Counter measurements of volume mean diameter, or volume percentages of particles between 2 and 5 mu m or 5 and 10 mu m diameter. Spectroturbidimetry (Emulsion Quality Analyser) and light microscopy analyses were found to be unsuitable for assessing the size of the fat particles. This study suggests that the fat globule size distribution as determined by the electrical impedance method (Coulter Counter) is the most useful for determining the efficiency of homogenisation and therefore for predicting the stability of the fat emulsion in UHT milk during storage.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Speciation of arsenic metabolites in the urine of occupational workers and experimental rats using an optimised hydride cold-trapping method

A hydride cold-trapping technique was developed and optimised for the measurement of urinary arsenic metabolites. The analytical precision of the method was found to be 6.1, 4.0 and 4.8% (n = 5) for inorganic arsenic (As-i), monomethylarsonate (MMA) and dimethylarsinate (DMA), respectively, with recoveries close to 100%, The detection limits were 1.0, 1.3 and 3 ng for As-i, MMA and DMA, respectively. The method was then used to analyse urine samples obtained from three groups of workers for occupational exposure in three companies where copper chrome arsenate was used for timber treatment. The results were compared with those for a normal control group of laboratory workers. Arsenic and its metabolites were also measured in experimental rats given 5 mg As kg(-1) body mass by oral gavage in the form of sodium arsenite, calcium arsenite or sodium arsenate. Occupational workers showed a significantly higher excretion of As-i, Up to two fold increases of urinary As-i excretion in rats compared with control rats were also observed in animals dosed with various forms of arsenicals. The method is suitable for the measurement of arsenic metabolites in urine of both humans and experimental animals.

A novel method for selection of chymotrypsin inhibitors from a phage peptide library

A novel screening strategy has been developed for the identification of alpha-chymotrypsin inhibitors from a phage peptide library. In this strategy, the standard affinity selection protocol was modified by adding a proteolytic cleavage period to avoid recovery of alpha-chymotrypsin substrates. After four cycles of selection and further activity assay, a group of related peptides were identified by DNA sequencing. These peptides share a consensus sequence motif as (S/T)RVPR(R/H). Then, a corresponding short peptide (Ac-ASRVPRRG-NH2) was synthesized chemically and proved to be an inhibitor of alpha-chymotrypsin. The present work provides a useful way for searching proteinase inhibitors without detailed knowledge of the molecular structure.

Solid-phase extraction method with high-performance liquid chromatography and electrochemical detection for the quantitative analysis of oxycodone in human plasma

A sensitive and reproducible solid-phase extraction (SPE) method for the quantification of oxycodone in human plasma was developed. Varian Certify SPE cartridges containing both C-8 and benzoic acid functional groups were the most suitable for the extraction of oxycodone and codeine (internal standard), with consistently high (greater than or equal to 80%) and reproducible recoveries. The elution mobile phase consisted of 1.2 ml of butyl chloride-isopropanol (80:20, v/v) containing 2% ammonia. The quantification limit for oxycodone was 5.3 pmol on-column. Within-day and inter-day coefficients of variation were 1.2% and 6.8% respectively for 284 nM oxycodone and 9.5% and 6.2% respectively for 28.4 nM oxycodone using 0.5-ml plasma aliquots. (C) 1998 Elsevier Science BN. All rights reserved.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.