A computational procedure for estimation of an unknown coefficient in an inverse boundary value problem
Introduction
Inverse problems play nowadays an important role in the solution of a number of practical problems. The use of inverse methods represent a new research direction, where the results obtained from numerical simulation and from experiments are not simply compared a posteriori, but a closed synergism exist between experimental and theoretical researches during a course of study, in order to obtain the maximum of information regarding the physical problem under consideration. Most of the methods for the solutions of inverse problems, which are currently in common use, were formalized in the last four decades in terms of their capabilities to treat ill-posed unstable problems. The basis of such formal methods resides on the idea of reformulating an inverse problem in terms of an approximate well-posed problem, by utilizing some kind of regularization techniques [1], [2], [3], [4], [5], [6], [7], [8], [9].
In [1], the in-cylinder pressure in an internal combustion engine was modeled as a function of crank angle. This model was stated as a boundary value problem of second-order and it was studied as an inverse problem. In this work, we investigate the use of conjugate gradient method together with an adjoint problem formulation for the estimation of an unknown function that appeared in the coefficients of a boundary value problem of second-order. The inverse problem considered in this work is solved by using a function estimation approach [4], [5], [6], where no information is a priori available regarding the forms of the unknown function, except for the functional space that they belong to. It is assumed that the unknown function belong to Hilbert space of square integrable function in the spatial domain of interest (L2[D]).
This article is organized as follows: (i) problem definition, (ii) variational formulation, (iii) iterative regularization procedure, (iv) stopping criterion, and (v) numerical results.
Section snippets
Problem definition
In this study we consider the following boundary value problem:where p(x) ∈ C1[a, b] and q(x),g(x) ∈ L2[a,b].
Notice that in the direct problem, all coefficients p(x), q(x), and g(x) are regarded as known quantities, so that a direct problem is concerned with the computation of y(x).
For the inverse problem of interest here, the function q(x) is regarded as unknown. Such function will be estimated by using measurements of y(x) taken at the locations xi
Variational formulation
For estimation of the unknown function q(x), we use the variational formulation of the inverse problem of interest. For this purpose we make use of minimization procedure involving the following objective functionalThe estimated dependent variable y(xi; q) is obtained from the solution of the direct problem at the measurement positions xi, i = 1, 2, …, N with the estimates for q(x).
The conjugate gradient method with an adjoint problem is used for the minimization of the
Iterative regularization procedure
For the estimation of q(x), the iterative procedure of the conjugate gradient method is written as follows [4], [5], [6]:where dn(x) is the direction of descent, β is the search step size, and n is the number of iterations.
For the iterative procedure, the direction of descent is obtained as a linear combination of the gradient direction with directions of descent of the previous iterations. The direction of descent for the conjugate gradient method can be written as
Stopping criterion
In Section 3, it was pointed out that the conjugate gradient method of function estimation belongs to the class of iterative regularization methods. In such class of methods, the stopping criterion for the computational procedure is used as a regularization parameter, so that sufficiently accurate and smooth solution is obtained for the unknown functions. For the stopping criterion, we illustrate in this work the use of the discrepancy principle [4], [5], [6], [8], [9].
With the use of the
Numerical results
The accuracy of the present solution approach is examined by using simulated measurements of y(x) containing random errors in the inverse analysis in the following examples.
For the test cases examined, the direct, sensitivity and adjoint problems were numerically solved by using MATLAB software. The simulated measurements of y(x) in the following examples contained random errors with standard deviation σ = 0.01Ymax, where Ymax is the maximum absolute value of the measured variable. Example 1 For the first
Conclusion
In this article, a function estimation approach based on the conjugate gradient method is applied for the estimation of the function that appeared in the coefficient of a boundary value problem of second kind. The accuracy of the proposed solution method is addressed by using simulated measurements containing random errors.
This method can be applied for a lot of industrial and engineering problems that deals with boundary value problems with unknown boundaries, such as the inverse problem that
References (9)
- et al.
Numerical study of in-cylinder pressure in an internal combustion engine
Appl. Math. Comput.
(2005) - et al.
Numerical approximation of solution of an inverse heat conduction problem based on Legendre polynomials
Appl. Math. Comput.
(2006) - et al.
Numerical solution of inverse heat conduction problem with nonstationary measurements
Appl. Math. Comput.
(2005) - et al.
Extreme Method for Solving Ill-posed Problems with Application to Inverse Heat Transfer Problems
(1995)
Cited by (8)
An effective 3-D inverse procedure to retrieve cooling conditions in an aluminium alloy continuous casting problem
2010, Applied Thermal EngineeringCitation Excerpt :Because of the difficulties mentioned above and due to high costs of the experimental measurements, computational tests are usually carried out with numerically simulated measurements, which is a typical procedure applied for this type of problems. For example, Fakhraie et al. [31] used measurements burdened with normally distributed random errors for an estimation of the unknown coefficients in the inverse boundary problem. The problem in question was a theoretical one and it was defined as a second-order differential equation.
A numerical procedure for estimation of the melt depth in laser material processing
2009, Optics and Laser TechnologyNumerical Analysis of Direct and Inverse Problems for a Fractional Parabolic Integro-Differential Equation
2023, Fractal and FractionalParameter Estimation Inspired by Temperature Measurements for a Chemotactic Model of Honeybee Thermoregulation
2023, Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics)Numerical Optimization Identification of a Keller-Segel Model for Thermoregulation in Honey Bee Colonies in Winter
2023, Communications in Computer and Information ScienceGradient Optimization in Reconstruction of the Diffusion Coefficient in a Time Fractional Integro-Differential Equation of Pollution in Porous Media
2023, Communications in Computer and Information Science