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.

Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

Protein fold recognition without Boltzmann statistics or explicit physical basis

We present a fast method for finding optimal parameters for a low-resolution (threading) force field intended to distinguish correct from incorrect folds for a given protein sequence. In contrast to other methods, the parameterization uses information from >10(7) misfolded structures as well as a set of native sequence-structure pairs. In addition to testing the resulting force field's performance on the protein sequence threading problem, results are shown that characterize the number of parameters necessary for effective structure recognition.

Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

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.

Gaze angle: a possible mechanism of visual stress in virtual reality headsets

It is known that some Virtual Reality (VR) head-mounted displays (HMDs) can cause temporary deficits in binocular vision. On the other hand, the precise mechanism by which visual stress occurs is unclear. This paper is concerned with a potential source of visual stress that has not been previously considered with regard to VR systems: inappropriate vertical gaze angle. As vertical gaze angle is raised or lowered the 'effort' required of the binocular system also changes. The extent to which changes in vertical gaze angle alter the demands placed upon the vergence eye movement system was explored. The results suggested that visual stress may depend, in part, on vertical gaze angle. The proximity of the display screens within an HMD means that a VR headset should be in the correct vertical location for any individual user. This factor may explain some previous empirical results and has important implications for headset design. Fortuitously, a reasonably simple solution exists.

Improved one-step solid-phase extraction method for morphine, morphine-3-glucuronide, and morphine-6-glucuronide from plasma and quantitation using high-performance liquid chromatography with electrochemical detection

This communication describes an improved one-step solid-phase extraction method for the recovery of morphine (M), morphine-3-glucuronide (M3G), and morphine-6-glucuronide (M6G) from human plasma with reduced coextraction of endogenous plasma constituents, compared to that of the authors' previously reported method. The magnitude of the peak caused by endogenous plasma components in the chromatogram that eluted immediately before the retention time of M3G has been reduced (similar to 80%) significantly (p < 0.01) while achieving high extraction efficiencies for the compounds of interest, viz morphine, M6G, and M3G (93.8 +/- 2.5, 91.7 +/- 1.7, and 93.1 +/- 2.2%, respectively). Furthermore, when the improved solid-phase extraction method was used, the extraction cartridge-derived late-eluting peak (retention time 90 to 100 minutes) reported in our previous method, was no longer present in the plasma extracts. Therefore the combined effect of reducing the recovery of the endogenous components of plasma that chromatographed just before the retention time of M3G and the removal of the late-eluting, extraction cartridge-derived peak has resulted in a decrease in the chromatographic run-time to 20 minutes, thereby increasing the sample throughput by up to 100%.

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.

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.

A reproducible method for automated extraction of brain volumes from 3D human head MR images

An automated method for extracting brain volumes from three commonly acquired three-dimensional (3D) MR images (proton density, T1 weighted, and T2-weighted) of the human head is described. The procedure is divided into four levels: preprocessing, segmentation, scalp removal, and postprocessing. A user-provided reference point is the sole operator-dependent input required, The method's parameters were first optimized and then fixed and applied to 30 repeat data sets from 15 normal older adult subjects to investigate its reproducibility. Percent differences between total brain volumes (TBVs) for the subjects' repeated data sets ranged from .5% to 2.2%. We conclude that the method is both robust and reproducible and has the potential for wide application.

A novel method for tailoring the pore-opening size of MCM-41 materials

A novel pore tailoring method is proposed by which the pore-opening sizes of MCM-41 materials can be finely tuned without significant loss in pore volume and surface area.