A multigrid method for distributed parameter estimation problems.
In this paper we consider the problem of detecting pollution in some non linear parabolic systems using the sentinel method. For this purpose we develop and analyze a new approach to the discretization which pays careful attention to the stability of the solution. To illustrate convergence properties we give some numerical results that present good properties and show new ways for building discrete sentinels.
We propose an adaptive finite element method for the solution of a coefficient inverse problem of simultaneous reconstruction of the dielectric permittivity and magnetic permeability functions in the Maxwell's system using limited boundary observations of the electric field in 3D. We derive a posteriori error estimates in the Tikhonov functional to be minimized and in the regularized solution of this functional, as well as formulate the corresponding adaptive algorithm. Our numerical experiments...
In specific fields of research such as preservation of historical buildings, medical imaging, geophysics and others, it is of particular interest to perform only a non-intrusive boundary measurements. The idea is to obtain comprehensive information about the material properties inside the considered domain while keeping the test sample intact. This paper is focused on such problems, i.e. synthesizing a physical model of interest with a boundary inverse value technique. The forward model is represented...
Development of engineering structures and technologies frequently works with advanced materials, whose properties, because of their complicated microstructure, cannot be predicted from experience, unlike traditional materials. The quality of computational modelling of relevant physical processes, based mostly on the principles of classical thermomechanics, is conditioned by the reliability of constitutive relations, coming from simplified experiments. The paper demonstrates some possibilities of...
In this article we transform a large class of parabolic inverse problems into a nonclassical parabolic equation whose coefficients consist of trace type functionals of the solution and its derivatives subject to some initial and boundary conditions. For this nonclassical problem, we study finite element methods and present an immediate analysis for global superconvergence for these problems, on basis of which we obtain a posteriori error estimators.
Inverse three-dimensional heat conduction problems devoted to heating source localization are ill posed. Identification can be performed using an iterative regularization method based on the conjugate gradient algorithm. Such a method is usually implemented off-line, taking into account observations (temperature measurements, for example). However, in a practical context, if the source has to be located as fast as possible (e.g., for diagnosis), the observation horizon has to be reduced. To this...
We discuss the numerical modeling of heat exchange between the infiltrated water and porous media matrix. An unsaturated-saturated flow is considered with boundary conditions reflecting the external driven forces. The developed numerical method is efficient and can be used for solving the inverse problems concerning determination of transmission coefficients for heat energy exchange inside and also on the boundary of porous media. Numerical experiments support our method.
The time-dependent intensity of a UV-photon source, located inside an interstellar cloud, is determined by formulating and solving an inverse problem for the integro-differential transport equation of photons in a one-dimensional slab. Starting from a discretizazion of the forward problem, an iterative procedure is used to compute the values of the source intensity at increasing values of the time.
The design of an experiment, e.g., the setting of initial conditions, strongly influences the accuracy of the process of determining model parameters from data. The key concept relies on the analysis of the sensitivity of the measured output with respect to the model parameters. Based on this approach we optimize an experimental design factor, the initial condition for an inverse problem of a model parameter estimation. Our approach, although case independent, is illustrated at the FRAP (Fluorescence...
We present the solution of some inverse problems for one-dimensional free boundary problems of oxygen consumption type, with a semilinear convection-diffusion-reaction parabolic equation. Using a fixed domain transformation (Landau’s transformation) the direct problem is reduced to a system of ODEs. To minimize the objective functionals in the inverse problems, we approximate the data by a finite number of parameters with respect to which automatic differentiation is applied.
In this paper, we consider a 2nd order semilinear parabolic initial boundary value problem (IBVP) on a bounded domain , with nonstandard boundary conditions (BCs). More precisely, at some part of the boundary we impose a Neumann BC containing an unknown additive space-constant , accompanied with a nonlocal (integral) Dirichlet side condition. We design a numerical scheme for the approximation of a weak solution to the IBVP and derive error estimates for the approximation of the solution and...
A class of variational data assimilation problems on reconstructing the initial-value functions is considered for the models governed by quasilinear evolution equations. The optimality system is reduced to the equation for the control function. The properties of the control equation are studied and the solvability theorems are proved for linear and quasilinear data assimilation problems. The iterative algorithms for solving the problem are formulated and justified.
A class of variational data assimilation problems on reconstructing the initial-value functions is considered for the models governed by quasilinear evolution equations. The optimality system is reduced to the equation for the control function. The properties of the control equation are studied and the solvability theorems are proved for linear and quasilinear data assimilation problems. The iterative algorithms for solving the problem are formulated and justified.