A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Tax simulation in the SHP

This paper describes a simulation package designed to estimate the annual income taxes paid by respondents of the Swiss Household Panel (SHP). In Switzerland, the 26 cantons have their own tax system. Additionally, tax levels vary between the over 2000 municipalities and over time. The simulation package takes account of this complexity by building on existing tables on tax levels which are provided by the Swiss Federal Tax Administration Office. Because these are limited to a few types of households and only 812 municipalities, they have to be extended to cover all households and municipalities. A further drawback of these tables is that they neglect several deductions. The tax simulation package fills this gap by taking additionally account of deductions for children, double-earner couples, third pillar and support for dependent persons according to cantonal legislation. The resulting variable on direct taxes not only serves to calculate household income net of taxes, but can also be a variable for analysis by its own account.

Phase 3 randomized trial on larynx preservation comparing sequential vs alternating chemotherapy and radiotherapy.

BACKGROUND: Both induction chemotherapy followed by irradiation and concurrent chemotherapy and radiotherapy have been reported as valuable alternatives to total laryngectomy in patients with advanced larynx or hypopharynx cancer. We report results of the randomized phase 3 trial 24954 from the European Organization for Research and Treatment of Cancer. METHODS: Patients with resectable advanced squamous cell carcinoma of the larynx (tumor stage T3-T4) or hypopharynx (T2-T4), with regional lymph nodes in the neck staged as N0-N2 and with no metastasis, were randomly assigned to treatment in the sequential (or control) or the alternating (or experimental) arm. In the sequential arm, patients with a 50% or more reduction in primary tumor size after two cycles of cisplatin and 5-fluorouracil received another two cycles, followed by radiotherapy (70 Gy total). In the alternating arm, a total of four cycles of cisplatin and 5-fluorouracil (in weeks 1, 4, 7, and 10) were alternated with radiotherapy with 20 Gy during the three 2-week intervals between chemotherapy cycles (60 Gy total). All nonresponders underwent salvage surgery and postoperative radiotherapy. The Kaplan-Meier method was used to obtain time-to-event data. RESULTS: The 450 patients were randomly assigned to treatment (224 to the sequential arm and 226 to the alternating arm). Median follow-up was 6.5 years. Survival with a functional larynx was similar in sequential and alternating arms (hazard ratio of death and/or event = 0.85, 95% confidence interval = 0.68 to 1.06), as were median overall survival (4.4 and 5.1 years, respectively) and median progression-free interval (3.0 and 3.1 years, respectively). Grade 3 or 4 mucositis occurred in 64 (32%) of the 200 patients in the sequential arm who received radiotherapy and in 47 (21%) of the 220 patients in the alternating arm. Late severe edema and/or fibrosis was observed in 32 (16%) patients in the sequential arm and in 25 (11%) in the alternating arm. CONCLUSIONS: Larynx preservation, progression-free interval, and overall survival were similar in both arms, as were acute and late toxic effects.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

The Sequential Organ Failure Assessment score and copeptin for predicting survival in ventilator-associated pneumonia.

INTRODUCTION: Ventilator-associated pneumonia remains the most common nosocomial infection in the critically ill and contributes to significant morbidity. Eventual decisions regarding withdrawal or maximal therapy are demanding and rely on physicians' experience. Additional objective tools for risk assessment may improve medical judgement. Copeptin, reflecting vasopressin release, as well as the Sequential Organ Failure Assessment (SOFA) score, reflecting the individual degree of organ dysfunction, might qualify for survival prediction in ventilator-associated pneumonia. We investigated the predictive value of the SOFA score and copeptin in ventilator-associated pneumonia. METHODS: One hundred one patients with ventilator-associated pneumonia were prospectively assessed. Death within 28 days after ventilator-associated pneumonia onset was the primary end point. RESULTS: The SOFA score and the copeptin levels at ventilator-associated pneumonia onset were significantly elevated in nonsurvivors (P = .002 and P = .017, respectively). Both markers had different time courses in survivors and nonsurvivors (P < .001 and P = .006). Mean SOFA (average SOFA of 10 days after VAP onset) was superior in predicting 28-day survival as compared with SOFA and copeptin at ventilator-associated pneumonia onset (area under the curve, 0.90 vs 0.73 and 0.67, respectively). CONCLUSIONS: The predictive value of serial-measured SOFA significantly exceeds those of single SOFA and copeptin measurements. Serial SOFA scores accurately predict outcome in ventilator-associated pneumonia.