No Arabic abstract
The estimate of coefficients of the Convection-Diffusion Equation (CDE) from experimental measurements belongs in the category of inverse problems, which are known to come with issues of ill-conditioning or singularity. Here we concentrate on a particular class that can be reduced to a linear algebraic problem, with explicit solution. Ill-conditioning of the problem corresponds to the vanishing of one eigenvalue of the matrix to be inverted. The comparison with algorithms based upon matching experimental data against numerical integration of the CDE sheds light on the accuracy of the parameter estimation procedures, and suggests a path for a more precise assessment of the profiles and of the related uncertainty. Several instances of the implementation of the algorithm to real data are presented.
The calculation of transport profiles from experimental measurements belongs in the category of inverse problems which are known to come with issues of ill-conditioning or singularity. A reformulation of the calculation, the matricial approach, is proposed for periodically modulated experiments, within the context of the standard advection-diffusion model where these issues are related to the vanishing of the determinant of a 2x2 matrix. This sheds light on the accuracy of calculations with transport codes, and provides a path for a more precise assessment of the profiles and of the related uncertainty.
We introduce non-trivial contributions to diffusion constant in generic many-body systems arising from quadratic fluctuations of ballistically propagating, i.e. convective, modes. Our result is obtained by expanding the current operator in the vicinity of equilibrium states in terms of powers of local and quasi-local conserved quantities. We show that only the second-order terms in this expansion carry a finite contribution to diffusive spreading. Our formalism implies that whenever there are at least two coupled modes with degenerate group velocities, the system behaves super-diffusively, in accordance with the non-linear fluctuating hydrodynamics theory. Finally, we show that our expression saturates the exact diffusion constants in quantum and classical interacting integrable systems, providing a general framework to derive these expressions.
This paper deals with an inverse problem applied to the field of building physics to experimentally estimate three sorption isotherm coefficients of a wood fiber material. First, the mathematical model, based on convective transport of moisture, the Optimal Experiment Design (OED) and the experimental set-up are presented. Then measurements of relative humidity within the material are carried out, after searching the OED, which is based on the computation of the sensitivity functions and a priori values of the unknown parameters employed in the mathematical model. The OED enables to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs, ensuring the best accuracy for the identification method and, thus, for the estimated parameter. Two experimental procedures were identified: i) single step of relative humidity from 10% to 75% and ii) multiple steps of relative humidity 10-75-33-75% with an 8-day duration period for each step. For both experiment designs, it has been shown that the sensor has to be placed near the impermeable boundary. After the measurements, the parameter estimation problem is solved using an interior point algorithm to minimize the cost function. Several tests are performed for the definition of the cost function, by using the L^2 or L^infty norm and considering the experiments separately or at the same time. It has been found out that the residual between the experimental data and the numerical model is minimized when considering the discrete Euclidean norm and both experiments separately. It means that two parameters are estimated using one experiment while the third parameter is determined with the other experiment. Two cost functions are defined and minimized for this approach. Moreover, the algorithm requires less than 100 computations of the direct model to obtain the solution. In addition, the OED sensitivity functions enable to capture an approximation of the probability distribution function of the estimated parameters. The determined sorption isotherm coefficients calibrate the numerical model to fit better the experimental data. However, some discrepancies still appear since the model does not take into account the hysteresis effects on the sorption capacity. Therefore, the model is improved proposing a second differential equation for the sorption capacity to take into account the hysteresis between the main adsorption and desorption curves. The OED approach is also illustrated for the estimation of five of the coefficients involved in the hysteresis model. To conclude, the prediction of the model with hysteresis are compared with the experimental observations to illustrate the improvement of the prediction.
Given $(M,g)$, a compact connected Riemannian manifold of dimension $d geq 2$, with boundary $partial M$, we consider an initial boundary value problem for a fractional diffusion equation on $(0,T) times M$, $T>0$, with time-fractional Caputo derivative of order $alpha in (0,1) cup (1,2)$. We prove uniqueness in the inverse problem of determining the smooth manifold $(M,g)$ (up to an isometry), and various time-independent smooth coefficients appearing in this equation, from measurements of the solution on a subset of $partial M$ at fixed time. In the flat case where $M$ is a compact subset of $mathbb R^d$, two out the three coefficients $rho$ (weight), $a$ (conductivity) and $q$ (potential) appearing in the equation $rho partial_t^alpha u-textrm{div}(a abla u)+ q u=0$ on $(0,T)times Omega$ are recovered simultaneously.
This paper investigates the identification of two coefficients in a coupled hyperbolic system with an observation on one component of the solution. Based on the the Carleman estimate for coupled wave equations a logarithmic type stability result is obtained by measurement data only in a suitably chosen subdomain under the assumption that the coefficients are given in a neighborhood of some subboundary.