Fredrik Berntsson

1. In many industrial applications, one wishes to determine the temperature on the surface of a body in situations where the surface itself is inaccessible to direct measurements (Beck et al., 1985...

In this paper, we are interested in finding an optimal control support design for controlling vibrations due to pedestrian-bridge interactions. Therefore, we derive the topological derivatives of the proposed functionals using the averaged adjoint approach. A numerical algorithm initialized by these sensitivities is used as a solution strategy. The algorithm is tested numerically for two different cases of initial conditions.

In this paper, an optimal actuator placement problem, with a linear wave equation as the constraint, is considered. In particular, this work presents the framework for finding the best location of actuators depending upon the given initial conditions, and where the dependence on the initial conditions is averaged out. The problem is motivated by the need to control vibrations induced by pedestrian-bridge interactions. An approach based on the shape optimization techniques is used to solve the problem. Specifically, the shape sensitivities involving a cost functional are determined using the averaged adjoint approach. A numerical algorithm based on these sensitivities is used as a solution strategy. Numerical tests illustrate the theoretical results.

We consider an inverse problem for the Helmholtz equation of reconstructing a solution from measurements taken on a segment inside a semi-infinite strip. Homogeneous Neumann conditions are prescribed on both side boundaries of the strip and an unknown Dirichlet condition on the remaining part of the boundary. Additional complexity is that the radiation condition at infinity is unknown. Our aim is to find the unknown function in the Dirichlet boundary condition and the radiation condition. Such problems appear in acoustics to determine acoustical sources and surface vibrations from acoustic field measurements. The problem is split into two sub-problems, a well-posed and an ill-posed problem. We analyse the theoretical properties of both problems; in particular, we show that the radiation condition is described by a stable non-linear problem. The second problem is ill-posed, and we use the Landweber iteration method together with the discrepancy principle to regularize it. Numerical tests show that the approach works well.

We consider the Cauchy problem for the Helmholtz equation with a domain in ℝ d {\mathbb{R}^{d}} , d ≥ 2 {d\geq 2} with N cylindrical outlets to infinity with bounded inclusions in ℝ d - 1 . {\mathbb{R}^{d-1}.} Cauchy data are prescribed on the boundary of the bounded domains and the aim is to find solution on the unbounded part of the boundary. In 1989, Kozlov and Maz’ya [14] proposed an alternating iterative method for solving Cauchy problems associated with elliptic, self-adjoint and positive-definite operators in bounded domains. Different variants of this method for solving Cauchy problems associated with Helmholtz-type operators exists. We consider the variant proposed by Berntsson, Kozlov, Mpinganzima and Turesson (2018) [4] for bounded domains and derive the necessary conditions for the convergence of the procedure in unbounded domains. For the numerical implementation, a finite difference method is used to solve the problem in a simple rectangular domain in ℝ 2 {\mathbb{R}^{2}} that represent a truncated infinite strip. The numerical results shows that by appropriate truncation of the domain and with appropriate choice of the Robin parameters μ 0 {\mu_{0}} and μ 1 {\mu_{1}} , the Robin–Dirichlet alternating iterative procedure is convergent.

We consider the Cauchy problem for the Helmholtz equation defined in a rectangular domain. The Cauchy data are prescribed on a part of the boundary and the aim is to find the solution in the entire domain. The problem occurs in applications related to acoustics and is illposed in the sense of Hadamard. In our work we consider regularizing the problem by introducing a bounded approximation of the second derivative by using Cubic smoothing splines. We derive a bound for the approximate derivative and show how to obtain stability estimates for the method. Numerical tests show that the method works well and can produce accurate results. We also demonstrate that the method can be extended to more complicated domains.

The Cauchy problem for the heat equation is a model of situation where one seeks to compute the temperature, or heat-flux, at the surface of a body by using interior measurements. The problem is well-known to be ill-posed, in the sense that measurement errors can be magnified and destroy the solution, and thus regularization is needed. In previous work it has been found that a method based on approximating the time derivative by a Fourier series works well [Berntsson F. A spectral method for solving the sideways heat equation. Inverse Probl. 1999;15:891–906; Eldén L, Berntsson F, Regińska T. Wavelet and Fourier methods for solving the sideways heat equation. SIAM J Sci Comput. 2000;21(6):2187–2205]. However, in our situation it is not resonable to assume that the temperature is periodic which means that additional techniques are needed to reduce the errors introduced by implicitly making the assumption that the solution is periodic in time. Thus, as an alternative approach, we instead approximate the time derivative by using a cubic smoothing spline. This means avoiding a periodicity assumption which leads to slightly smaller errors at the end points of the measurement interval. The spline method is also shown to satisfy similar stability estimates as the Fourier series method. Numerical simulations shows that both methods work well, and provide comparable accuracy, and also that the spline method gives slightly better results at the ends of the measurement interval.

We study a non-linear operator equation originating from a Cauchy problem for an elliptic equation. The problem appears in applications where surface measurements are used to calculate the temperature below the earth surface. The Cauchy problem is ill-posed and small perturbations to the used data can result in large changes in the solution. Since the problem is non-linear certain assumptions on the coefficients are needed. We reformulate the problem as an non-linear operator equation and show that under suitable assumptions the operator is well-defined. The proof is based on making a change of variables and removing the non-linearity from the leading term of the equation. As a part of the proof we obtain an iterative procedure that is convergent and can be used for evaluating the operator. Numerical results show that the proposed procedure converges faster than a simple fixed point iteration for the equation in the the original variables.