Stability of preconditioned finite volume schemes at low Mach numbers

In [4], Guillard and Viozat propose a finite volume method for the simulation of inviscid steady as well as unsteady flows at low Mach numbers, based on a preconditioning technique. The scheme satisfies the results of a single scale asymptotic analysis in a discrete sense and comprises the advantage that this can be derived by a slight modification of the dissipation term within the numerical flux function. Unfortunately, it can be observed by numerical experiments that the preconditioned approach combined with an explicit time integration scheme turns out to be unstable if the time step Dt does not satisfy the requirement to be O(M2) as the Mach number M tends to zero, whereas the corresponding standard method remains stable up to Dt=O(M), M to 0, which results from the well-known CFL-condition. We present a comprehensive mathematical substantiation of this numerical phenomenon by means of a von Neumann stability analysis, which reveals that in contrast to the standard approach, the dissipation matrix of the preconditioned numerical flux function possesses an eigenvalue growing like M-2 as M tends to zero, thus causing the diminishment of the stability region of the explicit scheme. Thereby, we present statements for both the standard preconditioner used by Guillard and Viozat [4] and the more general one due to Turkel [21]. The theoretical results are after wards confirmed by numerical experiments.

The effect of vessel diameter on time dependent gas hold-up variations in highly viscous impeller agitated liquids

Time dependent gas hold-up generated in the 0.3 and 0.6 m diameter vessels using high viscosity castor oil and carboxy methyl cellulose (CMC) solution was compared on the basis of impeller speed (N) and gas velocity (V-G). Two types of hold-up were distinguished-the hold-up due to tiny bubbles (epsilon(ft)) and total hold-up (epsilon(f)), which included large and tiny bubbles. It was noted that vessel diameter (i.e. the scale of operation) significantly influences (i) the trends and the values of epsilon(f) and epsilon(ft), and (ii) the values of tau (a constant reflecting the time dependency of hold-up). The results showed that a scale independent correlation for gas hold-up of the form epsilon(f) or epsilon(ft) = A(N or P-G/V)(a) (V-G)(b), where "a" and "b" are positive constants is not appropriate for viscous liquids. This warrants further investigations into the effect of vessel diameter on gas hold-up in impeller agitated high viscosity liquids (mu or mu(a) > 0.4 Pa s). (C) 2003 Elsevier B.V. All rights reserved.

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

Applicability of the linear delta expansion for the lambda phi(4) field theory at finite temperature in the symmetric and broken phases

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

HTLV-1/2 seroprevalence and coinfection rate in Brazilian first-time blood donors: an 11-year follow-up

The seroprevalence and geographic distribution of HTLV-1/2 among blood donors are extremely important to transfusion services. We evaluated the seroprevalence of HTLV-1/2 infection among first-time blood donor candidates in Ribeirão Preto city and region. From January 2000 to December 2010, 1,038,489 blood donations were obtained and 301,470 were first-time blood donations. All samples were screened with serological tests for HTLV-1/2 using enzyme immunoassay (EIA). In addition, the frequency of coinfection with hepatitis B virus (HBV), hepatitis C virus (HCV), human immunodeficiency virus (HIV), Chagas disease (CD) and syphilis was also determined. In-house PCR was used as confirmatory test for HTLV-1/2. A total of 296 (0.1%) first-time donors were serologically reactive for HTLV-1/2. Confirmatory PCR of 63 samples showed that 28 were HTLV-1 positive, 13 HTLV-2 positive, 19 negative and three indeterminate. Regarding HTLV coinfection rates, the most prevalent was with HBV (51.3%) and HCV (35.9%), but coinfection with HIV, CD and syphilis was also detected. The real number of HTLV-infected individual and coinfection rate in the population is underestimated and epidemiological studies like ours are very informative.

CT screening and follow-up of lung nodules: effects of tube current-time setting and nodule size and density on detectability and of tube current-time setting on apparent size

The purpose of the study was to quantify and compare the effect of CT dose and of size and density of nodules on the detectability of lung nodules and to quantify the influence of CT dose on the size of the nodules.

Time-lapse microscopy and classification of 2D human mesenchymal stem cells based on cell shape picks up myogenic from osteogenic and adipogenic differentiation

Current methods to characterize mesenchymal stem cells (MSCs) are limited to CD marker expression, plastic adherence and their ability to differentiate into adipogenic, osteogenic and chondrogenic precursors. It seems evident that stem cells undergoing differentiation should differ in many aspects, such as morphology and possibly also behaviour; however, such a correlation has not yet been exploited for fate prediction of MSCs. Primary human MSCs from bone marrow were expanded and pelleted to form high-density cultures and were then randomly divided into four groups to differentiate into adipogenic, osteogenic chondrogenic and myogenic progenitor cells. The cells were expanded as heterogeneous and tracked with time-lapse microscopy to record cell shape, using phase-contrast microscopy. The cells were segmented using a custom-made image-processing pipeline. Seven morphological features were extracted for each of the segmented cells. Statistical analysis was performed on the seven-dimensional feature vectors, using a tree-like classification method. Differentiation of cells was monitored with key marker genes and histology. Cells in differentiation media were expressing the key genes for each of the three pathways after 21 days, i.e. adipogenic, osteogenic and chondrogenic, which was also confirmed by histological staining. Time-lapse microscopy data were obtained and contained new evidence that two cell shape features, eccentricity and filopodia (= 'fingers') are highly informative to classify myogenic differentiation from all others. However, no robust classifiers could be identified for the other cell differentiation paths. The results suggest that non-invasive automated time-lapse microscopy could potentially be used to predict the stem cell fate of hMSCs for clinical application, based on morphology for earlier time-points. The classification is challenged by cell density, proliferation and possible unknown donor-specific factors, which affect the performance of morphology-based approaches. Copyright © 2012 John Wiley & Sons, Ltd.

Time-related, cross-sectional and prospective follow-up of pancreatic endocrine function after pancreas allograft transplantation in type 1 (insulin-dependent) diabetic patients

It has been established that successful pancreas transplantation in Type 1 (insulin-dependent) diabetic patients results in normal but exaggerated phasic glucose-induced insulin secretion, normal intravenous glucose disappearance rates, improved glucose recovery from insulin-induced hypoglycaemia, improved glucagon secretion during insulin-induced hypoglycaemia, but no alterations in pancreatic polypeptide responses to hypoglycaemia. However, previous reports have not segregated the data in terms of the length of time following successful transplantation and very little prospective data collected over time in individual patients has been published. This article reports that in general there are no significant differences in the level of improvement when comparing responses as early as three months post-operatively up to as long as two years post-operatively when examining the data cross-sectionally in patients who have successfully maintained their allografts. Moreover, this remarkable constancy in pancreatic islet function is also seen in a smaller group of patients who have been examined prospectively at various intervals post-operatively. It is concluded that successful pancreas transplantation results in remarkable improvements in Alpha and Beta cell but not PP cell function that are maintained for at least one to two years.