Construction of sundials via vectors

A method based on oblique projection is presented for construction of sundials. The derived formulas are classical, but usage of vectors and projections renders a coherent presentation rather than a number of special cases. The presented work is aimed to be useful for those taking a beginning module on vector algebra

Uniqueness and counterexamples in some inverse source problems

Uniqueness of a solution is investigated for some inverse source problems arising in linear parabolic equations. We prove new uniqueness results formulated in Theorems 3.1 and 3.2. We also show optimality of the conditions under which uniqueness holds by explicitly constructing counterexamples, that is by constructing more than one solution in the case when the conditions for uniqueness are violated

An iterative regularizing method for an incomplete boundary data problem for the biharmonic equation

An incomplete boundary data problem for the biharmonic equation is considered, where the displacement is known throughout the boundary of the solution domain whilst the normal derivative and bending moment are specified on only a portion of the boundary. For this inverse ill‐posed problem an iterative regularizing method is proposed for the stable data reconstruction on the underspecified boundary part. Convergence is proven by showing that the method can be written as a Landweber‐type procedure for an operator formulation of the incomplete data problem. This reformulation renders a stopping rule, the discrepancy principle, for terminating the iterations in the case of noisy data. Uniqueness of a solution to the considered problem is also shown

A boundary integral equation method for numerical solution of parabolic and hyperbolic Cauchy problems

We present a unified boundary integral approach for the stable numerical solution of the ill-posed Cauchy problem for the heat and wave equation. The method is based on a transformation in time (semi-discretisation) using either the method of Rothe or the Laguerre transform, to generate a Cauchy problem for a sequence of inhomogenous elliptic equations; the total entity of sequences is termed an elliptic system. For this stationary system, following a recent integral approach for the Cauchy problem for the Laplace equation, the solution is represented as a sequence of single-layer potentials invoking what is known as a fundamental sequence of the elliptic system thereby avoiding the use of volume potentials and domain discretisation. Matching the given data, a system of boundary integral equations is obtained for finding a sequence of layer densities. Full discretisation is obtained via a Nyström method together with the use of Tikhonov regularization for the obtained linear systems. Numerical results are included both for the heat and wave equation confirming the practical usefulness, in terms of accuracy and resourceful use of computational effort, of the proposed approach

On the numerical solution of a Cauchy problem in an elastostatic half-plane with a bounded inclusion

We propose an iterative procedure for the inverse problem of determining the displacement vector on the boundary of a bounded planar inclusion given the displacement and stress fields on an infinite (planar) line-segment. At each iteration step mixed boundary value problems in an elastostatic half-plane containing the bounded inclusion are solved. For efficient numerical implementation of the procedure these mixed problems are reduced to integral equations over the bounded inclusion. Well-posedness and numerical solution of these boundary integral equations are presented, and a proof of convergence of the procedure for the inverse problem to the original solution is given. Numerical investigations are presented both for the direct and inverse problems, and these results show in particular that the displacement vector on the boundary of the inclusion can be found in an accurate and stable way with small computational cost

On non-denseness for a method of fundamental solutions with source points fixed in time for parabolic equations

Linear combinations of fundamental solutions to the parabolic heat equation with source points fixed in time is investigated. The open problem whether these linear combinations generate a dense set in the space of square integrable functions on the lateral boundary of a space-time cylinder, is settled in the negative. Linear independence of the set of fundamental solutions is shown to hold. It is outlined at the end, for a particular example, that such linear combinations constitute a linearly independent and dense set in the space of square integrable functions on the upper top part (where time is fixed) of the boundary of this space-time cylinder

On non-denseness for a method of fundamental solutions with source points fixed in time for parabolic equations

Linear combinations of fundamental solutions to the parabolic heat equation with source points fixed in time is investigated. The open problem whether these linear combinations generate a dense set in the space of square integrable functions on the lateral boundary of a space-time cylinder, is settled in the negative. Linear independence of the set of fundamental solutions is shown to hold. It is outlined at the end, for a particular example, that such linear combinations constitute a linearly independent and dense set in the space of square integrable functions on the upper top part (where time is fixed) of the boundary of this space-time cylinder

Numerical simulations in 3-dimensions of reaction–diffusion models for brain tumour growth

We work with a well-known model of reaction–diffusion type for brain tumour growth and accomplish full 3-dimensional (3d) simulations of the tumour in time on two types of imaging data, the 3d Shepp–Logan head phantom image and an MRI T1-weighted brain scan from the Internet Brain Segmentation Repository. The source term is such that we have logistic growth. These simulations are obtained using standard finite difference approximations with novel calculations to increase speed and accuracy. Moreover, biological background to the model, its well-posedness together with a variational formulation are given. The variational formulation enable the feasibility of different derivations and modifications of the model

An iterative method for the Cauchy problem for second-order elliptic equations

The problem of reconstructing the solution to a second-order elliptic equation in a doubly-connected domain from knowledge of the solution and its normal derivative on the outer part of the boundary of the solution domain, that is from Cauchy data, is considered. An iterative method is given to generate a stable numerical approximation to this inverse ill-posed problem. The procedure is physically feasible in that boundary data is updated with data of the same type in the iterations, meaning that Dirichlet values is updated with Dirichlet values from the previous step and Neumann values by Neumann data. Proof of convergence and stability are given by showing that the proposed method is an extension of the Landweber method for an operator equation reformulation of the Cauchy problem. Connection with the alternating method is discussed. Numerical examples are included confirming the feasibility of the suggested approach

Properties of a method of fundamental solutions for the parabolic heat equation

We show that a set of fundamental solutions to the parabolic heat equation, with each element in the set corresponding to a point source located on a given surface with the number of source points being dense on this surface, constitute a linearly independent and dense set with respect to the standard inner product of square integrable functions, both on lateral- and time-boundaries. This result leads naturally to a method of numerically approximating solutions to the parabolic heat equation denoted a method of fundamental solutions (MFS). A discussion around convergence of such an approximation is included