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 175
- Downloads 0
The analysis of long-term borehole heat exchanger system exploitation
The key issue in designing borehole heat exchangers (BHE) is the long-term performance of the ground source heat pump (GSHP) systems. The performance directly reflects economic profitability and depends on a large number...
Numerical investigations of the convergence of a remeshing algorithm on an example of subsonic flow
The main goal of the paper is to analyze convergence of a remeshing scheme evaluated by the author [8] on the example of a potential flow around a profile. It is assumed that flow is stationary, irrotational, inviscid an...
Novelty detection based on elastic wave signals measured by different techniques
The paper discusses the results of laboratory experiments in which three independent measurement techniques were compared: a digital oscilloscope, phased array acquisition system, a laser vibrometer 3D. These techniques...
A block sparse shared-memory multifrontal finite element solver for problems of structural mechanics. (Received in the final form July 17, 2009)
The presented method is used in finite-element analysis software developed for multicore and multiprocessor shared-memory computers, or it can be used on single-processor personal computers under the operating systems Wi...
Eigenvalue analysis for high telecommunication towers with lognormal stiffness by the response function method and SFEM. (Received in the final form February 28, 2010)
The main aim of this paper is to demonstrate the application of the generalized stochastic perturbation technique to model the lognormal random variables in structural mechanics. This is done to study probabilistic chara...