THN

Resolution d un probleme simplfie et couple de Thermo-Hydrauliqte Neutronique
comparaison d'une approche numerique et de solutions analytiques qui utilisent une fonction elliptique complete du premier ordre
le modele numerique et la resolution analytique sont decrits en detail dans l article d'Olivier Lafitte et de l'auteur
    'Analytic solutions and numerical method for a coupled thermo-neutronic problem'
    https://arxiv.org/abs/2205.05945  et https://hal.science/hal-03663985
on se donne une fonction u -> Sigma(u) definie sur [0,1] a partir de trois valeurs  Sigma(0) = v00, Sigma(1/2) = v12,  Sigma(1) = v11
on cherche deux fonctions phi et u definies sur l intervalle [0,1]  et un nombre lambda > 0 appele zlambda dans le programme  et strictement positif de sorte que
   − d2phi / dx**2 + phi(x) = lambda Sigma(u(x)) phi(x)     pour 0 < x < 1
    dh / dx = phi(x),    pour 0 < x < 1
    u(0) = 0, u(1) = 1, phi(0) = phi(1) = 0   et  phi(x) > 0 pour  0 < x < 1. 
On doit construire la fonction  Sigma(u) par interpolation

Ce logiciel ne traite que quelques cas particuliers parametres par deux nombres entiers
  flag         : entre 0 et 12      choix des trois valeurs v00, v12, v11 
  interpol_ini : entre 0 et 5       type d'interpolation

Ce travail a ete effectue par l'auteur alors qu'il etait en delegation CNRS à l'International Research Laboratory au Centre de Recherches Mathematiques de Montreal de février 2021 a juillet 2021.