Analysis of a finite volume element method for the Stokes problem

In this paper we propose a stabilized conforming finite volume element method for the Stokes equations. On stating the convergence of the method, optimal a priori error estimates in different norms are obtained by establishing the adequate connection between the finite volume and stabilized finite element formulations. A superconvergence result is also derived by using a postprocessing projection method. In particular, the stabilization of the continuous lowest equal order pair finite volume element discretization is achieved by enriching the velocity space with local functions that do not necessarily vanish on the element boundaries. Finally, some numerical experiments that confirm the predicted behavior of the method are provided.

Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

Fixed rings and strictures in Eosinophilic Esophagitis develop due to continuing inflammation over time

Background and Aims: Two distinct e ndoscopic phenotypes of E osinophilic Esophagitis (EoE) h ave been identified: t he inflammatory (IP) a nd the stenosing (SP) p henotype. I t is not known whether these EoE-associated phenotypes are reflective of different phases during disease course. We aimed to assess the phenotype a t initial EoE p resentation and d iagnosis and to evaluate if SP increases over time. Methods: R etrospective a nalysis of t he Swiss EoE Database (SEED) extended b y a review of p atients charts, endoscopy and pathology records. Results: F orty-four E oE p atients were a nalyzed (33 males, mean age at index visit 41 ± 14 years, all Caucasians). Median follow-up t ime was 3.1 years (IQR 1-4, r ange 1 -18 years). Median diagnostic delay w as 5 y ears (IQR 2-16, range 0-34 years). A t first diagnosis, 3 2% ( 14/44) o f EoE patients h ad already presented w ith a stenosis. T he mean d iameter o f the stenoses w as 1 0 ± 2 mm, and the mean length was 2 .8 ± 2 .9 cm. Peak e osinophil count d id n ot c hange over t ime (48 ± 39 eos/HPF at index visit vs. 59 ± 41 eos/HPF at end of follow-up, n=44). The risk of the presence of a stenosis at index visit was 0% f or a d isease duration of 0 -4 y ears, 37% f or a d isease duration between 5-10 years and 67% f or a d isease duration >10 years (p = 0.0035, trend test). Conclusions: T he frequency of e sophageal stenoses i s proportional to the disease duration, whereas the inflammatory activity does n ot s ignificantly c hange over t ime. O ur f indings underscore the necessity to reduce diagnostic delay in EoE and to control the underlying inflammatory processes to prevent esophageal remodeling.

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.

Random mating with a finite number of matings.

Random mating is the null model central to population genetics. One assumption behind random mating is that individuals mate an infinite number of times. This is obviously unrealistic. Here we show that when each female mates a finite number of times, the effective size of the population is substantially decreased.

RadA protein from Archaeoglobus fulgidus forms rings, nucleoprotein filaments and catalyses homologous recombination.

Proteins that catalyse homologous recombination have been identified in all living organisms and are essential for the repair of damaged DNA as well as for the generation of genetic diversity. In bacteria homologous recombination is performed by the RecA protein, whereas in the eukarya a related protein called Rad51 is required to catalyse recombination and repair. More recently, archaeal homologues of RecA/Rad51 (RadA) have been identified and isolated. In this work we have cloned and purified the RadA protein from the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus and characterised its in vitro activities. We show that (i) RadA protein forms ring structures in solution and binds single- but not double-stranded DNA to form nucleoprotein filaments, (ii) RadA is a single-stranded DNA-dependent ATPase at elevated temperatures, and (iii) RadA catalyses efficient D-loop formation and strand exchange at temperatures of 60-70 degrees C. Finally, we have used electron microscopy to visualise RadA-mediated joint molecules, the intermediates of homologous recombination. Intriguingly, RadA shares properties of both the bacterial RecA and eukaryotic Rad51 recombinases.

Development and experimental validation of a finite element model of total ankle replacement.

Total ankle replacement remains a less satisfactory solution compared to other joint replacements. The goal of this study was to develop and validate a finite element model of total ankle replacement, for future testing of hypotheses related to clinical issues. To validate the finite element model, an experimental setup was specifically developed and applied on 8 cadaveric tibias. A non-cemented press fit tibial component of a mobile bearing prosthesis was inserted into the tibias. Two extreme anterior and posterior positions of the mobile bearing insert were considered, as well as a centered one. An axial force of 2kN was applied for each insert position. Strains were measured on the bone surface using digital image correlation. Tibias were CT scanned before implantation, after implantation, and after mechanical tests and removal of the prosthesis. The finite element model replicated the experimental setup. The first CT was used to build the geometry and evaluate the mechanical properties of the tibias. The second CT was used to set the implant position. The third CT was used to assess the bone-implant interface conditions. The coefficient of determination (R-squared) between the measured and predicted strains was 0.91. Predicted bone strains were maximal around the implant keel, especially at the anterior and posterior ends. The finite element model presented here is validated for future tests using more physiological loading conditions.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).