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.

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

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.

Apports de la médecine nucléaire en cardiologie préventive: l'exemple du syndrome métabolique [Value of nuclear medicine in preventive cardiology: the metabolic syndrome example]

Metabolic syndrome represents a grouping of risk factors closely linked to cardiovascular diseases and diabetes. At first, nuclear medicine has no direct application in cardiology at the level of primary prevention, but positron emission tomography is a non invasive imaging technique that can assess myocardial perfusion as well as the endothelium-dependent coronary vasomotion--a surrogate marker of cardiovascular event rate--thus finding an application in studying coronary physiopathology. As the prevalence of the metabolic syndrome is still unknown in Switzerland, we will estimate it from data available in the frame of a health promotion program. Based on the deleterious effect on the endothelium already observed with two components, we will estimate the number of persons at risk in Switzerland.

Changes in volume, surface estimate, three-dimensional shape and total number of neurons of the human primary visual cortex from midgestation until old age.

Macroscopic features such as volume, surface estimate, thickness and caudorostral length of the human primary visual cortex (Brodman's area 17) of 46 human brains between midgestation and 93 years were studied by means of camera lucida drawings from serial frontal sections. Individual values were best fitted by a logistic function from midgestation to adulthood and by a regression line between adulthood and old age. Allometric functions were calculated to study developmental relationships between all the features. The three-dimensional shape of area 17 was also reconstructed from the serial sections in 15 cases and correlated with the sequence of morphological events. The sulcal pattern of area 17 begins to develop around 21 weeks of gestation but remains rather simple until birth, while it becomes more convoluted, particularly in the caudal part, during the postnatal period. Until birth, a large increase in cortical thickness (about 83% of its mean adult value) and caudorostral length (69%) produces a moderate increase in cortical volume (31%) and surface estimate (40%) of area 17. After birth, the cortical volume and surface undergo their maximum growth rate, in spite of a rather small increase in cortical thickness and caudorostral length. This is due to the development of the pattern of gyrification within and around the calcarine fissure. All macroscopic features have reached the mean adult value by the end of the first postnatal year. With aging, the only features to undergo significant regression are the cortical surface estimate and the caudorostral length. The total number of neurons in area 17 shows great interindividual variability at all ages. No decrease in the postnatal period or in aging could be demonstrated.

A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.

The risk factors for fractures and trabecular bone-score value in patients with endogenous Cushing's syndrome.

In a cohort study of 182 consecutive patients with active endogenous Cushing's syndrome, the only predictor of fracture occurrence after adjustment for age, gender bone mineral density (BMD) and trabecular bone score (TBS) was 24-h urinary free cortisol (24hUFC) levels with a threshold of 1472 nmol/24 h (odds ratio, 3.00 (95 % confidence interval (CI), 1.52-5.92); p = 0.002). INTRODUCTION: The aim was to estimate the risk factors for fracture in subjects with endogenous Cushing's syndrome (CS) and to evaluate the value of the TBS in these patients. METHODS: All enrolled patients with CS (n = 182) were interviewed in relation to low-traumatic fractures and underwent lateral X-ray imaging from T4 to L5. BMD measurements were performed using a DXA Prodigy device (GEHC Lunar, Madison, Wisconsin, USA). The TBS was derived retrospectively from existing BMD scans, blinded to clinical outcome, using TBS iNsight software v2.1 (Medimaps, Merignac, France). Urinary free cortisol (24hUFC) was measured by immunochemiluminescence assay (reference range, 60-413 nmol/24 h). RESULTS: Among enrolled patients with CS (149 females; 33 males; mean age, 37.8 years (95 % confidence interval, 34.2-39.1); 24hUFC, 2370 nmol/24 h (2087-2632), fractures were confirmed in 81 (44.5 %) patients, with 70 suffering from vertebral fractures, which were multiple in 53 cases; 24 patients reported non-vertebral fractures. The mean spine TBS was 1.207 (1.187-1.228), and TBS Z-score was -1.86 (-2.07 to -1.65); area under the curve (AUC) was used to predict fracture (mean spine TBS) = 0.548 (95 % CI, 0.454-0.641)). In the final regression model, the only predictor of fracture occurrence was 24hUFC levels (p = 0.001), with an increase of 1.041 (95 % CI, 1.019-1.063), calculated for every 100 nmol/24-h cortisol elevation (AUC (24hUFC) = 0.705 (95 % CI, 0.629-0.782)). CONCLUSIONS: Young patients with CS have a low TBS. However, the only predictor of low traumatic fracture is the severity of the disease itself, indicated by high 24hUFC levels.

New measurements of energy expenditure and physical activity in chronic kidney disease.

The accurate estimation of total daily energy expenditure (TEE) in chronic kidney patients is essential to allow the provision of nutritional requirements; however, it remains a challenge to collect actual physical activity and resting energy expenditure in maintenance dialysis patients. The direct measurement of TEE by direct calorimetry or doubly labeled water cannot be used easily so that, in clinical practice, TEE is usually estimated from resting energy expenditure and physical activity. Prediction equations may also be used to estimate resting energy expenditure; however, their use has been poorly documented in dialysis patients. Recently, a new system called SenseWear Armband (BodyMedia, Pittsburgh, PA) was developed to assess TEE, but so far no data have been published in chronic kidney disease patients. The aim of this review is to describe new measurements of energy expenditure and physical activity in chronic kidney disease patients.

Screening primary-care patients forgoing health care for economic reasons.

BACKGROUND: Growing social inequities have made it important for general practitioners to verify if patients can afford treatment and procedures. Incorporating social conditions into clinical decision-making allows general practitioners to address mismatches between patients' health-care needs and financial resources. OBJECTIVES: Identify a screening question to, indirectly, rule out patients' social risk of forgoing health care for economic reasons, and estimate prevalence of forgoing health care and the influence of physicians' attitudes toward deprivation. DESIGN: Multicenter cross-sectional survey. PARTICIPANTS: Forty-seven general practitioners working in the French-speaking part of Switzerland enrolled a random sample of patients attending their private practices. MAIN MEASURES: Patients who had forgone health care were defined as those reporting a household member (including themselves) having forgone treatment for economic reasons during the previous 12 months, through a self-administered questionnaire. Patients were also asked about education and income levels, self-perceived social position, and deprivation levels. KEY RESULTS: Overall, 2,026 patients were included in the analysis; 10.7% (CI95% 9.4-12.1) reported a member of their household to have forgone health care during the 12 previous months. The question "Did you have difficulties paying your household bills during the last 12 months" performed better in identifying patients at risk of forgoing health care than a combination of four objective measures of socio-economic status (gender, age, education level, and income) (R(2) = 0.184 vs. 0.083). This question effectively ruled out that patients had forgone health care, with a negative predictive value of 96%. Furthermore, for physicians who felt powerless in the face of deprivation, we observed an increase in the odds of patients forgoing health care of 1.5 times. CONCLUSION: General practitioners should systematically evaluate the socio-economic status of their patients. Asking patients whether they experience any difficulties in paying their bills is an effective means of identifying patients who might forgo health care.

3-T direct MR arthrography of the wrist: value of finger trap distraction to assess intrinsic ligament and triangular fibrocartilage complex tears.

PURPOSE: To determine the value of applying finger trap distraction during direct MR arthrography of the wrist to assess intrinsic ligament and triangular fibrocartilage complex (TFCC) tears. MATERIALS AND METHODS: Twenty consecutive patients were prospectively investigated by three-compartment wrist MR arthrography. Imaging was performed with 3-T scanners using a three-dimensional isotropic (0.4 mm) T1-weighted gradient-recalled echo sequence, with and without finger trap distraction (4 kg). In a blind and independent fashion, two musculoskeletal radiologists measured the width of the scapholunate (SL), lunotriquetral (LT) and ulna-TFC (UTFC) joint spaces. They evaluated the amount of contrast medium within these spaces using a four-point scale, and assessed SL, LT and TFCC tears, as well as the disruption of Gilula's carpal arcs. RESULTS: With finger trap distraction, both readers found a significant increase in width of the SL space (mean Δ = +0.1mm, p ≤ 0.040), and noticed more contrast medium therein (p ≤ 0.035). In contrast, the differences in width of the LT (mean Δ = +0.1 mm, p ≥ 0.057) and UTFC (mean Δ = 0mm, p ≥ 0.728) spaces, as well as the amount of contrast material within these spaces were not statistically significant (p = 0.607 and ≥ 0.157, respectively). Both readers detected more SL (Δ = +1, p = 0.157) and LT (Δ = +2, p = 0.223) tears, although statistical significance was not reached, and Gilula's carpal arcs were more frequently disrupted during finger trap distraction (Δ = +5, p = 0.025). CONCLUSION: The application of finger trap distraction during direct wrist MR arthrography may enhance both detection and characterisation of SL and LT ligament tears by widening the SL space and increasing the amount of contrast within the SL and LT joint spaces.

Is the formula of Traub still up to date in antemortem blood glucose level estimation?

According to the hypothesis of Traub, also known as the 'formula of Traub', postmortem values of glucose and lactate found in the cerebrospinal fluid or vitreous humor are considered indicators of antemortem blood glucose levels. However, because the lactate concentration increases in the vitreous and cerebrospinal fluid after death, some authors postulated that using the sum value to estimate antemortem blood glucose levels could lead to an overestimation of the cases of glucose metabolic disorders with fatal outcomes, such as diabetic ketoacidosis. The aim of our study, performed on 470 consecutive forensic cases, was to ascertain the advantages of the sum value to estimate antemortem blood glucose concentrations and, consequently, to rule out fatal diabetic ketoacidosis as the cause of death. Other biochemical parameters, such as blood 3-beta-hydroxybutyrate, acetoacetate, acetone, glycated haemoglobin and urine glucose levels, were also determined. In addition, postmortem native CT scan, autopsy, histology, neuropathology and toxicology were performed to confirm diabetic ketoacidosis as the cause of death. According to our results, the sum value does not add any further information for the estimation of antemortem blood glucose concentration. The vitreous glucose concentration appears to be the most reliable marker to estimate antemortem hyperglycaemia and, along with the determination of other biochemical markers (such as blood acetone and 3-beta-hydroxybutyrate, urine glucose and glycated haemoglobin), to confirm diabetic ketoacidosis as the cause of death.

Impact of low-level viremia on clinical and virological outcomes in treated HIV-1-infected patients.

BACKGROUND: The goal of antiretroviral therapy (ART) is to reduce HIV-related morbidity and mortality by suppressing HIV replication. The prognostic value of persistent low-level viremia (LLV), particularly for clinical outcomes, is unknown. OBJECTIVE: Assess the association of different levels of LLV with virological failure, AIDS event, and death among HIV-infected patients receiving combination ART. METHODS: We analyzed data from 18 cohorts in Europe and North America, contributing to the ART Cohort Collaboration. Eligible patients achieved viral load below 50 copies/ml within 3-9 months after ART initiation. LLV50-199 was defined as two consecutive viral loads between 50 and 199 copies/ml and LLV200-499 as two consecutive viral loads between 50 and 499 copies/ml, with at least one between 200 and 499 copies/ml. We used Cox models to estimate the association of LLV with virological failure (two consecutive viral loads at least 500 copies/ml or one viral load at least 500 copies/ml, followed by a modification of ART) and AIDS event/death. RESULTS: Among 17 902 patients, 624 (3.5%) experienced LLV50-199 and 482 (2.7%) LLV200-499. Median follow-up was 2.3 and 3.1 years for virological and clinical outcomes, respectively. There were 1903 virological failure, 532 AIDS events and 480 deaths. LLV200-499 was strongly associated with virological failure [adjusted hazard ratio (aHR) 3.97, 95% confidence interval (CI) 3.05-5.17]. LLV50-199 was weakly associated with virological failure (aHR 1.38, 95% CI 0.96-2.00). LLV50-199 and LLV200-499 were not associated with AIDS event/death (aHR 1.19, 95% CI 0.78-1.82; and aHR 1.11, 95% CI 0.72-1.71, respectively). CONCLUSION: LLV200-499 was strongly associated with virological failure, but not with AIDS event/death. Our results support the US guidelines, which define virological failure as a confirmed viral load above 200 copies/ml.