Monolithic multigrid method for the coupled stokes flow and deformable porous medium system

The interaction between fluid flow and a deformable porous medium is a complicated multi-physics problem, which can be described by a coupled model based on the Stokes and poroelastic equations. A monolithic multigrid method together with either a coupled Vanka smoother or a decoupled Uzawa smoother is employed as an efficient numerical technique for the linear discrete system obtained by finite volumes on staggered grids. A specialty in our modeling approach is that at the interface of the fluid and poroelastic medium, two unknowns from the different subsystems are defined at the same grid point. We propose a special discretization at and near the points on the interface, which combines the approximation of the governing equations and the considered interface conditions. In the decoupled Uzawa smoother, Local Fourier Analysis (LFA) helps us to select optimal values of the relaxation parameter appearing. To implement the monolithic multigrid method, grid partitioning is used to deal with the interface updates when communication is required between two subdomains. Numerical experiments show that the proposed numerical method has an excellent convergence rate. The efficiency and robustness of the method are confirmed in numerical experiments with typically small realistic values of the physical coefficients

On a multigrid method for the coupled Stokes and porous media flow problem

The multigrid solution of coupled porous media and Stokes flow problems is considered. The Darcy equation as the saturated porous medium model is coupled to the Stokes equations by means of appropriate interface conditions. We focus on an efficient multigrid solution technique for the coupled problem, which is discretized by finite volumes on staggered grids, giving rise to a saddle point linear system. Special treatment is required regarding the discretization at the interface. An Uzawa smoother is employed in multigrid, which is a decoupled procedure based on symmetric Gauss-Seidel smoothing for velocity components and a simple Richardson iteration for the pressure field. Since a relaxation parameter is part of a Richardson iteration, Local Fourier Analysis (LFA) is applied to determine the optimal parameters. Highly satisfactory multigrid convergence is reported, and, moreover, the algorithm performs well for small values of the hydraulic conductivity and fluid viscosity, that are relevant for applications

On an Uzawa smoother in multigrid for poroelasticity equations

A poroelastic saturated medium can be modeled by means of Biot''s theory of consolidation. It describes the time-dependent interaction between the deformation of porous material and the fluid flow inside of it. Here, for the efficient solution of the poroelastic equations, a multigrid method is employed with an Uzawa-type iteration as the smoother. The Uzawa smoother is an equation-wise procedure. It shall be interpreted as a combination of the symmetric Gauss-Seidel smoothing for displacements, together with a Richardson iteration for the Schur complement in the pressure field. The Richardson iteration involves a relaxation parameter which affects the convergence speed, and has to be carefully determined. The analysis of the smoother is based on the framework of local Fourier analysis (LFA) and it allows us to provide an analytic bound of the smoothing factor of the Uzawa smoother as well as an optimal value of the relaxation parameter. Numerical experiments show that our upper bound provides a satisfactory estimate of the exact smoothing factor, and the selected relaxation parameter is optimal. In order to improve the convergence performance, the acceleration of multigrid by iterant recombination is taken into account. Numerical results confirm the efficiency and robustness of the acceleration scheme

B→DsπB \to D_s \pi and the tree amplitude in B→π+π−B \to \pi^+ \pi^-

The recently-observed decay $B^0 \to D_s^+ \pi^-$ is expected to proceed mainly by means of a tree amplitude in the factorization limit: $B^0 \to \pi^- {(W^+)}^*$, ${(W^+)}^* \to D_s^+$. Under this assumption, we predict the corresponding contribution of the tree amplitude to $B^0 \to \pi^+ \pi^-$. We indicate the needed improvements in data that will allow a useful estimate of this amplitude with errors comparable to those accompanying other methods. Since the factorization hypothesis for this process goes beyond that proved in most approaches, we also discuss independent tests of this hypothesis.Comment: 7 pages, LaTeX, 1 figure, to be submitted to Phys. Rev. D (Brief Reports

Multigrid method for nonlinear poroelasticity equations

In this study, a nonlinear multigrid method is applied for solving the system of incompressible poroelasticity equations considering nonlinear hydraulic conductivity. For the unsteady problem, an additional artificial term is utilized to stabilize the solutions when the equations are discretized on collocated grids. We employ two nonlinear multigrid methods, i.e. the “full approximation scheme” and “Newton multigrid” for solving the corresponding system of equations arising after discretization. For the steady case, both homogeneous and heterogeneous cases are solved and two different smoothers are examined to search for an efficient multigrid method. Numerical results show a good convergence performance for all the strategies

Uzawa smoother in multigrid for the coupleD porous medium and stokes flow system

The multigrid solution of coupled porous media and Stokes flow problems is considered. The Darcy equation as the saturated porous medium model is coupled to the Stokes equations by means of appropriate interface conditions. We focus on an efficient multigrid solution technique for the coupled problem, which is discretized by finite volumes on staggered grids, giving rise to a saddle point linear system. Special treatment is required regarding the discretization at the interface. An Uzawa smoother is employed in multigrid, which is a decoupled procedure based on symmetric Gauss–Seidel smoothing for velocity components and a simple Richardson iteration for the pressure field. Since a relaxation parameter is part of a Richardson iteration, local Fourier analysis is applied to determine the optimal parameters. Highly satisfactory multigrid convergence is reported, and, moreover, the algorithm performs very well for small values of the hydraulic conductivity and fluid viscosity, which are relevant for applications

A multigrid multilevel Monte Carlo method for transport in the Darcy–Stokes system

A multilevel Monte Carlo (MLMC) method for Uncertainty Quantification (UQ) of advection-dominated contaminant transport in a coupled Darcy–Stokes flow system is described. In particular, we focus on high-dimensional epistemic uncertainty due to an unknown permeability field in the Darcy domain that is modelled as a lognormal random field. This paper explores different numerical strategies for the subproblems and suggests an optimal combination for the MLMC estimator. We propose a specific monolithic multigrid algorithm to efficiently solve the steady-state Darcy–Stokes flow with a highly heterogeneous diffusion coefficient. Furthermore, we describe an Alternating Direction Implicit (ADI) based time-stepping for the flux-limited quadratic upwinding discretization for the transport problem. Numerical experiments illustrating the multigrid convergence and cost of the MLMC estimator with respect to the smoothness of permeability field are presented

Crystallization and preliminary X-ray diffraction analysis of PAT, an acetyltransferase from Sulfolobus solfataricus

PAT is an acetyltransferase from the archaeon Sulfolobus solfataricus that specifically acetylates the chromatin protein Alba. The enzyme was expressed, purified and subsequently crystallized using the sitting-drop vapour-diffusion technique. Native diffraction data were collected to 1.70 angstrom resolution on the BL13C1 beamline of NSRRC from a flash-frozen crystal at 100 K. The crystals belonged to space group P2(1)2(1)2(1), with unit-cell parameters a = 44.30, b = 46.59, c = 68.39 angstrom

Direct CP Violation Asymmetries in Exclusive B Decays in a Bethe-Salpeter Approach

CP asymmetry in some exclusive decays of charged B meson are calculated in a Bethe-Salpeter approach. Hadronic final state interactions are ignored. Complex decay amplitudes are assumed to arise entirely from perturbative quark-antiquark loops. Calculations are done both with and without the gluon quark-antiquark vacuum polarization loops. The effects of neglecting the imaginary parts arising from the diagonal quark-antiquark loops are also studied.Comment: 15 pages, Latex, 2 eps-figures. Minor revisions. To be published in Phys.Lett.

A Lanczos algorithm for linear response

An iterative algorithm is presented for solving the RPA equations of linear response. The method optimally computes the energy-weighted moments of the strength function, allowing one to match the computational effort to the intrinsic accuracy of the basic mean-field approximation, avoiding the problem of solving very large matrices. For local interactions, the computational effort for the method scales with the number of particles N_p as O(N_p^3).Comment: 12 pages including 3 figures; Late