A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

A Remark on the A Posteriori Error Analysis of Discontinuous Galerkin Methods for the Obstacle Problem

We revisit the a posteriori error analysis of discontinuous Galerkin methods for the obstacle problem derived in 25]. Under a mild assumption on the trace of obstacle, we derive a reliable a posteriori error estimator which does not involve min/max functions. A key in this approach is an auxiliary problem with discrete obstacle. Applications to various discontinuous Galerkin finite element methods are presented. Numerical experiments show that the new estimator obtained in this article performs better.

Two player game variant of the Erdos-Szekeres problem

The classical Erdos-Szekeres theorem states that a convex k-gon exists in every sufficiently large point set. This problem has been well studied and finding tight asymptotic bounds is considered a challenging open problem. Several variants of the Erdos-Szekeres problem have been posed and studied in the last two decades. The well studied variants include the empty convex k-gon problem, convex k-gon with specified number of interior points and the chromatic variant. In this paper, we introduce the following two player game variant of the Erdos-Szekeres problem: Consider a two player game where each player playing in alternate turns, place points in the plane. The objective of the game is to avoid the formation of the convex k-gon among the placed points. The game ends when a convex k-gon is formed and the player who placed the last point loses the game. In our paper we show a winning strategy for the player who plays second in the convex 5-gon game and the empty convex 5-gon game by considering convex layer configurations at each step. We prove that the game always ends in the 9th step by showing that the game reaches a specific set of configurations.

A Variant of the Hadwiger-Debrunner (p, q)-Problem in the Plane

Let X be a convex curve in the plane (say, the unit circle), and let be a family of planar convex bodies such that every two of them meet at a point of X. Then has a transversal of size at most . Suppose instead that only satisfies the following ``(p, 2)-condition'': Among every p elements of , there are two that meet at a common point of X. Then has a transversal of size . For comparison, the best known bound for the Hadwiger-Debrunner (p, q)-problem in the plane, with , is . Our result generalizes appropriately for if is, for example, the moment curve.

A posteriori error estimates of discontinuous Galerkin methods for the Signorini problem

A reliable and efficient a posteriori error estimator is derived for a class of discontinuous Galerkin (DG) methods for the Signorini problem. A common property shared by many DG methods leads to a unified error analysis with the help of a constraint preserving enriching map. The error estimator of DG methods is comparable with the error estimator of the conforming methods. Numerical experiments illustrate the performance of the error estimator. (C) 2015 Elsevier B.V. All rights reserved.

Exact Phase Retrieval in Principal Shift-Invariant Spaces

We address the problem of phase retrieval from Fourier transform magnitude spectrum for continuous-time signals that lie in a shift-invariant space spanned by integer shifts of a generator kernel. The phase retrieval problem for such signals is formulated as one of reconstructing the combining coefficients in the shift-invariant basis expansion. We develop sufficient conditions on the coefficients and the bases to guarantee exact phase retrieval, by which we mean reconstruction up to a global phase factor. We present a new class of discrete-domain signals that are not necessarily minimum-phase, but allow for exact phase retrieval from their Fourier magnitude spectra. We also establish Hilbert transform relations between log-magnitude and phase spectra for this class of discrete signals. It turns out that the corresponding continuous-domain counterparts need not satisfy a Hilbert transform relation; notwithstanding, the continuous-domain signals can be reconstructed from their Fourier magnitude spectra. We validate the reconstruction guarantees through simulations for some important classes of signals such as bandlimited signals and piecewise-smooth signals. We also present an application of the proposed phase retrieval technique for artifact-free signal reconstruction in frequency-domain optical-coherence tomography (FDOCT).

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.

A general design method for electric machines using magnetic circuit model considering the flux saturation problem

This paper proposes an analytical approach that is generalized for the design of various types of electric machines based on a physical magnetic circuit model. Conventional approaches have been used to predict the behavior of electric machines but have limitations in accurate flux saturation analysis and hence machine dimensioning at the initial design stage. In particular, magnetic saturation is generally ignored or compensated by correction factors in simplified models since it is difficult to determine the flux in each stator tooth for machines with any slot-pole combinations. In this paper, the flux produced by stator winding currents can be calculated accurately and rapidly for each stator tooth using the developed model, taking saturation into account. This aids machine dimensioning without the need for a computationally expensive finite element analysis (FEA). A 48-slot machine operated in induction and doubly-fed modes is used to demonstrate the proposed model. FEA is employed for verification.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.