Verifying data refinements using a model checker

In this paper, we consider how refinements between state-based specifications (e.g., written in Z) can be checked by use of a model checker. Specifically, we are interested in the verification of downward and upward simulations which are the standard approach to verifying refinements in state-based notations. We show how downward and upward simulations can be checked using existing temporal logic model checkers. In particular, we show how the branching time temporal logic CTL can be used to encode the standard simulation conditions. We do this for both a blocking, or guarded, interpretation of operations (often used when specifying reactive systems) as well as the more common non-blocking interpretation of operations used in many state-based specification languages (for modelling sequential systems). The approach is general enough to use with any state-based specification language, and we illustrate how refinements between Z specifications can be checked using the SAL CTL model checker using a small example.

Pore characterization of carbonaceous materials by DFT and GCMC simulations: A review

A review is given of the pore characterization of carbonaceous materials, including activated carbon, carbon fibres, carbon nanotubes, etc., using adsorption techniques. Since the pores of carbon media are mostly of molecular dimensions, the appropriate modem tools for the analysis of adsorption isotherms are grand canonical Monte Carlo (GCMC) simulations and density functional theory (DFT). These techniques are presented and applications of such tools in the derivation of pore-size distribution highlighted.

CFD simulations of hydrocyclones with an air core - Comparison between large eddy simulations and a second moment closure

CFD simulations of the 75 mm, hydrocyclone of Hsieh (1988) have been conducted using Fluent TM. The simulations used 3-dimensional body fitted grids. The simulations were two phase simulations where the air core was resolved using the mixture (Manninen et al., 1996) and VOF (Hirt and Nichols, 1981) models. Velocity predictions from large eddy simulations (LES), using the Smagorinsky-Lilly sub grid scale model (Smagorinsky, 1963; Lilly, 1966) and RANS simulations using the differential Reynolds stress turbulence model (Launder et al., 1975) were compared with Hsieh's experimental velocity data. The LES simulations gave very good agreement with Hsieh's data but required very fine grids to predict the velocities correctly in the bottom of the apex. The DRSM/RANS simulations under predicted tangential velocities, and there was little difference between the velocity predictions using the linear (Launder, 1989) and quadratic (Speziale et al., 1991) pressure strain models. Velocity predictions using the DRSM turbulence model and the linear pressure strain model could be improved by adjusting the pressure strain model constants.

Toposequential effects on water balance and productivity in rainfed lowland rice ecosystem in Southern Laos

In the Mekong region, most paddies in rainfed lowland rice (Oryza sativa L.) lie in a sequence on gentle sloping land, and grain yield (GY) often depends on the toposequence position. There is, however, lack of information on toposequential effects on field water supply in rainfed lowland rice and how that influences GY. A total of eight field experiments were carried out on sandy, coarse-textured soils in Southern Laos (Champassak Province and Savannakhet Province) over three wet seasons (2000-2002). Components of the water balance, including downward and lateral water movement (D and L, respectively), were quantified at three different positions along toposequences (top, middle and bottom). GY, days-to-flower (DTF) and rainfall were measured, and the water productivity (WP) was determined. In most experiments, standing water disappeared first in the top position and gradually in lower positions. This was associated with the observation that when there was standing water in the field, the higher position had larger D in both the provinces and also larger L in Champassak Province. However, in one experiment, water loss appeared later in the higher position, as the result of lower L, apparently due to some water inputs other than rainfall occurring at this position. Despite larger D plus L at the top position, seasonal sum of D and L were not much affected by the toposequence position, as the daily rate of D plus L became minimal when the standing water was lost earlier in the top position. Lower GY was associated with earlier disappearance of standing water from the field. Relatively low GY was expected in the top toposequence position. This was clearly shown in the toposequence of Phonthong, Champassak Province, as the timing of standing water disappearance relative to flowering was earlier in the top position. Variation in GY across the toposequence positions was coupled with the WP variation, and both GY and WP tended to decline with increased DTF. Therefore, variation in productivity of rainfed lowland rice across toposequence positions depends mainly on the field water status around flowering time. (c) 2005 Elsevier B.V. All rights reserved.

The v-MFG test: Investigating maternal, offspring and maternal-fetal genetic incompatibility effects on disease and viability

The MFG test is a family-based association test that detects genetic effects contributing to disease in offspring, including offspring allelic effects, maternal allelic effects and MFG incompatibility effects. Like many other family-based association tests, it assumes that the offspring survival and the offspring-parent genotypes are conditionally independent provided the offspring is affected. However, when the putative disease-increasing locus can affect another competing phenotype, for example, offspring viability, the conditional independence assumption fails and these tests could lead to incorrect conclusions regarding the role of the gene in disease. We propose the v-MFG test to adjust for the genetic effects on one phenotype, e.g., viability, when testing the effects of that locus on another phenotype, e.g., disease. Using genotype data from nuclear families containing parents and at least one affected offspring, the v-MFG test models the distribution of family genotypes conditional on offspring phenotypes. It simultaneously estimates genetic effects on two phenotypes, viability and disease. Simulations show that the v-MFG test produces accurate genetic effect estimates on disease as well as on viability under several different scenarios. It generates accurate type-I error rates and provides adequate power with moderate sample sizes to detect genetic effects on disease risk when viability is reduced. We demonstrate the v-MFG test with HLA-DRB1 data from study participants with rheumatoid arthritis (RA) and their parents, we show that the v-MFG test successfully detects an MFG incompatibility effect on RA while simultaneously adjusting for a possible viability loss.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Toxin insights into nicotinic acetylcholine receptors

Venomous species have evolved cocktails of bioactive peptides to facilitate prey capture. Given their often exquisite potency and target selectivity, venom peptides provide unique biochemical tools for probing the function of membrane proteins at the molecular level. in the field of the nicotinic acetylcholine receptors (nAChRs), the subtype specific snake alpha-neurotoxins and cone snail alpha-conotoxins have been widely used to probe receptor structure and function in native tissues and recombinant systems. However, only recently has it been possible to generate an accurate molecular view of these nAChR-toxin interactions. Crystal structures of AChBP, a homologue of the nAChR ligand binding domain, have now been solved in complex with alpha-cobratoxin, alpha-conotoxin PnIA and alpha-conotoxin Iml. The orientation of all three toxins in the ACh binding site confirms many of the predictions obtained from mutagenesis and docking simulations on homology models of mammalian nAChR. The precise understanding of the molecular determinants of these complexes is expected to contribute to the development of more selective nAChR modulators. In this commentary, we review the structural data on nAChR-toxin interactions and discuss their implications for the design of novel ligands acting at the nAChR. (c) 2006 Elsevier Inc. All rights reserved.

Robustness of mathematical models for biological systems

The robustness of mathematical models for biological systems is studied by sensitivity analysis and stochastic simulations. Using a neural network model with three genes as the test problem, we study robustness properties of synthesis and degradation processes. For single parameter robustness, sensitivity analysis techniques are applied for studying parameter variations and stochastic simulations are used for investigating the impact of external noise. Results of sensitivity analysis are consistent with those obtained by stochastic simulations. Stochastic models with external noise can be used for studying the robustness not only to external noise but also to parameter variations. For external noise we also use stochastic models to study the robustness of the function of each gene and that of the system.