Моделирование осесимметричной теплопроводности в компактных изделиях керамического ядерного топлива с учетом температурных зависимостей теплофизических характеристик

Abstract

Обсуждаются возможности применения различных математических формулировок для моделирования осесимметрично теплопроводности компактных изделий керамического ядерного топлива. Показано, что применение уравнения теплопроводности, записанного относительно температуры, может вносить погрешности, связанные с неопределенностью исходных данных о производной температурной зависимости коэффициента теплопроводности, которая доступна исключительно в табулированной форме для значений температуры, разделенных достаточно большим шагом. Это обстоятельство является существенным для моделирования осесимметричной теплопроводности компактных изделий керамического ядерного топлива, поскольку его теплопроводность существенно зависит от температуры, уменьшаясь в 2,5 раза при изменении температуры от 323 K до 1073 K. Показано, что для изучения осесимметричной теплопроводности компактных изделий керамического ядерного топлива наибольший интерес представляет смешанная математическая формулировка задачи теплопроводности относительно полей температуры и вектора теплового потока, поскольку соответствующие дифференциальные уравнения не содержат производной температурной зависимости коэффициента теплопроводности. При этом погрешности, вносимые аппроксимацией значений коэффициента теплопроводности по имеющимся табличным данным, будут ограничены погрешностью аппроксимации значений коэффициента теплопроводности, которая легко контролируется по имеющимся табличным данным о значениях коэффициента теплопроводности. Для решения задачи теплопроводности, формулированной в смешанной форме относительно полей температуры и вектора теплового потока, предлагается использовать метод полу-дискретизации, который сводит рассматриваемую задачу к определению временных зависимостей искомых величин в отдельных точках исследуемой области компактного изделия керамического ядерного топлива. Для этого предлагается конечно-разностными формулами заменять производные только по пространственным координатам, что позволит получить обыкновенные дифференциальные уравнения с начальными условиями для определения узловых значений искомых величин).

The possibilities of using of various mathematical formulations for modelling the axisymmetric heat conduction of compact products of ceramic nuclear fuel are discussed. It is shown that the application of the heat equation wrote with respect to temperature can introduce errors due to the uncertainty of the initial data on the derivative of the temperature dependence of the thermal conductivity, which is available exclusively in a tabulated form for temperature values separated by a sufficiently large step. This circumstance is essential for modelling the axisymme tric heat conduction of compact products of ceramic nuclear fuel, since its thermal conductivity depends significantly on temperature, decreasing by a factor of 2.5 with a temperature change from 323 K to 1073 K. It is shown that to study of the axisymmetric heat conduction of compact products of ceramic nuclear fuel, the most interesting is the mixed mathematical formulation of the heat conduction problem with respect to the temperature fields and the heat flux vector, since the corresponding differential equations do not contain the derivative of the temperature dependence of the thermal conductivity. In this case, the approximations of the values the thermal conductivity coefficient from the available tabular data will be limited by the error in approximating the values of the thermal conductivity, which is easily monitored from the available tabular data on the values of the thermal conductivity coefficient. To solve the heat conduction problem, formulated in mixed form with respect to the temperature fields and the heat flux vector, it is proposed to use the semi-discretisation method, which reduces the problem under consideration to determining the time dependences of the unknown quantities at individual points of the investigated region of a compact product of ceramic nuclear fuel. To this end, it is proposed to replace finite-difference formulas by finite-difference formulas only in terms of spatial coordinates, which will allow us to obtain ordinary differential equations with initial conditions for determining the nodal values of the unknown quantities.