A numerical study of variable coefficient elliptic Cauchy problem via projection method

In this paper, we investigate a numerical method for the solution of an inverse problem of recovering lacking data on some part of the boundary of a domain from the Cauchy data on other part for a variable coefficient elliptic Cauchy problem. In the process, the Cauchy problem is transformed into the problem of solving a compact linear operator equation. As a remedy to the ill-posedness of the problem, we use a projection method which allows regularization solely by discretization. The discretization level plays the role of regularization parameter in the case of projection method. The balancing principle is used for the choice of an appropriate discretization level. Several numerical examples show that the method produces a stable good approximate solution.

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.

Construction of Singer Subgroup Orbit Codes Based on Cyclic Difference Sets

A recent approach for the construction of constant dimension subspace codes, designed for error correction in random networks, is to consider the codes as orbits of suitable subgroups of the general linear group. In particular, a cyclic orbit code is the orbit of a cyclic subgroup. Hence a possible method to construct large cyclic orbit codes with a given minimum subspace distance is to select a subspace such that the orbit of the Singer subgroup satisfies the distance constraint. In this paper we propose a method where some basic properties of difference sets are employed to select such a subspace, thereby providing a systematic way of constructing cyclic orbit codes with specified parameters. We also present an explicit example of such a construction.

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.

Observation of correlated spin-orbit order in a strongly anisotropic quantum wire system

Quantum wires with spin-orbit coupling provide a unique opportunity to simultaneously control the coupling strength and the screened Coulomb interactions where new exotic phases of matter can be explored. Here we report on the observation of an exotic spin-orbit density wave in Pb-atomic wires on Si(557) surfaces by mapping out the evolution of the modulated spin-texture at various conditions with spin-and angle-resolved photoelectron spectroscopy. The results are independently quantified by surface transport measurements. The spin polarization, coherence length, spin dephasing rate and the associated quasiparticle gap decrease simultaneously as the screened Coulomb interaction decreases with increasing excess coverage, providing a new mechanism for generating and manipulating a spin-orbit entanglement effect via electronic interaction. Despite clear evidence of spontaneous spin-rotation symmetry breaking and modulation of spin-momentum structure as a function of excess coverage, the average spin polarization over the Brillouin zone vanishes, indicating that time-reversal symmetry is intact as theoretically predicted.

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.

Edge states, spin transport, and impurity-induced local density of states in spin-orbit coupled graphene

We study graphene, which has both spin-orbit coupling (SOC), taken to be of the Kane-Mele form, and a Zeeman field induced due to proximity to a ferromagnetic material. We show that a zigzag interface of graphene having SOC with its pristine counterpart hosts robust chiral edge modes in spite of the gapless nature of the pristine graphene; such modes do not occur for armchair interfaces. Next we study the change in the local density of states (LDOS) due to the presence of an impurity in graphene with SOC and Zeeman field, and demonstrate that the Fourier transform of the LDOS close to the Dirac points can act as a measure of the strength of the spin-orbit coupling; in addition, for a specific distribution of impurity atoms, the LDOS is controlled by a destructive interference effect of graphene electrons which is a direct consequence of their Dirac nature. Finally, we study transport across junctions, which separates spin-orbit coupled graphene with Kane-Mele and Rashba terms from pristine graphene both in the presence and absence of a Zeeman field. We demonstrate that such junctions are generally spin active, namely, they can rotate the spin so that an incident electron that is spin polarized along some direction has a finite probability of being transmitted with the opposite spin. This leads to a finite, electrically controllable, spin current in such graphene junctions. We discuss possible experiments that can probe our theoretical predictions.

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.

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.

First-principles study of the electronic structure and dynamical electronic excitations of strong spin-orbit metal Pb: bulk, surface and nanosized films

En la presente tesis se ha realizado el estudio de primeros principios (esto es, sinhacer uso de parámetros ajustables) de la estructura electrónica y la dinámica deexcitaciones electrónicas en plomo, tanto en volumen como en superficie y en formade películas de espesor nanométrico. Al presentar el plomo un número atómico alto(82), deben tenerse en cuenta los efectos relativistas. Con este fin, el doctorando haimplementado el acoplo espín-órbita en los códigos computacionales que hanrepresentado la principal herramienta de trabajo.En volumen, se han encontrado fuertes efectos relativistas asi como de lalocalización de los electrones, tanto en la respuesta dieléctrica (excitacioneselectrónicas colectivas) como en el tiempo de vida de electrones excitados. Lacomparación de nuestros resultados con medidas experimentales ha ayudado aprofundizar en dichos efectos.En el estudio de las películas a escala nanométrica se han hallado fuertes efectoscuánticos debido al confinamiento de los estados electrónicos. Dichos efectos semanifiestan tanto en el estado fundamental (en acuerdo con estudiosexperimentales), como en la respuesta dieléctrica a través de la aparición y dinámicade plasmones de diversas características. Los efectos relativistas, a pesar de no serimportantes en la estructura electrónica de las películas, son los responsables de ladesaparación del plasmón de baja energía en nuestros resultados.