Download rayo-axi-3.dgibi
Back to the list
* fichier : rayo-axi-3.dgibi
*************************************************************
* *
* Calcul d'un cylindre infini (de rayon L) soumis à de *
* la convection et du rayonnement. *
* *
* Modélisation axisymétrique. *
* *
* Auteurs : Michel Bulik & Nadia Coulon *
* *
* Date : Octobre 1995 *
* *
* Références : *
* ------------ *
* [1] Klaus-Jürgen Bathe & Mohammad R. Khoshgoftaar, *
* Finite element formulation and solution of *
* non-linear heat transfer, Nuclear Engineering and *
* Design, v. 51 (1979), pp. 389-401 *
* *
* [2] J. Joly, Cas tests non linéaires de validation pour *
* DELFINE, Note technique EMT.SMTS.TTMF 84/29 *
* *
* [3] James SUCEC & Ashok KUMAR, Transient Cooling of a *
* Solid Cylinder by Combined Convection ans Radiation *
* at its Surface, International Journal for Numerical *
* Methods in Engineering, vol. 6 (1973), pp. 297-312 *
* *
*************************************************************
*
* La solution de référence a été obtenue par une méthode de
* différences finies décrite dans [3] et que réalise le programme
* suivant en C. Après sa compilation il convient de demander
* l'éditions de liens avec la subroutine gauscl de Castem2000.
*
*
* #include
* #include
*
* #define N 26 /* Nombre de points */
*
* double A[N*N], phi_n[N], phi_n1[N] ;
*
* void main()
* {
* double L = 1. ;
* double eta = 0.55 ;
* double Bi = 1.0 ;
* double Gamma = 1.0 ;
* double alpha = 1.0 ;
*
* double delta_r, delta_R, delta_t, M ;
* int i, j, k, dimension, iprint ;
*
* delta_t = 0.001 ;
* delta_r = L / (N - 1) ;
* delta_R = delta_r / L ;
* M = alpha * delta_t / (delta_r*delta_r) ;
* dimension = N ;
* iprint = 0 ;
*
* /* État initial */
* for(i=0;i maxerr) ;
maxerr = errrela ;
finsi ;
valref = extr listcent &surtst ;
valcalc = extr tcalc T p1 ;
valcalcr = (valcalc - Tempext) / (Tempini - Tempext) ;
errrel = valcalcr - valref ;
mess 'Erreur au centre à t =' instant 'est de' errrel ;
errrela = abs errrel ;
si(errrela > maxerr) ;
maxerr = errrela ;
finsi ;
fin surtst ;
si(maxerr > 0.02) ;
erre 5 ;
finsi ;
*** Bye ...
fin ;