Hierarchical error estimates for the energy functional in obstacle problems

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

Quasi-optimal nonconforming methods for symmetric elliptic problems. I—abstract theory

We consider nonconforming methods for symmetric elliptic problems and characterize their quasi-optimality in terms of suitable notions of stability and consistency. The quasi-optimality constant is determined, and the possible impact of nonconformity on its size is quantified by means of two alternative consistency measures. Identifying the structure of quasi-optimal methods, we show that their construction reduces to the choice of suitable linear operators mapping discrete functions to conforming ones. Such smoothing operators are devised in the forthcoming parts of this work for various finite element spaces

Quasi-optimal nonconforming methods for symmetric elliptic problems. III-discontinuous Galerkin and other interior penalty methods

We devise new variants of the following nonconforming finite element methods: discontinuous Galerkin methods of fixed arbitrary order for the Poisson problem, the Crouzeix-Raviart interior penalty method for linear elasticity, and the quadratic C0 interior penalty method for the biharmonic problem. Each variant differs from the original method only in the discretization of the right-hand side. Before applying the load functional, a linear operator transforms nonconforming discrete test functions into conforming functions such that stability and consistency are improved. The new variants are thus quasi-optimal with respect to an extension of the energy norm. Furthermore, their quasi-optimality constants are uniformly bounded for shape regular meshes and tend to 1 as the penalty parameter increases

An efficient and reliable residual-type a posteriori error estimator for the Signorini problem

We derive a new a posteriori error estimator for the Signorini problem. It generalizes the standard residual-type estimators for unconstrained problems in linear elasticity by additional terms at the contact boundary addressing the non-linearity. Remarkably these additional contact-related terms vanish in the case of so-called full-contact. We prove reliability and efficiency for two- and three-dimensional simplicial meshes. Moreover, we address the case of non-discrete gap functions. Numerical tests for different obstacles and starting grids illustrate the good performance of the a posteriori error estimator in the two- and three-dimensional case, for simplicial as well as for unstructured mixed meshes

A posteriori error estimates with point sources in fractional sobolev spaces

We consider Poisson's equation with a finite number of weighted Dirac masses as a source term, together with its discretization by means of conforming finite elements. For the error in fractional Sobolev spaces, we propose residual-type a posteriori estimators with a specifically tailored oscillation and show that, on two-dimensional polygonal domains, they are reliable and locally efficient. In numerical tests, their use in an adaptive algorithm leads to optimal error decay rates

Suitability of Magnesium Oxide as a Visar Window

Impedance matching of a velocity interferometer for any reflector (VISAR) window to a material under study helps simplify a shock experiment by effectively allowing one to measure an in situ particle velocity. The shock impedance of magnesium oxide (MgO) falls roughly midway between those of sapphire and LiF, two of the most frequently used VISAR window materials. A series of symmetric impact experiments was performed to characterize the suitability of single crystal, (100) oriented magnesium oxide as a VISAR window material. These experiments yielded good results and show the viability of MgO as a VISAR window up to 23 GPa. Results were used to determine window correction factors and, subsequently, to estimate the pressure induced change in index of refraction. In many of the shots in this work we exceeded the Hugoniot elastic limit (HEL) of MgO, and both elastic and plastic waves are evident in the velocity profiles. The presence of both waves within the VISAR window complicates the typical VISAR window correction analysis. Preliminary analysis of the elastic and plastic contributions to the window correction is presented

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

We introduce an iterative method for computing the first eigenpair $(\lambda_{p},e_{p})$ for the $p$-Laplacian operator with homogeneous Dirichlet data as the limit of $(\mu_{q,}u_{q})$ as $q\rightarrow p^{-}$, where $u_{q}$ is the positive solution of the sublinear Lane-Emden equation $-\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 $u_{q}$ to $e_{p}$ is in the $C^{1}$-norm and the rate of convergence of $\mu_{q}$ to $\lambda_{p}$ is at least $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

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

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

Free-Surface Optical Scattering as an Indicator of the Shock-Induced Solid-Liquid Phase Transition in Tin

When highly polished metal surfaces melt upon release after shock loading, they exhibit features that suggest significant surface changes accompany the phase transition. The reflection of light from such surfaces changes from specular (pre-shock) to diffuse upon melting. A familiar manifestation of this phenomenon is the loss of signal light in VISAR measurements, which occurs at pressures high enough to melt the free surface. Unlike many other potential material phase-sensitive diagnostics (e.g., reflectometry, conductivity) that show relatively small (1%–10%) changes, the specularity of reflection provides a more sensitive and definitive indication of the solid-liquid phase transition. Data will be presented that support the hypothesis that specularity changes indicate melt in a way that can be measured easily and unambiguously