Convergence in C([0,T];L2(Ω))C(\lbrack0,T\rbrack;L^2(\Omega)) of weak solutions to perturbed doubly degenerate parabolic equations

We study the behaviour of solutions to a class of nonlinear degenerate parabolic problems when the data are perturbed. The class includes the Richards equation, Stefan problem and the parabolic $p$-Laplace equation. We show that, up to a subsequence, weak solutions of the perturbed problem converge uniformly-in-time to weak solutions of the original problem as the perturbed data approach the original data. We do not assume uniqueness or additional regularity of the solution. However, when uniqueness is known, our result demonstrates that the weak solution is uniformly temporally stable to perturbations of the data. Beginning with a proof of temporally-uniform, spatially-weak convergence, we strengthen the latter by relating the unknown to an underlying convex structure that emerges naturally from energy estimates on the solution. The double degeneracy --- shown to be equivalent to a maximal monotone operator framework --- is handled with techniques inspired by a classical monotonicity argument and a simple variant of the compensated compactness phenomenon.Comment: J. Differential Equations, 201

Unified convergence analysis of numerical schemes for a miscible displacement problem

This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in $L^\infty(0,T; L^2(\Omega))$ of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes are compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion

How much larger quantum correlations are than classical ones

Considering as distance between two two-party correlations the minimum number of half local results one party must toggle in order to turn one correlation into the other, we show that the volume of the set of physically obtainable correlations in the Einstein-Podolsky-Rosen-Bell scenario is (3 pi/8)^2 = 1.388 larger than the volume of the set of correlations obtainable in local deterministic or probabilistic hidden-variable theories, but is only 3 pi^2/32 = 0.925 of the volume allowed by arbitrary causal (i.e., no-signaling) theories.Comment: REVTeX4, 6 page

Finite volume approximation of a diffusion-dissolution model and application to nuclear waste storage

International audienceThe study of two phase flow in porous media under high capillary pressures, in the case where one phase is incompressible and the other phase is gaseous, shows complex phenomena. We present in this paper a numerical approximation method, based on a two pressures formulation in the case where both phases are miscible, which is shown to also handle the limit case of immiscible phases. The space discretization is performed using a finite volume method, which can handle general grids. The efficiency of the formulation is shown on three numerical examples related tounderground waste disposal situations

Probabilistic analysis of the upwind scheme for transport

We provide a probabilistic analysis of the upwind scheme for multi-dimensional transport equations. We associate a Markov chain with the numerical scheme and then obtain a backward representation formula of Kolmogorov type for the numerical solution. We then understand that the error induced by the scheme is governed by the fluctuations of the Markov chain around the characteristics of the flow. We show, in various situations, that the fluctuations are of diffusive type. As a by-product, we prove that the scheme is of order 1/2 for an initial datum in BV and of order 1/2-a, for all a>0, for a Lipschitz continuous initial datum. Our analysis provides a new interpretation of the numerical diffusion phenomenon

Aircraft-based observations of air-sea fluxes over Denmark Strait and the Irminger sea during high wind speed conditions

The impact of targeted sonde observations on the 1-3 day forecasts for northern Europe is evaluated using the Met Office four-dimensional variational data assimilation scheme and a 24 km gridlength limited-area version of the Unified Model (MetUM). The targeted observations were carried out during February and March 2007 as part of the Greenland Flow Distortion Experiment, using a research aircraft based in Iceland. Sensitive area predictions using either total energy singular vectors or an ensemble transform Kalman filter were used to predict where additional observations should be made to reduce errors in the initial conditions of forecasts for northern Europe. Targeted sonde data was assimilated operationally into the MetUM. Hindcasts show that the impact of the sondes was mixed. Only two out of the five cases showed clear forecast improvement; the maximum forecast improvement seen over the verifying region was approximately 5% of the forecast error 24 hours into the forecast. These two cases are presented in more detail: in the first the improvement propagates into the verification region with a developing polar low; and in the second the improvement is associated with an upper-level trough. The impact of cycling targeted data in the background of the forecast (including the memory of previous targeted observations) is investigated. This is shown to cause a greater forecast impact, but does not necessarily lead to a greater forecast improvement. Finally, the robustness of the results is assessed using a small ensemble of forecasts

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.Comment: Computational Geosciences (2009) http://www.springerlink.com/content/933p55085742m203/?p=db14bb8c399b49979ba8389a3cae1b0f&pi=1

Finite volume scheme based on cell-vertex reconstructions for anisotropic diffusion problems with discontinuous coefficients

We propose a new second-order finite volume scheme for non-homogeneous and anisotropic diffusion problems based on cell to vertex reconstructions involving minimization of functionals to provide the coefficients of the cell to vertex mapping. The method handles complex situations such as large preconditioning number diffusion matrices and very distorted meshes. Numerical examples are provided to show the effectiveness of the method

Mutation update and genotype-phenotype correlations of novel and previously described mutations in TPM2 and TPM3 causing congenital myopathies

Mutations affecting skeletal muscle isoforms of the tropomyosin genes may cause nemaline myopathy, cap myopathy, core-rod myopathy, congenital fiber-type disproportion, distal arthrogryposes, and Escobar syndrome. We correlate the clinical picture of these diseases with novel (19) and previously reported (31) mutations of the TPM2 and TPM3 genes. Included are altogether 93 families: 53 with TPM2 mutations and 40 with TPM3 mutations. Thirty distinct pathogenic variants of TPM2 and 20 of TPM3 have been published or listed in the Leiden Open Variant Database (http://www.dmd.nl/). Most are heterozygous changes associated with autosomal-dominant disease. Patients with TPM2 mutations tended to present with milder symptoms than those with TPM3 mutations, DA being present only in the TPM2 group. Previous studies have shown that five of the mutations in TPM2 and one in TPM3 cause increased Ca2+ sensitivity resulting in a hypercontractile molecular phenotype. Patients with hypercontractile phenotype more often had contractures of the limb joints (18/19) and jaw (6/19) than those with nonhypercontractile ones (2/22 and 1/22), whereas patients with the non-hypercontractile molecular phenotype more often (19/22) had axial contractures than the hypercontractile group (7/19). Our in silico predictions show that most mutations affect tropomyosin–actin association or tropomyosin head-to-tail binding