09.07.2015 Views

Couplages multi-physiques et applications aux ... - gedepeon

Couplages multi-physiques et applications aux ... - gedepeon

Couplages multi-physiques et applications aux ... - gedepeon

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

POURQUOI UN COUPLAGE DES CODES? (3/3) Pour le fonctionnement normal du réacteur, le couplage n/th est effectué via lessections efficaces qui sont réévaluées pour chaque maille spatiale <strong>et</strong> en fonction duchangement des propriétés combustible <strong>et</strong> modérateur ; celles-ci sont en généralévaluées via un modèle stationnaire intégré dans le code de cœur.Atelier GEDEPEON 11-12 juill<strong>et</strong> 20114


L’OUTIL HEMERAOUTIL D’ANALYSE METHODOLOGIQUE ET DE CONCEPTIONEN PARTIE CO-DEVELOPPE AVEC L’IRSNDesign cœur : REP (900, 1300, N4, EPR)• Souplesse sur la gestion du combustible <strong>et</strong> le pilotage du cœurScénarios de fonctionnement: normal à accidentelAccidents traités: RTV, REAscenario accidentel1 2 3NominaloperatingconditionscalculationAccidentinitiatedconditionscalculationTransientcalculationAtelier GEDEPEON 11-12 juill<strong>et</strong> 20118


SCENARIOS ACCIDENTELSAccident REAAccident RTVInitiateur: rupture cartergrappe de contrôle⇒Ejection grappe en 0.1s (force d’envol):Libération réactivité (ρ g) = poids neutroniquede la grappe (fonction de son enfoncement)1.Puissance cœur ↑ - zone chaude auvoisinage de la grappe éjectée2.Puissance cœur ↓ par contre-réactionsDoppler (α D- 1 er ordre) + modérateur(α m-2 ème ordre )3.Puissance cœur ↓ par chute grappes arrêtd’urgence (ρ arr<strong>et</strong>)ρ(t) =" ρ −αDΔTcomb−αmΔTg 43 mod14243142C.R DopplerC.R moderateurDurée du transitoire: qq secondesBesoins en performances: pas de temps1 ms sur 2s (2000 steps de calculs),environ 24 heures CPU-ρP(x, y,z, t) =arr<strong>et</strong>"P (t)123comportement intégralInitiateur: rupture guillotine instantanée<strong>et</strong> doublement débattue de la tuyauterievapeur en sortie d’un générateur devapeur <strong>et</strong> en amont de la vanned’isolement vapeurLe cœur est initialement sous-critique (ρ 0)Dépressurisation secondaire1. Transfert thermique GV ↑ Température caloporteur ↓3. Puissance cœur ↑ par eff<strong>et</strong> modérateur (α m)4. Ralentissement puissance cœur par eff<strong>et</strong>Doppler (α D)5. Puissance cœur ↓ par Injection boreρ(t) = " ρ0 + αmΔTmod−αDΔTcomb−ρ1424314243C.R moderateurDurée du transitoire: qq 100 secondesBesoins en performances: pas de temps20s (2000 steps de calculs), environ 3 joursCPUcoeur⋅ t)F(x, y,z,14243comportement localAtelier GEDEPEON 11-12 juill<strong>et</strong> 2011C.R Dopplerbore"9


SCENARIOS ACCIDENTELSAccident REA⇒Transitoire violent (ρ>β)(temps double puissance: qq ms)Accident RTV⇒ Transitoire "lent" (ρ autour de β)(temps double puissance: ∼100 ms)654Pcoeur32100 50 100 150 200 250 300Temps [s]Atelier GEDEPEON 11-12 juill<strong>et</strong> 201110


Que signifie HEMERA ?“Highly Evolutionary M<strong>et</strong>hods for Extensive Reactor Analyses”mais aussiDivinité grecque dont le rôle est de « dissiper les ténèbreschaque matin <strong>et</strong> de révéler la lumière divine <strong>aux</strong> hommes »Qu’est-ce que l’outil HEMERA ?3 codes couplés : CATHARE2, FLICA4, CRONOS2 + 1 superviseur• CATHARE2 simule la thermohydraulique des circuitsprimaire <strong>et</strong> secondaire en 1D à l’exception du cœur• FLICA4 simule la thermohydraulique du cœur en 3D• CRONOS2 simule la neutronique du cœur en 3DAtelier GEDEPEON 11-12 juill<strong>et</strong> 201111


CONTENU D’ HEMERACode CRONOS2– Code neutronique 3D modulaire– Résolution de l’équation de ladiffusion <strong>et</strong> du transport– Calculs statiques <strong>et</strong> cinétiques– Evolution isotopique <strong>et</strong> du burn-up– Alimenté en amont en sectionsefficaces par le code APOLLO2APOLLO2 <strong>et</strong> CRONOS2 ouverts à l’Europe(via NURISP) <strong>et</strong> <strong>aux</strong> universités via un accordbilatéral (accès à la documentation, <strong>aux</strong>sessions formation)Atelier GEDEPEON 11-12 juill<strong>et</strong> 201112


CONTENU D’ HEMERACode FLICA4• Code de thermohydraulique 3D diphasique• Modèle à 4 équations : 3 pour le mélange (masse, moment<strong>et</strong> énergie) <strong>et</strong> 1 pour la conservation de la masse vapeur• 2 lois de ferm<strong>et</strong>ure (drift + vapeur à saturation)• Module thermique 1D• FLICA4 ouvert à l’Europe(viaNURISP) <strong>et</strong> <strong>aux</strong> universités(accès à la documentation <strong>et</strong><strong>aux</strong> sessions de formation)Atelier GEDEPEON 11-12 juill<strong>et</strong> 201113


CONTENU D’ HEMERACode CATHARE2• Code système “best-estimate” développé depuisplus de 30 ans dans le cadre d’une collaborationCEA/EDF/AREVA/IRSN• Modèle à 6 équations <strong>et</strong> 2 phases• Modules 0D, 1D <strong>et</strong> 3D• Possibilité d’un accord delicence validé par lespartenaires industrielsAtelier GEDEPEON 11-12 juill<strong>et</strong> 201114


CONTENU D’ HEMERAHEMERA V2: la plate-forme SALOMEPlate-forme de de développement générique‣ Pre-processing‣ Post-processing‣ Modélisations <strong>multi</strong><strong>physiques</strong>• CATHARE2• FLICA4• CRONOS2Atelier GEDEPEON 11-12 juill<strong>et</strong> 201115


LA MODELISATIONAPOLLO2Sections efficaces Vrai couplage (utilisation de la sortie de FLICA en entrée de CATHARE) Pas de recouvrement de domaine (thermohydraulique du cœur simuléeseulement par FLICA en 3D)Atelier GEDEPEON 11-12 juill<strong>et</strong> 201116


ACCES AUX GRANDEURS D’INTERËTAnalyse de l’état neutronique du cœur─ IntégraleLes grandeurs d’intérêt de la modélisation neutroniqueperm<strong>et</strong>tent une analyse intégrale <strong>et</strong> locale du comportement du cœuren situations de fonctionnement stationnaires <strong>et</strong> transitoires- Réactivité totale du cœur <strong>et</strong> par poste (grappes, Doppler, modérateur,<strong>et</strong> bore)- Facteurs de forme, F q , F 3D , F xy , F z, Axial Offs<strong>et</strong> (AO)─ Locale- Distributions de BU 3D, radiale (2D) <strong>et</strong> axiale (1D) cœur, par lotd’assemblages, par type de combustible- Distributions de puissance 3D, radiale (2D) <strong>et</strong> axiale (1D) cœur <strong>et</strong>assemblage- Accès à la puissance locale par méthode de reconstruction fine depuissanceAnalyse de l’état thermique <strong>et</strong> thermohydraulique du cœur/système─ Large gamme de sorties disponibles dans FLICA4 <strong>et</strong> CATHAREAtelier GEDEPEON 11-12 juill<strong>et</strong> 201117


EXEMPLE DE MODELISATIONCATHARE : circuitssecondaire <strong>et</strong> primairesauf le coeurCRONOS <strong>et</strong> FLICA : Coeur 3Dau niveau assemblageCRONOS <strong>et</strong> FLICA : Assemblagechaud au niveau crayonGrappe bloquéeBoucle affectéeAtelier GEDEPEON 11-12 juill<strong>et</strong> 201118


APPLICATIONSExemple 1 : la Rupture de Tuyauterie Vapeur (RTV)Rupture1- brèche sur une tuyauterie vapeur2- refroidissement du primaire3- insertion de réactivité4- r<strong>et</strong>our potentiel en puissance avecrisque d’intégrité de la premièrebarrière Phénomènes thermohydrauliques<strong>et</strong> neutroniques fortement couplés<strong>et</strong> tridimensionnelsAtelier GEDEPEON 11-12 juill<strong>et</strong> 201119


APPLICATIONSExemple 1 : la Rupture de Tuyauterie Vapeur (RTV)654Pcoeu32100 50 100 150 200 250 300Temps [s]Atelier GEDEPEON 11-12 juill<strong>et</strong> 201120


APPLICATIONSExemple 2 : L’Ejection D’une Grappe (EDG)Ejection ⇒ montée en puissance T comb ⇒ Contre-réaction Doppler T mod ⇒ Contre-réaction modérateur Arrêt d’urgence (chute des grappes de commande)⇒Transitoire rapide de puissance (pulse)⇒ Dépôt d’énergie fortement localisé àproximité de la grappe éjectéeModélisation cinétique tridimensionnel<strong>et</strong> <strong>multi</strong> physiqueAtelier GEDEPEON 11-12 juill<strong>et</strong> 201121


LES ANALYSES DE SENSIBILITELes paramètres clés de la modélisation neutronique sont les paramètres quiconditionnent le comportement du cœur en conditions de fonctionnementnormales <strong>et</strong> accidentelles1.Les données de base de la modélisationSections efficacesParamètres cinétiques ( fraction de neutrons r<strong>et</strong>ardés)2.Les paramètres de définition du point de fonctionnementT<strong>aux</strong> de combustionTempérature combustibleDensité/température modérateurConcentration en boreLes leviers mis en œuvre dans HEMERA perm<strong>et</strong>tent une action sur cesparamètres:En espace:Modification globale sur tout le cœurModification locale dans un groupe d’assemblages, un assemblageune maille précise du cœur, …Atelier GEDEPEON 11-12 juill<strong>et</strong> 201123


CONCLUSIONS ET PERSPECTIVES HEMERA est un outil pertinent en regard des objectifs fixés :études en support à l’expertise, montée en compétence,compréhension des phénomènes, analyse de concepts … LES PERSPECTIVESExtension périmètre d’utilisationConsolidation architecture logicielle pour <strong>applications</strong> à d’autres réacteurs,passage à APOLLO3 (outil <strong>multi</strong>-filières) en remplacementd’APOLLO2/CRONOS2, introduction nouve<strong>aux</strong> scénarios accidentelsAnalyse d’incertitudesMise en place d’une méthodologie d’analyse d’incertitudes couvrantl’ensemble des paramètres sensibles (neutronique – thermohydraulique –thermique) dans le cadre des accidents RIA <strong>et</strong> RTVModélisation thermique du combustibleAmélioration du calcul des températures dans le crayon au regard de sonévolution en cours d’irradiation <strong>et</strong> en situation accidentelle (couplage avecl’outil CEA/DEN PLEIADES)Supervision (SALOME)Souplesse dans le contrôle des simulations <strong>et</strong> la mise en œuvre de moyensd’analyse (plan d’expériences, analyse d’incertitudes)Atelier GEDEPEON 11-12 juill<strong>et</strong> 201124

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!