Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Evaluation of statistical and geostatistical models of digital soil properties mapping in tropical mountain regions

Soil properties have an enormous impact on economic and environmental aspects of agricultural production. Quantitative relationships between soil properties and the factors that influence their variability are the basis of digital soil mapping. The predictive models of soil properties evaluated in this work are statistical (multiple linear regression-MLR) and geostatistical (ordinary kriging and co-kriging). The study was conducted in the municipality of Bom Jardim, RJ, using a soil database with 208 sampling points. Predictive models were evaluated for sand, silt and clay fractions, pH in water and organic carbon at six depths according to the specifications of the consortium of digital soil mapping at the global level (GlobalSoilMap). Continuous covariates and categorical predictors were used and their contributions to the model assessed. Only the environmental covariates elevation, aspect, stream power index (SPI), soil wetness index (SWI), normalized difference vegetation index (NDVI), and b3/b2 band ratio were significantly correlated with soil properties. The predictive models had a mean coefficient of determination of 0.21. Best results were obtained with the geostatistical predictive models, where the highest coefficient of determination 0.43 was associated with sand properties between 60 to 100 cm deep. The use of a sparse data set of soil properties for digital mapping can explain only part of the spatial variation of these properties. The results may be related to the sampling density and the quantity and quality of the environmental covariates and predictive models used.

Rebuild Iowa Office Transition Strategy, December 2010

The closing of the RIO does not mean that the recovery process is complete for Iowa families and communities. The recovery process will continue for many years to come, and the Rebuild Iowa Transition Strategy has been drafted to provide a comprehensive set of recommended action steps to help the state complete long-term recovery efforts while better preparing the state for future disasters. This report begins with a review of the twelve major Rebuild Iowa Advisory Commission (RIAC) recommendations which have guided RIO’s work, followed by a summary of the major accomplishments toward each recommendation. Complete, detailed information on all the work that has been accomplished toward the RIAC recommendations can be found in the RIO’s Quarterly Reports. The identification of remaining needs and issues serves as the basis for the transition strategy.

Kernel-based mapping of orographic rainfall enhancement in the Swiss Alps as detected by weatherradar

In this paper, we develop a data-driven methodology to characterize the likelihood of orographic precipitation enhancement using sequences of weather radar images and a digital elevation model (DEM). Geographical locations with topographic characteristics favorable to enforce repeatable and persistent orographic precipitation such as stationary cells, upslope rainfall enhancement, and repeated convective initiation are detected by analyzing the spatial distribution of a set of precipitation cells extracted from radar imagery. Topographic features such as terrain convexity and gradients computed from the DEM at multiple spatial scales as well as velocity fields estimated from sequences of weather radar images are used as explanatory factors to describe the occurrence of localized precipitation enhancement. The latter is represented as a binary process by defining a threshold on the number of cell occurrences at particular locations. Both two-class and one-class support vector machine classifiers are tested to separate the presumed orographic cells from the nonorographic ones in the space of contributing topographic and flow features. Site-based validation is carried out to estimate realistic generalization skills of the obtained spatial prediction models. Due to the high class separability, the decision function of the classifiers can be interpreted as a likelihood or susceptibility of orographic precipitation enhancement. The developed approach can serve as a basis for refining radar-based quantitative precipitation estimates and short-term forecasts or for generating stochastic precipitation ensembles conditioned on the local topography.

How to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index.

Cortical folding (gyrification) is determined during the first months of life, so that adverse events occurring during this period leave traces that will be identifiable at any age. As recently reviewed by Mangin and colleagues(2), several methods exist to quantify different characteristics of gyrification. For instance, sulcal morphometry can be used to measure shape descriptors such as the depth, length or indices of inter-hemispheric asymmetry(3). These geometrical properties have the advantage of being easy to interpret. However, sulcal morphometry tightly relies on the accurate identification of a given set of sulci and hence provides a fragmented description of gyrification. A more fine-grained quantification of gyrification can be achieved with curvature-based measurements, where smoothed absolute mean curvature is typically computed at thousands of points over the cortical surface(4). The curvature is however not straightforward to comprehend, as it remains unclear if there is any direct relationship between the curvedness and a biologically meaningful correlate such as cortical volume or surface. To address the diverse issues raised by the measurement of cortical folding, we previously developed an algorithm to quantify local gyrification with an exquisite spatial resolution and of simple interpretation. Our method is inspired of the Gyrification Index(5), a method originally used in comparative neuroanatomy to evaluate the cortical folding differences across species. In our implementation, which we name local Gyrification Index (lGI(1)), we measure the amount of cortex buried within the sulcal folds as compared with the amount of visible cortex in circular regions of interest. Given that the cortex grows primarily through radial expansion(6), our method was specifically designed to identify early defects of cortical development. In this article, we detail the computation of local Gyrification Index, which is now freely distributed as a part of the FreeSurfer Software (http://surfer.nmr.mgh.harvard.edu/, Martinos Center for Biomedical Imaging, Massachusetts General Hospital). FreeSurfer provides a set of automated reconstruction tools of the brain's cortical surface from structural MRI data. The cortical surface extracted in the native space of the images with sub-millimeter accuracy is then further used for the creation of an outer surface, which will serve as a basis for the lGI calculation. A circular region of interest is then delineated on the outer surface, and its corresponding region of interest on the cortical surface is identified using a matching algorithm as described in our validation study(1). This process is repeatedly iterated with largely overlapping regions of interest, resulting in cortical maps of gyrification for subsequent statistical comparisons (Fig. 1). Of note, another measurement of local gyrification with a similar inspiration was proposed by Toro and colleagues(7), where the folding index at each point is computed as the ratio of the cortical area contained in a sphere divided by the area of a disc with the same radius. The two implementations differ in that the one by Toro et al. is based on Euclidian distances and thus considers discontinuous patches of cortical area, whereas ours uses a strict geodesic algorithm and include only the continuous patch of cortical area opening at the brain surface in a circular region of interest.

Neuromodulation by oxytocin and vasopressin in the central nervous system as a basis for their rapid behavioral effects.

The last several years have seen an increasing number of studies that describe effects of oxytocin and vasopressin on the behavior of animals or humans. Studies in humans have reported behavioral changes and, through fMRI, effects on brain function. These studies are paralleled by a large number of reports, mostly in rodents, that have also demonstrated neuromodulatory effects by oxytocin and vasopressin at the circuit level in specific brain regions. It is the scope of this review to give a summary of the most recent neuromodulatory findings in rodents with the aim of providing a potential neurophysiological basis for their behavioral effects. At the same time, these findings may point to promising areas for further translational research towards human applications.

Molecular basis of leukocyte rolling on PSGL-1. Predominant role of core-2 O-glycans and of tyrosine sulfate residue 51.

Interactions between the leukocyte adhesion receptor L-selectin and P-selectin glycoprotein ligand-1 play an important role in regulating the inflammatory response by mediating leukocyte tethering and rolling on adherent leukocytes. In this study, we have examined the effect of post-translational modifications of PSGL-1 including Tyr sulfation and presentation of sialylated and fucosylated O-glycans for L-selectin binding. The functional importance of these modifications was determined by analyzing soluble L-selectin binding and leukocyte rolling on CHO cells expressing various glycoforms of PSGL-1 or mutant PSGL-1 targeted at N-terminal Thr or Tyr residues. Simultaneous expression of core-2 beta1,6-N-acetylglucosaminyltransferase and fucosyltransferase VII was required for optimal L-selectin binding to PSGL-1. Substitution of Thr-57 by Ala but not of Thr-44, strongly decreased L-selectin binding and leukocyte rolling on PSGL-1. Substitution of Tyr by Phe revealed that PSGL-1 Tyr-51 plays a predominant role in mediating L-selectin binding and leukocyte rolling whereas Tyr-48 has a minor role, an observation that contrasts with the pattern seen for the interactions between PSGL-1 and P-selectin where Tyr-48 plays a key role. Molecular modeling analysis of L-selectin and P-selectin interactions with PSGL-1 further supported these observations. Additional experiments showed that core-2 O-glycans attached to Thr-57 were also of critical importance in regulating the velocity and stability of leukocyte rolling. These observations pinpoint the structural characteristics of PSGL-1 that are required for optimal interactions with L-selectin and may be responsible for the specific kinetic and mechanical bond properties of the L-selectin-PSGL-1 adhesion receptor-counterreceptor pair.

Genetic basis and variable phenotypic expression of Kallmann syndrome: towards a unifying theory.

Idiopathic hypogonadotropic hypogonadism (IHH) is defined by absent or incomplete puberty and characterised biochemically by low levels of sex steroids, with low or inappropriately normal gonadotropin hormones. IHH is frequently accompanied by non-reproductive abnormalities, most commonly anosmia, which is present in 50-60% of cases and defines Kallmann syndrome. The understanding of IHH has undergone rapid evolution, both in respect of genetics and breadth of phenotype. Once considered in monogenic Mendelian terms, it is now more coherently understood as a complex genetic condition. Oligogenic and complex genetic-environmental interactions have now been identified, with physiological and environmental factors interacting in genetically susceptible individuals to alter the clinical course and phenotype. These potentially link IHH to ancient evolutionary pressures on the ancestral human genome.

Functional dissection of the ash2 and ash1 transcriptomes provides insights into the transcriptional basis of wing phenotypes and reveals conserved protein interactions

Background: The trithorax group (trxG) genes absent, small or homeotic discs 1 (ash1) and 2 (ash2) were isolated in a screen for mutants with abnormal imaginal discs. Mutations in either gene cause homeotic transformations but Hox genes are not their only targets. Although analysis of double mutants revealed that ash2 and ash1 mutations enhance each other's phenotypes, suggesting they are functionally related, it was shown that these proteins are subunits of distinct complexes.Results: The analysis of wing imaginal disc transcriptomes from ash2 and ash1 mutants showed that they are highly similar. Functional annotation of regulated genes using Gene Ontology allowed identification of severely affected groups of genes that could be correlated to the wing phenotypes observed. Comparison of the differentially expressed genes with those from other genome-wide analyses revealed similarities between ASH2 and Sin3A, suggesting a putative functional relationship. Coimmunoprecipitation studies and immunolocalization on polytene chromosomes demonstrated that ASH2 and Sin3A interact with HCF (host-cell factor). The results of nucleosome western blots and clonal analysis indicated that ASH2 is necessary for trimethylation of the Lys4 on histone 3 (H3K4).Conclusion: The similarity between the transcriptomes of ash2 and ash1 mutants supports a model in which the two genes act together to maintain stable states of transcription. Like in humans, both ASH2 and Sin3A bind HCF. Finally, the reduction of H3K4 trimethylation in ash2 mutants is the first evidence in Drosophila regarding the molecular function of this trxG gene.

Monitoring Disease Activity and Progression in Crohn's Disease. A Swiss Perspective on the IBD Ahead 'Optimised Monitoring' Recommendations.

BACKGROUND AND AIMS: The structured IBD Ahead 'Optimised Monitoring' programme was designed to obtain the opinion, insight and advice of gastroenterologists on optimising the monitoring of Crohn's disease activity in four settings: (1) assessment at diagnosis, (2) monitoring in symptomatic patients, (3) monitoring in asymptomatic patients, and (4) the postoperative follow-up. For each of these settings, four monitoring methods were discussed: (a) symptom assessment, (b) endoscopy, (c) laboratory markers, and (d) imaging. Based on literature search and expert opinion compiled during an international consensus meeting, recommendations were given to answer the question 'which diagnostic method, when, and how often'. The International IBD Ahead Expert Panel advised to tailor this guidance to the healthcare system and the special prerequisites of each country. The IBD Ahead Swiss National Steering Committee proposes best-practice recommendations adapted for Switzerland. METHODS: The IBD Ahead Steering Committee identified key questions and provided the Swiss Expert Panel with a structured literature research. The expert panel agreed on a set of statements. During an international expert meeting the consolidated outcome of the national meetings was merged into final statements agreed by the participating International and National Steering Committee members - the IBD Ahead 'Optimized Monitoring' Consensus. RESULTS: A systematic assessment of symptoms, endoscopy findings, and laboratory markers with special emphasis on faecal calprotectin is deemed necessary even in symptom-free patients. The choice of recommended imaging methods is adapted to the specific situation in Switzerland and highlights the importance of ultrasonography and magnetic resonance imaging besides endoscopy. CONCLUSION: The recommendations stress the importance of monitoring disease activity on a regular basis and by objective parameters, such as faecal calprotectin and endoscopy with detailed documentation of findings. Physicians should not rely on symptoms only and adapt the monitoring schedule and choice of options to individual situations. © 2014 S. Karger AG, Basel.

Iowa DOT Quick Tour, February 2012

The Iowa DOT, created by the Iowa Legislature, began operating July 1, 1975. The newly formed Department of Transportation had a staff of 4,568 full-time equivalents (FTEs). These employees and their responsibilities came from the Highway Commission, Reciprocity Board, Aeronautics Commission, and parts of the Commerce Commission, Department of Revenue and Finance, Department of Public Safety and the Energy Policy Council. The Iowa DOT continually looks at ways to become more efficient in serving our customers. As a result, the current workforce equals approximately 2,818 full-time employees. The seven transportation commissioners set policy for the department and are appointed by the governor, with no more than four from one political party. Membership on the commission is gender-balanced. Commissioners are confirmed by the Iowa Senate and serve on a staggered basis for four-year terms. Commission meetings are usually held monthly in Ames. The meetings are open to the public and follow a printed agenda.

Apprentissage automatique de la morphologie : le cas des structures racine-schème

RÉSUMÉ Cette thèse porte sur le développement de méthodes algorithmiques pour découvrir automatiquement la structure morphologique des mots d'un corpus. On considère en particulier le cas des langues s'approchant du type introflexionnel, comme l'arabe ou l'hébreu. La tradition linguistique décrit la morphologie de ces langues en termes d'unités discontinues : les racines consonantiques et les schèmes vocaliques. Ce genre de structure constitue un défi pour les systèmes actuels d'apprentissage automatique, qui opèrent généralement avec des unités continues. La stratégie adoptée ici consiste à traiter le problème comme une séquence de deux sous-problèmes. Le premier est d'ordre phonologique : il s'agit de diviser les symboles (phonèmes, lettres) du corpus en deux groupes correspondant autant que possible aux consonnes et voyelles phonétiques. Le second est de nature morphologique et repose sur les résultats du premier : il s'agit d'établir l'inventaire des racines et schèmes du corpus et de déterminer leurs règles de combinaison. On examine la portée et les limites d'une approche basée sur deux hypothèses : (i) la distinction entre consonnes et voyelles peut être inférée sur la base de leur tendance à alterner dans la chaîne parlée; (ii) les racines et les schèmes peuvent être identifiés respectivement aux séquences de consonnes et voyelles découvertes précédemment. L'algorithme proposé utilise une méthode purement distributionnelle pour partitionner les symboles du corpus. Puis il applique des principes analogiques pour identifier un ensemble de candidats sérieux au titre de racine ou de schème, et pour élargir progressivement cet ensemble. Cette extension est soumise à une procédure d'évaluation basée sur le principe de la longueur de description minimale, dans- l'esprit de LINGUISTICA (Goldsmith, 2001). L'algorithme est implémenté sous la forme d'un programme informatique nommé ARABICA, et évalué sur un corpus de noms arabes, du point de vue de sa capacité à décrire le système du pluriel. Cette étude montre que des structures linguistiques complexes peuvent être découvertes en ne faisant qu'un minimum d'hypothèses a priori sur les phénomènes considérés. Elle illustre la synergie possible entre des mécanismes d'apprentissage portant sur des niveaux de description linguistique distincts, et cherche à déterminer quand et pourquoi cette coopération échoue. Elle conclut que la tension entre l'universalité de la distinction consonnes-voyelles et la spécificité de la structuration racine-schème est cruciale pour expliquer les forces et les faiblesses d'une telle approche. ABSTRACT This dissertation is concerned with the development of algorithmic methods for the unsupervised learning of natural language morphology, using a symbolically transcribed wordlist. It focuses on the case of languages approaching the introflectional type, such as Arabic or Hebrew. The morphology of such languages is traditionally described in terms of discontinuous units: consonantal roots and vocalic patterns. Inferring this kind of structure is a challenging task for current unsupervised learning systems, which generally operate with continuous units. In this study, the problem of learning root-and-pattern morphology is divided into a phonological and a morphological subproblem. The phonological component of the analysis seeks to partition the symbols of a corpus (phonemes, letters) into two subsets that correspond well with the phonetic definition of consonants and vowels; building around this result, the morphological component attempts to establish the list of roots and patterns in the corpus, and to infer the rules that govern their combinations. We assess the extent to which this can be done on the basis of two hypotheses: (i) the distinction between consonants and vowels can be learned by observing their tendency to alternate in speech; (ii) roots and patterns can be identified as sequences of the previously discovered consonants and vowels respectively. The proposed algorithm uses a purely distributional method for partitioning symbols. Then it applies analogical principles to identify a preliminary set of reliable roots and patterns, and gradually enlarge it. This extension process is guided by an evaluation procedure based on the minimum description length principle, in line with the approach to morphological learning embodied in LINGUISTICA (Goldsmith, 2001). The algorithm is implemented as a computer program named ARABICA; it is evaluated with regard to its ability to account for the system of plural formation in a corpus of Arabic nouns. This thesis shows that complex linguistic structures can be discovered without recourse to a rich set of a priori hypotheses about the phenomena under consideration. It illustrates the possible synergy between learning mechanisms operating at distinct levels of linguistic description, and attempts to determine where and why such a cooperation fails. It concludes that the tension between the universality of the consonant-vowel distinction and the specificity of root-and-pattern structure is crucial for understanding the advantages and weaknesses of this approach.

The set of undominated imputations and the core: an axiomatic approach

This paper provides an axiomatic framework to compare the D-core (the set of undominatedimputations) and the core of a cooperative game with transferable utility. Theorem1 states that the D-core is the only solution satisfying projection consistency, reasonableness (from above), (*)-antimonotonicity, and modularity. Theorem 2 characterizes the core replacing (*)-antimonotonicity by antimonotonicity. Moreover, these axioms alsocharacterize the core on the domain of convex games, totally balanced games, balancedgames, and superadditive games

Coalitionally monotonic set-solutions for cooperative TU games

A static comparative study on set-solutions for cooperative TU games is carried out. The analysis focuses on studying the compatibility between two classical and reasonable properties introduced by Young (1985) in the context of single valued solutions, namely core-selection and coalitional monotonicity. As the main result, it is showed that coalitional monotonicity is not only incompatible with the core-selection property but also with the bargaining-selection property. This new impossibility result reinforces the tradeoff between these kinds of interesting and intuitive economic properties. Positive results about compatibility between desirable economic properties are given replacing the core selection requirement by the core-extension property.