Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Heterogeneity of human monocytes: an optimized four-color flow cytometry protocol for analysis of monocyte subsets.

Monocytes are central mediators in the development of atherosclerotic plaques. They circulate in blood and eventually migrate into tissue including the vessel wall where they give rise to macrophages and dendritic cells. The existence of monocyte subsets with distinct roles in homeostasis and inflammation suggests specialization of function. These subsets are identified based on expression of the CD14 and CD16 markers. Routinely applicable protocols remain elusive, however. Here, we present an optimized four-color flow cytometry protocol for analysis of human blood monocyte subsets using a specific PE-Cy5-conjugated monoclonal antibody (mAb) to HLA-DR, a PE-Cy7-conjugated mAb to CD14, a FITC-conjugated mAb to CD16, and PE-conjugated mAbs to additional markers relevant to monocyte function. Classical CD14(+)CD16(-) monocytes (here termed "Mo1" subset) expressed high CCR2, CD36, CD64, and CD62L, but low CX(3)CR1, whereas "nonclassical" CD14(lo)CD16(+) monocytes (Mo3) essentially showed the inverse expression pattern. CD14(+)CD16(+) monocytes (Mo2) expressed high HLA-DR, CD36, and CD64. In patients with stable coronary artery disease (n = 13), classical monocytes were decreased, whereas "nonclassical" monocytes were increased 90% compared with healthy subjects with angiographically normal coronary arteries (n = 14). Classical monocytes from CAD patients expressed higher CX(3)CR1 and CCR2 than controls. Thus, stable CAD is associated with expansion of the nonclassical monocyte subset and increased expression of inflammatory markers on monocytes. Flow cytometric analysis of monocyte subsets and marker expression may provide valuable information on vascular inflammation. This may translate into the identification of monocyte subsets as selective therapeutic targets, thus avoiding adverse events associated with indiscriminate monocyte inhibition.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Increased flow of fatty acids toward beta-oxidation in developing seeds of Arabidopsis deficient in diacylglycerol acyltransferase activity or synthesizing medium-chain-length fatty acids.

Synthesis of polyhydroxyalkanoates (PHAs) from intermediates of fatty acid beta-oxidation was used as a tool to study fatty acid degradation in developing seeds of Arabidopsis. Transgenic plants expressing a peroxisomal PHA synthase under the control of a napin promoter accumulated PHA in developing seeds to a final level of 0. 06 mg g(-1) dry weight. In plants co-expressing a plastidial acyl-acyl carrier protein thioesterase from Cuphea lanceolata and a peroxisomal PHA synthase, approximately 18-fold more PHA accumulated in developing seeds. The proportion of 3-hydroxydecanoic acid monomer in the PHA was strongly increased, indicating a large flow of capric acid toward beta-oxidation. Furthermore, expression of the peroxisomal PHA synthase in an Arabidopsis mutant deficient in the enzyme diacylglycerol acyltransferase resulted in a 10-fold increase in PHA accumulation in developing seeds. These data indicate that plants can respond to the inadequate incorporation of fatty acids into triacylglycerides by recycling the fatty acids via beta-oxidation and that a considerable flow toward beta-oxidation can occur even in a plant tissue primarily devoted to the accumulation of storage lipids.

Oligoclonal expansions in the CD8(+)CD28(-) T cells largely explain the shorter telomeres detected in this subset: analysis by flow FISH.

We have previously reported that CD8(+)CD28(-) T cells have relatively shorter telomeres compared with CD8(+)CD28(+) T cells. Oligoclonal expansion is a common feature of CD8(+) T cells in human peripheral blood, and these expansions predominantly occur in the CD57(+)/CD28(-) population. We studied the telomere length in subsets of CD8(+) T cells using quantitative fluorescence in situ hybridization and flow cytometry (flow FISH). Our results confirm that CD8(+)CD28(-) T cells have shorter telomeres as compared with their CD28(+) counterpart cells. In addition, the oligoclonally expanded cells within the CD8(+)CD28(-) T cell subset generally have even shorter telomeres than the CD28(-) subset as a whole. We conclude that the presence of clonal expansions in the CD8(+)CD28(-) T cell population largely explain the shorter telomeres in this subset. These clonally expanded CD8(+)CD28(-) T cells generally have characteristics of terminally differentiated effector cells. Nevertheless, there is considerable individual variation in the degree of telomere shortening in these cells, which may reflect host genetic factors as well as the type and timing of the antigenic exposure.

Fluid flow regulates stromal cell organization and CCL21 expression in a tissue-engineered lymph node microenvironment.

In the paracortex of the lymph node (LN), T zone fibroblastic reticular cells (TRCs) orchestrate an immune response by guiding lymphocyte migration both physically, by creating three-dimensional (3D) cell networks, and chemically, by secreting the chemokines CCL19 and CCL21 that direct interactions between CCR7-expressing cells, including mature dendritic cells and naive T cells. TRCs also enwrap matrix-based conduits that transport fluid from the subcapsular sinus to high endothelial venules, and fluid flow through the draining LN rapidly increases upon tissue injury or inflammation. To determine whether fluid flow affects TRC organization or function within a 3D network, we regenerated the 3D LN T zone stromal network by culturing murine TRC clones within a macroporous polyurethane scaffold containing type I collagen and Matrigel and applying slow interstitial flow (1-23 microm/min). We show that the 3D environment and slow interstitial flow are important regulators of TRC morphology, organization, and CCL21 secretion. Without flow, CCL21 expression could not be detected. Furthermore, when flow through the LN was blocked in mice in vivo, CCL21 gene expression was down-regulated within 2 h. These results highlight the importance of lymph flow as a homeostatic regulator of constitutive TRC activity and introduce the concept that increased lymph flow may act as an early inflammatory cue to enhance CCL21 expression by TRCs, thereby ensuring efficient immune cell trafficking, lymph sampling, and immune response induction.

Cytostatic lung perfusion results in heterogeneous spatial regional blood flow and drug distribution: evaluation of different cytostatic lung perfusion techniques in a porcine model.

OBJECTIVES: Comparison of doxorubicin uptake, leakage and spatial regional blood flow, and drug distribution was made for antegrade, retrograde, combined antegrade and retrograde isolated lung perfusion, and pulmonary artery infusion by endovascular inflow occlusion (blood flow occlusion), as opposed to intravenous administration in a porcine model. METHODS: White pigs underwent single-pass lung perfusion with doxorubicin (320 mug/mL), labeled 99mTc-microspheres, and Indian ink. Visual assessment of the ink distribution and perfusion scintigraphy of the perfused lung was performed. 99mTc activity and doxorubicin levels were measured by gamma counting and high-performance liquid chromatography on 15 tissue samples from each perfused lung at predetermined localizations. RESULTS: Overall doxorubicin uptake in the perfused lung was significantly higher (P = .001) and the plasma concentration was significantly lower (P < .0001) after all isolated lung perfusion techniques, compared with intravenous administration, without differences between them. Pulmonary artery infusion (blood flow occlusion) showed an equally high doxorubicin uptake in the perfused lung but a higher systemic leakage than surgical isolated lung perfusion (P < .0001). The geometric coefficients of variation of the doxorubicin lung tissue levels were 175%, 279%, 226%, and 151% for antegrade, retrograde, combined antegrade and retrograde isolated lung perfusion, and pulmonary artery infusion by endovascular inflow occlusion (blood flow occlusion), respectively, compared with 51% for intravenous administration (P = .09). 99mTc activity measurements of the samples paralleled the doxorubicin level measurements, indicating a trend to a more heterogeneous spatial regional blood flow and drug distribution after isolated lung perfusion and blood flow occlusion compared with intravenous administration. CONCLUSIONS: Cytostatic lung perfusion results in a high overall doxorubicin uptake, which is, however, heterogeneously distributed within the perfused lung.

Debris flow hazard modelling on medium scale: Valtellina di Tirano, Italy

Debris flow hazard modelling at medium (regional) scale has been subject of various studies in recent years. In this study, hazard zonation was carried out, incorporating information about debris flow initiation probability (spatial and temporal), and the delimitation of the potential runout areas. Debris flow hazard zonation was carried out in the area of the Consortium of Mountain Municipalities of Valtellina di Tirano (Central Alps, Italy). The complexity of the phenomenon, the scale of the study, the variability of local conditioning factors, and the lacking data limited the use of process-based models for the runout zone delimitation. Firstly, a map of hazard initiation probabilities was prepared for the study area, based on the available susceptibility zoning information, and the analysis of two sets of aerial photographs for the temporal probability estimation. Afterwards, the hazard initiation map was used as one of the inputs for an empirical GIS-based model (Flow-R), developed at the University of Lausanne (Switzerland). An estimation of the debris flow magnitude was neglected as the main aim of the analysis was to prepare a debris flow hazard map at medium scale. A digital elevation model, with a 10 m resolution, was used together with landuse, geology and debris flow hazard initiation maps as inputs of the Flow-R model to restrict potential areas within each hazard initiation probability class to locations where debris flows are most likely to initiate. Afterwards, runout areas were calculated using multiple flow direction and energy based algorithms. Maximum probable runout zones were calibrated using documented past events and aerial photographs. Finally, two debris flow hazard maps were prepared. The first simply delimits five hazard zones, while the second incorporates the information about debris flow spreading direction probabilities, showing areas more likely to be affected by future debris flows. Limitations of the modelling arise mainly from the models applied and analysis scale, which are neglecting local controlling factors of debris flow hazard. The presented approach of debris flow hazard analysis, associating automatic detection of the source areas and a simple assessment of the debris flow spreading, provided results for consequent hazard and risk studies. However, for the validation and transferability of the parameters and results to other study areas, more testing is needed.

Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale: some case studies

Every year, debris flows cause huge damage in mountainous areas. Due to population pressure in hazardous zones, the socio-economic impact is much higher than in the past. Therefore, the development of indicative susceptibility hazard maps is of primary importance, particularly in developing countries. However, the complexity of the phenomenon and the variability of local controlling factors limit the use of processbased models for a first assessment. A debris flow model has been developed for regional susceptibility assessments using digital elevation model (DEM) with a GIS-based approach.. The automatic identification of source areas and the estimation of debris flow spreading, based on GIS tools, provide a substantial basis for a preliminary susceptibility assessment at a regional scale. One of the main advantages of this model is its workability. In fact, everything is open to the user, from the data choice to the selection of the algorithms and their parameters. The Flow-R model was tested in three different contexts: two in Switzerland and one in Pakistan, for indicative susceptibility hazard mapping. It was shown that the quality of the DEM is the most important parameter to obtain reliable results for propagation, but also to identify the potential debris flows sources.

A flow cytometry based oligotrophic pollutant exposure test to detect bacterial growth inhibition and cell injury.

Toxicity of chemical pollutants in aquatic environments is often addressed by assays that inquire reproductive inhibition of test microorganisms, such as algae or bacteria. Those tests, however, assess growth of populations as a whole via macroscopic methods such as culture turbidity or colony-forming units. Here we use flow cytometry to interrogate the fate of individual cells in low-density populations of the bacterium Pseudomonas fluorescens SV3 exposed or not under oligotrophic conditions to a number of common pollutants, some of which derive from oil contamination. Cells were stained at regular time intervals during the exposure assay with fluorescent dyes that detect membrane injury (i.e., live-dead assay). Reduction of population growth rates was observed upon toxicant insult and depended on the type of toxicant. Modeling and cell staining indicate that population growth rate decrease is a combined effect of an increased number of injured cells that may or may not multiply, and live cells dividing at normal growth rates. The oligotrophic assay concept presented here could be a useful complement for existing biomarker assays in compliance with new regulations on chemical effect studies or, more specifically, for judging recovery after exposure to fluctuating toxicant conditions.

Using flow cytometry for in situ monitoring of antimicrobial compound production in the biocontrol bacterium Pseudomonas fluorescens CHA0.

Pseudomonas fluorescens strain CHA0 is able to protect plants against a variety of pathogens, notably by producing the two antimicrobial compounds 2,4-diacetylphloroglucinol (DAPG) and pyoluteorin (PLT). The regulation of the expression of these compounds is affected by many biotic factors, such as fungal pathogens, rhizosphere bacteria as well as plant species. Therefore, the influence of some plant phenolic compounds on the expression of DAPG and PLT biosynthetic genes has been tested using GFP-based reporter, monitored by standard fluometry and flow cytometry. In situ experiments were also performed with cucumber plants. We found that several plant metabolites such as IAA and umbelliferone are able to modify significantly the expression of DAPG and PLT. The use of flow cytometry with autofluorescents proteins seems to be a promising method to study rhizobacteria-plant interactions.

Artificial neural network study of whole-cell bacterial bioreporter response determined using fluorescence flow cytometry.

Genetically engineered bioreporters are an excellent complement to traditional methods of chemical analysis. The application of fluorescence flow cytometry to detection of bioreporter response enables rapid and efficient characterization of bacterial bioreporter population response on a single-cell basis. In the present study, intrapopulation response variability was used to obtain higher analytical sensitivity and precision. We have analyzed flow cytometric data for an arsenic-sensitive bacterial bioreporter using an artificial neural network-based adaptive clustering approach (a single-layer perceptron model). Results for this approach are far superior to other methods that we have applied to this fluorescent bioreporter (e.g., the arsenic detection limit is 0.01 microM, substantially lower than for other detection methods/algorithms). The approach is highly efficient computationally and can be implemented on a real-time basis, thus having potential for future development of high-throughput screening applications.

Antigenic heterogeneity of clones and subclones from human melanoma cell lines demonstrated by a panel of monoclonal antibodies and flow microfluorometry analysis.

Cells from two melanoma cell lines, Me43 and GLL-19, were cloned in methylcellulose cultures and 20 randomly selected colonies from each line were picked up by micromanipulation, expanded in liquid cultures, and considered as clones of the original cell lines. The antigenic cell surface phenotype of these clones defined by panel of 12 monoclonal antibodies (MAb) was analyzed by flow microfluorometry (FMF) using a fluorescence-activated cell sorter (FACS II) and compared with the known stable phenotype of the parent cell line. The antibody panel consisted of eight MAb against melanoma-associated antigens, two MAb against monomorphic determinants of HLA-DR (la) and HLA-ABC, respectively, one MAb against the common acute lymphoblastic leukemia antigen (CALLA) and one MAb against carcinoembryonic antigen used as control. A remarkable heterogeneity in terms of qualitative and quantitative expression of the cell surface antigens studied was observed among and within the different clones. The single-cell origin of the clones was assessed by comparing the clonogenic cell frequency, determined by limiting dilutions in microculture plates, with the cloning efficiency observed in Petri dishes. Both techniques using methylcellulose medium gave the same percentages of growing colonies. Cells from four Me43 clones were recloned in methylcellulose and the phenotype of five randomly selected subclones from each clone was analysed using the same panel of monoclonal antibodies. Each subclone also displayed heterogeneity with individual phenotypes different from that of the original clone and from the parental Me43 cell line. The antigen expression by individual cells in situ within clones was analyzed on frozen sections from colonies using the same panel of MAb and a biotin-avidin immunoperoxidase method. The results confirmed the marked heterogeneity of antigen expression within and among colonies, as indicated by the FMF analysis.