A POSTERIORI ERROR CONTROL OF DISCONTINUOUS GALERKIN METHODS FOR ELLIPTIC OBSTACLE PROBLEMS

In this article, we derive an a posteriori error estimator for various discontinuous Galerkin (DG) methods that are proposed in (Wang, Han and Cheng, SIAM J. Numer. Anal., 48: 708-733, 2010) for an elliptic obstacle problem. Using a key property of DG methods, we perform the analysis in a general framework. The error estimator we have obtained for DG methods is comparable with the estimator for the conforming Galerkin (CG) finite element method. In the analysis, we construct a non-linear smoothing function mapping DG finite element space to CG finite element space and use it as a key tool. The error estimator consists of a discrete Lagrange multiplier associated with the obstacle constraint. It is shown for non-over-penalized DG methods that the discrete Lagrange multiplier is uniformly stable on non-uniform meshes. Finally, numerical results demonstrating the performance of the error estimator are presented.

Third post-Newtonian gravitational waveforms for compact binary systems in general orbits: Instantaneous terms

We compute the instantaneous contributions to the spherical harmonic modes of gravitational waveforms from compact binary systems in general orbits up to the third post-Newtonian (PN) order. We further extend these results for compact binaries in quasielliptical orbits using the 3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric orbits. Using the multipolar post-Minkowskian formalism, starting from the different mass and current-type multipole moments, we compute the spin-weighted spherical harmonic decomposition of the instantaneous part of the gravitational waveform. These are terms which are functions of the retarded time and do not depend on the history of the binary evolution. Together with the hereditary part, which depends on the binary's dynamical history, these waveforms form the basis for construction of accurate templates for the detection of gravitational wave signals from binaries moving in quasielliptical orbits.

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 FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.

Anisotropic thermopiezoelectric solids with an elliptic inclusion or a hole under uniform heat flow

The two-dimensional problem of a thermopiezoelectric material containing an elliptic inclusion or a hole subjected to a remote uniform heat flow is studied. Based on the extended Lekhnitskii formulation for thermopiezoelectricity, conformal mapping and Laurent series expansion, the explicit and closed-form solutions are obtained both inside and outside the inclusion (or hole). For a hole problem, the exact electric boundary conditions on the hole surface are used. The results show that the electroelastic fields inside the inclusion or the electric field inside the hole are linear functions of the coordinates. When the elliptic hole degenerates into a slit crack, the electroelastic fields and the intensity factors are obtained. The effect of the heat how direction and the dielectric constant of air inside the crack on the thermal electroelastic fields are discussed. Comparison is made with two special cases of which the closed solutions exist and it is shown that our results are valid.

Elastic fields of interacting elliptic inhomogeneities

This paper presents a method for the calculation of two-dimensional elastic fields in a solid containing any number of inhomogeneities under arbitrary far field loadings. The method called 'pseudo-dislocations method', is illustrated for the solution of interacting elliptic inhomogeneities. It reduces the interacting inhomogeneities problem to a set of linear algebraic equations. Numerical results are presented for a variety of elliptic inhomogeneity arrangements, including the special cases of elliptic holes, cracks and circular inhomogeneities. All these complicated problems can be solved with high accuracy and efficiency.

Rotation Numbers Of Invariant Manifolds Around Unstable Periodic Orbits For The Diamagnetic Kepler Problem

In this paper, a method to construct topological template in terms of symbolic dynamics for the diamagnetic Kepler problem is proposed. To confirm the topological template, rotation numbers of invariant manifolds around unstable periodic orbits in a phase space are taken as an object of comparison. The rotation numbers are determined from the definition and connected with symbolic sequences encoding the periodic orbits in a reduced Poincare section. Only symbolic codes with inverse ordering in the forward mapping can contribute to the rotation of invariant manifolds around the periodic orbits. By using symbolic ordering, the reduced Poincare section is constricted along stable manifolds and a topological template, which preserves the ordering of forward sequences and can be used to extract the rotation numbers, is established. The rotation numbers computed from the topological template are the same as those computed from their original definition.

A method to find unstable periodic orbits for the diamagnetic Kepler problem

A method to determine the admissibility of symbolic sequences and to find the unstable periodic orbits corresponding to allowed symbolic sequences for the diamagnetic Kepler problem is proposed by using the ordering of stable and unstable manifolds. By investigating the unstable periodic orbits up to length 6, a one to one correspondence between the unstable periodic orbits and their corresponding symbolic sequences is shown under the system symmetry decomposition.

Resonant mode analysis in toroidal cavities with elliptic cross-sections

Resonant cavity modes in a torus with elliptical cross section are studied by means of a direct variational method. The nonlinear effects of toroidicity and ellipticity on the frequency of the basic mode are analyzed simply and systematically without the restriction of linear theory. It is shown that the toroidicity effect on the m = 0 transverse magnetic mode is less-than-or-equal-to 11%. The frequency of the mode shifts approximately 11-29% when the elongation of the cross section changes from 1 to 2. The effects of toroidicity and ellipticity differ for each resonant mode.

Magnetohydrodynamics equilibrium of a self-confined elliptic plasma ball

A variational principle is applied to the problem of magnetohydrodynamics (MHD) equilibrium of a self-contained elliptical plasma ball, such as elliptical ball lightning. The principle is appropriate for an approximate solution of partial differential equations with arbitrary boundary shape. The method reduces the partial differential equation to a series of ordinary differential equations and is especially valuable for treating boundaries with nonlinear deformations. The calculations conclude that the pressure distribution and the poloidal current are more uniform in an oblate self-confined plasma ball than that of an elongated plasma ball. The ellipticity of the plasma ball is obviously restricted by its internal pressure, magnetic field, and ambient pressure. Qualitative evidence is presented for the absence of sighting of elongated ball lightning.

Integration property of a wedge-elliptic cone waverider with scramjet engine

The hypersonic waverider forebody is designed in this paper. For the present waverider, the undersurface is carved out as a stream surface of a hypersonic inviscid flow field around wedge-elliptic cone, and the upper surface is assumed to be a freestream surface. A finite-volume code is used to generate the three-dimensional flow field. The leading edge is determined by satisfying the condition that the lip is situated at the intersection line of shocks.

Estimating the girth of fish by applying an elliptic model

A new method for estimating the maximum girth of roundfish is proposed as illustrated; this is based on an elliptic approximation of the cross sections of the fish body. Results derived from a small sample horse mackerel, Trachurus trachurus, suggest that maximum girth estimates based on the elliptic model are more precise than the values estimated by applying a conventional method.