Explicit Discrete Solution for Some Optimization Problems and Estimations with Respect to the Exact Solution
This work addresses incremental improvements in numerical methods for heat conduction optimization, relevant for computational physics and engineering applications.
The paper tackles optimization problems for steady-state heat conduction systems by deriving explicit discrete solutions using finite difference schemes, proving convergence and error estimates as discretization parameters approach limits, and showing that a three-point approximation improves convergence order from O(h) to O(h^2).
We consider two steady-state heat conduction systems called, $S$ and $S_α$, in a multidimensional bounded domain $D$ for the Poisson equation with source energy $g$. In one system, we impose mixed boundary conditions (temperature $b$ on the boundary $Î_1$, heat flux $q$ on $Î_2$ and an adiabatic condition on $Î_3$). In the other system, the condition on $Î_1$ is replaced by a convective heat flux condition with coefficient $α$. For each of these systems, we consider three associated optimization problems $(P_{i})$ and $(P_{iα})$, $i=1,2,3$, where the variable is the source energy $g$, the heat flux $q$ and the environmental temperature $b$, respectively. In the particular case where $D$ is a rectangle, the explicit continuous optimization variables and the corresponding state of the systems are known. In the present work, by using a finite difference scheme, we obtain the discrete systems $({S^h})$ and ${(S^h_α)}$ and discrete optimization problems ${(P^h_i)}$ and ${(P^h_{i α})}$, $i=1,2,3$, where $h$ is the space step in the discretization. Explicit discrete solutions are found, and convergence and estimation errors results are proved when $h$ goes to zero and when $α$ goes to infinity. Moreover, some numerical simulations are provided in order to test theoretical results. Finally, we note that the use of a three-point finite-difference approximation for the Neumann or Robin boundary condition at the boundary improves the global order of convergence from $O(h)$ to $O(h^2)$.