Singular Solutions of Protter’s Problem for a Class of 3-D Hyperbolic Equations

Недю Иванов Попиванов, Алексей Йорданов Николов - През 1952 г. М. Протър формулира нови гранични задачи за вълновото уравнение, които са тримерни аналози на задачите на Дарбу в равнината. Задачите са разгледани в тримерна област, ограничена от две характеристични конуса и равнина. Сега, след като са минали повече от 50 години, е добре известно, че за безброй гладки функции в дясната страна на уравнението тези задачи нямат класически решения, а обобщеното решение има силна степенна особеност във върха на характеристичния конус, която е изолирана и не се разпространява по конуса. Тук ние разглеждаме трета гранична задача за вълновото уравнение с младши членове и дясна страна във формата на тригонометричен полином. Дадена е по-нова от досега известната априорна оценка за максимално възможната особеност на решенията на тази задача. Оказва се, че при по-общото уравнение с младши членове възможната сингулярност е от същия ред като при чисто вълновото уравнение.

Polyharmonic Hardy Spaces on the Klein-Dirac Quadric with Application to Polyharmonic Interpolation and Cubature Formulas

2010 Mathematics Subject Classification: Primary 65D30, 32A35, Secondary 41A55.

Interior Boundaries for Degenerate Elliptic Equations of Second Order Some Theory and Numerical Observations

2010 Mathematics Subject Classification: Primary 35J70; Secondary 35J15, 35D05.

Micromechanical insight into the undrained instability of granular materials: deformation characteristics of geomaterials.

There is no agreement between experimental researchers whether the point where a granular material responds with a large change of stresses, strains or excess pore water pressure given a prescribed small input of some of the same variables defines a straight line or a curve in the stress space. This line, known as the instability line, may also vary in shape and position if the onset of instability is measured from drained or undrained triaxial tests. Failure of granular materials, which might be preceded by the onset of instability, is a subject that the geotechnical engineers have to deal with in the daily practice, and generally speaking it is associated to different phenomena observed not only in laboratory tests but also in the field. Examples of this are the liquefaction of loose sands subjected to undrained loading conditions and the diffuse instability under drained loading conditions. This research presents results of DEM simulations of undrained triaxial tests with the aim of studying the influence of stress history and relative density on the onset of instability in granular materials. Micro-mechanical analysis including the evolution of coordination numbers and fabric tensors is performed aiming to gain further insight on the particle-scale interactions that underlie the occurrence of this instability. In addition to provide a greater understanding, the results presented here may be useful as input for macro-scale constitutive models that enable the prediction of the onset of instability in boundary value problems.

Discontinuous Galerkin Methods for the Biharmonic Problem

This work is concerned with the design and analysis of hp-version discontinuous Galerkin (DG) finite element methods for boundary-value problems involving the biharmonic operator. The first part extends the unified approach of Arnold, Brezzi, Cockburn & Marini (SIAM J. Numer. Anal. 39, 5 (2001/02), 1749-1779) developed for the Poisson problem, to the design of DG methods via an appropriate choice of numerical flux functions for fourth order problems; as an example we retrieve the interior penalty DG method developed by Suli & Mozolevski (Comput. Methods Appl. Mech. Engrg. 196, 13-16 (2007), 1851-1863). The second part of this work is concerned with a new a-priori error analysis of the hp-version interior penalty DG method, when the error is measured in terms of both the energy-norm and L2-norm, as well certain linear functionals of the solution, for elemental polynomial degrees $p\ge 2$. Also, provided that the solution is piecewise analytic in an open neighbourhood of each element, exponential convergence is also proven for the p-version of the DG method. The sharpness of the theoretical developments is illustrated by numerical experiments.

The unsteady flow field about a right circular cone in unsteady flight.

"These studies were conducted by the General Electric Company, Reentry Systems Department, for the Stability and Control Section of the Flight Dynamics Laboratory of the Air Force Research and Technology Division."

Construcción y evaluación de un portafolio estructural, para obligatorias retiro programado, utilizando la metodología de frontera eficiente de Markowitz y la Teoría de Momentum

Este proyecto de investigación construye y evalúa la asignación de activos para el portafolio de Pensión Obligatoria de Retiro Programado, el cual atiende los retiros a los que un pensionado tiene derecho a través de su mesada pensional, utilizando el modelo de Frontera Eficiente de Markowitz, en combinación con la teoría de Momentum -- Para la ejecución del modelo se determinaron los activos de inversión admisibles en el Régimen de Inversiones -- Posteriormente, se construyen las matrices de rentabilidades, de restricciones, de varianzas y covarianzas, las cuales constituyen los insumos para ejecutar el modelo de optimización de portafolios de Markowitz -- A continuación, se realiza la selección de los portafolios obtenidos, teniendo en cuenta el nivel de volatilidad que el portafolio de Obligatorias Retiro Programado debe presentar; lo anterior, con el fin de cumplir con el objetivo de preservación del capital en la cuenta individual del pensionado, de manera que se pueda atender, de acuerdo a su esperanza de vida y la de sus beneficiarios, el pago de las mesadas pensionales que le correspondan -- El resultado obtenido corresponde a una asignación, en gran parte, en activos de Renta Fija expedidos por el Gobierno Nacional (TES), tanto en tasa fija como en tasa indexada a la UVR -- Adicionalmente, el modelo de optimización sugiere participaciones en activos de renta variable y, particularmente, no asigna recursos representativos en títulos de deuda privada indexados al IPC -- Esta investigación puede ser útil al momento de diseñar un portafolio base para Obligatorias Retiro Programado que, bajo una administración pasiva, permita cumplir el objetivo de otorgar a los pensionados una mesada para satisfacer las necesidades básicas

Introducing the Logarithmic finite element method: a geometrically exact planar Bernoulli beam element

We propose a novel finite element formulation that significantly reduces the number of degrees of freedom necessary to obtain reasonably accurate approximations of the low-frequency component of the deformation in boundary-value problems. In contrast to the standard Ritz–Galerkin approach, the shape functions are defined on a Lie algebra—the logarithmic space—of the deformation function. We construct a deformation function based on an interpolation of transformations at the nodes of the finite element. In the case of the geometrically exact planar Bernoulli beam element presented in this work, these transformation functions at the nodes are given as rotations. However, due to an intrinsic coupling between rotational and translational components of the deformation function, the formulation provides for a good approximation of the deflection of the beam, as well as of the resultant forces and moments. As both the translational and the rotational components of the deformation function are defined on the logarithmic space, we propose to refer to the novel approach as the “Logarithmic finite element method”, or “LogFE” method.

Eigenvalue problems associated with Korn's inequalities in the theory of elasticity

Interest in the possible applications of a priori inequalities in linear elasticity theory motivated the present investigation. Korn's inequality under various side conditions is considered, with emphasis on the Korn's constant. In the "second case" of Korn's inequality, a variational approach leads to an eigenvalue problem; it is shown that, for simply-connected two-dimensional regions, the problem of determining the spectrum of this eigenvalue problem is equivalent to finding the values of Poisson's ratio for which the displacement boundary-value problem of linear homogeneous isotropic elastostatics has a non-unique solution.

Previous work on the uniqueness and non-uniqueness issue for the latter problem is examined and the results applied to the spectrum of the Korn eigenvalue problem. In this way, further information on the Korn constant for general regions is obtained.

A generalization of the "main case" of Korn's inequality is introduced and the associated eigenvalue problem is a gain related to the displacement boundary-value problem of linear elastostatics in two dimensions.

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.

A high-wavenumber boundary-element method for an acoustic scattering problem

In this paper we show stability and convergence for a novel Galerkin boundary element method approach to the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data. This problem models, for example, outdoor sound propagation over inhomogeneous flat terrain. To achieve a good approximation with a relatively low number of degrees of freedom we employ a graded mesh with smaller elements adjacent to discontinuities in impedance, and a special set of basis functions for the Galerkin method so that, on each element, the approximation space consists of polynomials (of degree $\nu$) multiplied by traces of plane waves on the boundary. In the case where the impedance is constant outside an interval $[a,b]$, which only requires the discretization of $[a,b]$, we show theoretically and experimentally that the $L_2$ error in computing the acoustic field on $[a,b]$ is ${\cal O}(\log^{\nu+3/2}|k(b-a)| M^{-(\nu+1)})$, where $M$ is the number of degrees of freedom and $k$ is the wavenumber. This indicates that the proposed method is especially commendable for large intervals or a high wavenumber. In a final section we sketch how the same methodology extends to more general scattering problems.

Numerical solutions of elliptic inverse problems via the equation error method

To estimate a parameter in an elliptic boundary value problem, the method of equation error chooses the value that minimizes the error in the PDE and boundary condition (the solution of the BVP having been replaced by a measurement). The estimated parameter converges to the exact value as the measured data converge to the exact value, provided Tikhonov regularization is used to control the instability inherent in the problem. The error in the estimated solution can be bounded in an appropriate quotient norm; estimates can be derived for both the underlying (infinite-dimensional) problem and a finite-element discretization that can be implemented in a practical algorithm. Numerical experiments demonstrate the efficacy and limitations of the method.

Problems for P-Monge-Ampere Equations

2010 Mathematics Subject Classification: 35A23, 35B51, 35J96, 35P30, 47J20, 52A40.

Monte Carlo Algorithms for Linear Problems

MSC Subject Classification: 65C05, 65U05.

Riemann-Hilbert problems for poly-Hardy space on the unit ball

In this paper, we focus on a Riemann–Hilbert boundary value problem (BVP) with a constant coefficients for the poly-Hardy space on the real unit ball in higher dimensions. We first discuss the boundary behaviour of functions in the poly-Hardy class. Then we construct the Schwarz kernel and the higher order Schwarz operator to study Riemann–Hilbert BVPs over the unit ball for the poly- Hardy class. Finally, we obtain explicit integral expressions for their solutions. As a special case, monogenic signals as elements in the Hardy space over the unit sphere will be reconstructed in the case of boundary data given in terms of functions having values in a Clifford subalgebra. Such monogenic signals represent the generalization of analytic signals as elements of the Hardy space over the unit circle of the complex plane.