This article presents the principal results of the doctoral thesis “Direct Operational Methods in the Environment of a Computer Algebra System” by Margarita Spiridonova (Institute of mathematics and Informatics, BAS), successfully defended before the Specialised Academic Council for Informatics and Mathematical Modelling on 23 March, 2009.The presented research is related to the operational calculus approach and its representative applications. Operational methods are considered, as well as their...

This paper is devoted to the numerical solution of nonlinear elliptic partial differential equations. Such problems describe various phenomena in science. An approach that exploits Hilbert space theory in the numerical study of elliptic PDEs is the idea of preconditioning operators. In this survey paper we briefly summarize the main lines of this theory with various applications.

Let ${𝒯}_h$ be a triangulation of a bounded polygonal domain $\Omega \subset {\Re }^2$, ${ℒ}_h$ the space of the functions from $C\left(\overline{\Omega }\right)$ linear on the triangles from ${𝒯}_h$ and ${\Pi }_h$ the interpolation operator from $C\left(\overline{\Omega }\right)$ to ${ℒ}_h$. For a unit vector $z$ and an inner vertex $a$ of ${𝒯}_h$, we describe the set of vectors of coefficients such that the related linear combinations of the constant derivatives $\partial {\Pi }_h\left(u\right)/\partial z$ on the triangles surrounding $a$ are equal to $\partial u/\partial z\left(a\right)$ for all polynomials $u$ of the total degree less than or equal to two. Then we prove that, generally, the values of the...

In this paper, we examine a particular class of singularly perturbed convection-diffusion problems with a discontinuous coefficient of the convective term. The presence of a discontinuous convective coefficient generates a solution which mimics flow moving in opposing directions either side of some flow source. A particular transmission condition is imposed to ensure that the differential operator is stable. A piecewise-uniform Shishkin mesh is combined with a monotone finite difference operator...

A pair trade is a portfolio consisting of a long position in one asset and a short position in another, and it is a widely used investment strategy in the financial industry. Recently, Ekström, Lindberg, and Tysk studied the problem of optimally closing a pair trading strategy when the difference of the two assets is modelled by an Ornstein-Uhlenbeck process. In the present work the model is generalized to also include jumps. More precisely, we assume that the difference between the assets is an...

We propose a general approach for the numerical approximation of optimal control problems governed by a linear advection–diffusion equation, based on a stabilization method applied to the lagrangian functional, rather than stabilizing the state and adjoint equations separately. This approach yields a coherently stabilized control problem. Besides, it allows a straightforward a posteriori error estimate in which estimates of higher order terms are needless. Our a posteriori estimates stems from splitting...

We propose a general approach for the numerical approximation of optimal control problems governed by a linear advection–diffusion equation, based on a stabilization method applied to the Lagrangian functional, rather than stabilizing the state and adjoint equations separately. This approach yields a coherently stabilized control problem. Besides, it allows a straightforward a posteriori error estimate in which estimates of higher order terms are needless. Our a posteriori estimates stems from...

We study the numerical aspect of the optimal control of problems governed by a linear elliptic partial differential equation (PDE). We consider here the gas flow in porous media. The observed variable is the flow field we want to maximize in a given part of the domain or its boundary. The control variable is the pressure at one part of the boundary or the discharges of some wells located in the interior of the domain. The objective functional is a balance between the norm of the flux in the observation...

We consider an optimal control problem for the three-dimensional non-linear Primitive Equations of the ocean in a vertically bounded and horizontally periodic domain. We aim to reconstruct the initial state of the ocean from Lagrangian observations. This inverse problem is formulated as an optimal control problem which consists in minimizing a cost function representing the least square error between Lagrangian observations and their model counterpart, plus a regularization term. This paper proves...

The one-dimensional steady-state convection-diffusion problem for the unknown temperature $y\left(x\right)$ of a medium entering the interval $(a,b)$ with the temperature $y_{min}$ and flowing with a positive velocity $v\left(x\right)$ is studied. The medium is being heated with an intensity corresponding to $y_{max}-y\left(x\right)$ for a constant $y_{max}>y_{min}$. We are looking for a velocity $v\left(x\right)$ with a given average such that the outflow temperature $y\left(b\right)$ is maximal and discuss the influence of the boundary condition at the point $b$ on the “maximizing” function $v\left(x\right)$.

We consider the original DG method for solving the advection-reaction equations with arbitrary velocity in $d$ space dimensions. For triangulations satisfying the flow condition, we first prove that the optimal convergence rate is of order $k+1$ in the $L_2$-norm if the method uses polynomials of order $k$. Then, a very simple derivative recovery formula is given to produce an approximation to the derivative in the flow direction which superconverges with order $k+1$. Further we consider a residual-based a posteriori...