Inverse problem in anomalous diffusion with uncertainty propagation
Journal Title: Computer Assisted Methods in Engineering and Science - Year 2014, Vol 21, Issue 3
Abstract
Recently, Bevilacqua, Galeão and co-workers have developed a new analytical formulation for the simulation of diffusion with retention phenomena. This new formulation aims at the reduction of all diffusion processes with retention to a unifying model that can adequately simulate the retention effect. This model may have relevant applications in a number of different areas such as population spreading with partial hold up of the population to guarantee territorial domain chemical reactions inducing adsorption processes and multiphase flow through porous media. In this new formulation a discrete approach is firstly formulated taking into account a control parameter which represents the fraction of particles that are able to diffuse. The resulting governing equation for the modelling of diffusion with retention in a continuum medium requires a fourth-order differential term. Specific experimental techniques, together with an appropriate inverse analysis, need to be determined to characterize complementary parameters. The present work investigates an inverse problem which does not allow for simultaneous estimation of all model parameter. In addition a two-step characterization procedure is proposed: in the first step the diffusion coefficient is estimated and in the second one the complementary parameters are estimated. In this paper, it is assumed that the first step is already completed and the diffusion coefficient is known with a certain degree of reliability. Therefore, this work is aimed at investigating the confidence intervals of the complementary parameters estimates considering both the uncertainties due to measurement errors in the experimental data and due to the uncertainty propagation of the estimated value of the diffusion coefficient. The inverse problem solution is carried out through the maximum likelihood approach, with the minimization problem solved with the Levenberg-Marquardt method, and the estimation of the confidence intervals is carried out through the Monte Carlo analysis.
Authors and Affiliations
Luciano G. Silva, Diego C. Knupp, Luiz Bevilacqua, Antônio J. Silva Neto, Augusto C. N. R. Galeão
Keywords
Related Articles
- EP ID EP73503
- DOI -
- Views 183
- Downloads 0
Numerical estimation of internal stress relieving in destructive test. (Received in the final form January 18, 2010)
The paper deals with the analysis of residual stress fields in the riveted joint and the estimation of the internal stress magnitude releasing by partial and complete removing of the rivet material. Stress relieving caus...
Problems of the equilibrium of a rigid body and mechanical systems. (Received in the final form August 12, 2009)
In this article one of the greatest generalized methods for establishing the equilibrium equations of a rigid body and the set of rigid bodies is proposed. It is related to six equations of moments of force about six the...
Numerical analysis of radiant warmer
The main objective of this paper is to recognise the heat transfer phenomena between an infant placed in a radiant warmer and the surrounding environment. The influence of the parameters in the computational time during...
Application of sensitivity analysis in microscale heat transfer
In the paper, the thin metal film subjected to the ultrashort laser pulse has been analyzed. The heat conduction in the domain considered has been described by two-temperature model consisting of the system of two couple...
SVD as a preconditioner in nonlinear optimization
Finding a solution of nonlinear constrained optimization problem may be very computer resources consuming, regardless of solution method adopted. A conceptually simple preconditioning procedure, based on singular value d...