87 research outputs found

    Hierarchical error estimates for the energy functional in obstacle problems

    Get PDF
    We present a hierarchical a posteriori error analysis for the minimum value of the energy functional in symmetric obstacle problems. The main result is that the error in the energy minimum is, up to oscillation terms, equivalent to an appropriate hierarchical estimator. The proof does not invoke any saturation assumption. We even show that small oscillation implies a related saturation assumption. In addition, we prove efficiency and reliability of an a posteriori estimate of the discretization error and thereby cast some light on the theoretical understanding of previous hierarchical estimators. Finally, we illustrate our theoretical results by numerical computations

    Computing the first eigenpair of the p-Laplacian via inverse iteration of sublinear supersolutions

    Full text link
    We introduce an iterative method for computing the first eigenpair (λp,ep)(\lambda_{p},e_{p}) for the pp-Laplacian operator with homogeneous Dirichlet data as the limit of (μq,uq)(\mu_{q,}u_{q}) as qpq\rightarrow p^{-}, where uqu_{q} is the positive solution of the sublinear Lane-Emden equation Δpuq=μquqq1-\Delta_{p}u_{q}=\mu_{q}u_{q}^{q-1} with same boundary data. The method is shown to work for any smooth, bounded domain. Solutions to the Lane-Emden problem are obtained through inverse iteration of a super-solution which is derived from the solution to the torsional creep problem. Convergence of uqu_{q} to epe_{p} is in the C1C^{1}-norm and the rate of convergence of μq\mu_{q} to λp\lambda_{p} is at least O(pq)O(p-q). Numerical evidence is presented.Comment: Section 5 was rewritten. Jed Brown was added as autho

    Spallation reactions. A successful interplay between modeling and applications

    Get PDF
    The spallation reactions are a type of nuclear reaction which occur in space by interaction of the cosmic rays with interstellar bodies. The first spallation reactions induced with an accelerator took place in 1947 at the Berkeley cyclotron (University of California) with 200 MeV deuterons and 400 MeV alpha beams. They highlighted the multiple emission of neutrons and charged particles and the production of a large number of residual nuclei far different from the target nuclei. The same year R. Serber describes the reaction in two steps: a first and fast one with high-energy particle emission leading to an excited remnant nucleus, and a second one, much slower, the de-excitation of the remnant. In 2010 IAEA organized a worskhop to present the results of the most widely used spallation codes within a benchmark of spallation models. If one of the goals was to understand the deficiencies, if any, in each code, one remarkable outcome points out the overall high-quality level of some models and so the great improvements achieved since Serber. Particle transport codes can then rely on such spallation models to treat the reactions between a light particle and an atomic nucleus with energies spanning from few tens of MeV up to some GeV. An overview of the spallation reactions modeling is presented in order to point out the incomparable contribution of models based on basic physics to numerous applications where such reactions occur. Validations or benchmarks, which are necessary steps in the improvement process, are also addressed, as well as the potential future domains of development. Spallation reactions modeling is a representative case of continuous studies aiming at understanding a reaction mechanism and which end up in a powerful tool.Comment: 59 pages, 54 figures, Revie

    A posteriori error estimates for the virtual element method

    Get PDF
    An a posteriori error analysis for the virtual element method (VEM) applied to general elliptic problems is presented. The resulting error estimator is of residual-type and applies on very general polygonal/polyhedral meshes. The estimator is fully computable as it relies only on quantities available from the VEM solution, namely its degrees of freedom and element-wise polynomial projection. Upper and lower bounds of the error estimator with respect to the VEM approximation error are proven. The error estimator is used to drive adaptive mesh refinement in a number of test problems. Mesh adaptation is particularly simple to implement since elements with consecutive co-planar edges/faces are allowed and, therefore, locally adapted meshes do not require any local mesh post-processing

    Convergence of Adaptive Finite Element Methods with Error-Dominated Oscillation

    No full text
    Recently, we devised an approach to a posteriori error analysis, which clarifies the role of oscillation and where oscillation is bounded in terms of the current approximation error. Basing upon this approach, we derive plain convergence of adaptive linear finite elements approximating the Poisson problem. The result covers arbritray H1H^{-1}-data and characterizes convergent marking strategies
    corecore