This paper presents a numerical solution to the inverse problem consisting of recovering time-dependent thermal conductivity and heat source coefficients in the one-dimensional parabolic heat equation. This mathematical formulation ensures that the inverse problem has a unique solution. However, the problem is still ill-posed since small errors in the input data lead to a drastic amount of errors in the output coefficients. The finite difference method with the Crank-Nicolson scheme is adopted as a direct solver of the problem in a fixed domain. The inverse problem is solved subjected to both exact and noisy measurements by using the MATLAB optimization toolbox routine lsqnonlin , which is also applied to minimize the nonlinear Tikhonov regularization functional. The thermal conductivity and heat source coefficients are reconstructed using heat flux measurements. The root mean squares error is used to assess the accuracy of the approximate solutions of the problem. A couple of numerical examples are presented to verify the accuracy and stability of the solutions.
In this article, an inverse problem of finding timewise-dependent thermal conductivity has been investigated numerically. Numerical solution of forward (direct) problem has been solved by finite-difference method (FDM). Whilst, the inverse (indirect) problem solved iteratively using Lsqnonlin routine from MATLAB. Initial guess for unknown coefficient expressed by explicit relation based on nonlocal overdetermination conditions and intial input data .The obtained numrical results are presented and discussed in several figures and tables. These results are accurate and stable even in the presense of noisy data.
This paper examines the finding of spacewise dependent heat source function in pseudoparabolic equation with initial and homogeneous Dirichlet boundary conditions, as well as the final time value / integral specification as additional conditions that ensure the uniqueness solvability of the inverse problem. However, the problem remains ill-posed because tiny perturbations in input data cause huge errors in outputs. Thus, we employ Tikhonov’s regularization method to restore this instability. In order to choose the best regularization parameter, we employ L-curve method. On the other hand, the direct (forward) problem is solved by a finite difference scheme while the inverse one is reformulated as an optimization problem. The
... Show MoreExperimental study of heat transfer coefficients in air-liquid-solid fluidized beds were carried out by measuring the heat rate and the overall temperature differences across the heater at different operating conditions. The experiments were carried out in Q.V.F. glass column of 0.22 m inside diameter and 2.25 m height with an axially mounted cylindrical heater of 0.0367 m diameter and 0.5 m height. The fluidizing media were water as a continuous phase and air as a dispersed phase. Low density (Ploymethyl-methacrylate, 3.17 mm size) and high density (Glass beads, 2.31 mm size) particles were used as solid phase. The bed temperature profiles were measured axially and radially in the bed for different positions. Thermocouples were connecte
... Show MorePhase change materials (PCMs) such as paraffin wax can be used to store or release large amount of energy at certain temperature at which their solid-liquid phase changes occurs. Paraffin wax that used in latent heat thermal energy storage (LHTES) has low thermal conductivity. In this study, the thermal conductivity of paraffin wax has been enhanced by adding different mass concentration (1wt.%, 3wt.%, 5wt.%) of (TiO2) nano-particles with about (10nm) diameter. It is found that the phase change temperature varies with adding (TiO2) nanoparticles in to the paraffin wax. The thermal conductivity of the composites is found to decrease with increasing temperature. The increase in thermal conductivity ha
... Show MoreThis paper investigates the simultaneous recovery for two time-dependent coefficients for heat equation under Neumann boundary condition. This problem is considered under extra conditions of nonlocal type. The main issue with this problem is the solution unstable to small contamination of noise in the input data. The Crank-Nicolson finite difference method is utilized to solve the direct problem whilst the inverse problem is viewed as nonlinear optimization problem. The later problem is solved numerically using optimization toolbox from MATLAB. We found that the numerical results are accurate and stable.