INTEGRATED MULTI-PHYSICS DESIGN TOOL FOR FUSION BREEDING BLANKET SYSTEMS - DEVELOPMENT AND VALIDATION

The Breeding Blanket (BB) of the DEMO reactor represents a harsh system in a dangerous environment. It has to satisfy engineering requirements and constraints that are of nuclear, thermo-structural, material and safety kind. For these reasons, the application of advanced simulation tools, based on a multi-physics approach, is required for its comprehensive design. These tools have to simultaneously perform different kind of analyses among which three, and namely nuclear, thermofluid-dynamic and thermo-mechanical, can be prioritized and considered as propaedeutic for the investigation of all the other issues related to the BB. In this dissertation, a multi-physic approach, covering the three pillars of the BB design (the neutronics, thermal-hydraulics and thermo-mechanics), is proposed. These analyses have to be conducted in a strongly integrated way, allowing a holistic assessment of volumetric heat loads, thermal performances of coolant and structures as well as their stress and deformation states. The strategy, followed for the achievement of this challenge, consists of creating a CAD-centric and loosely-coupled procedure for the BB concepts design adopting a sub-modelling technique, named Multi-physics Approach for Integrated Analysis (MAIA). The MAIA procedure bases its architecture on the use of validated codes and on the minimisation of their number. It is articulated in 10 main steps that go from the decomposition of generic CAD in a format suitable for neutron/photon transport analysis to the nuclear analysis for the assessment of volumetric heating, from the assessment of temperature and velocity fields within coolant and structure to the evaluation of their displacement, deformation and stress fields, from the evaluation of nitrogen isotopes production rates from water oxygen activation to the calculation of their concentration spatial distribution taking into account the effects of passive convective transport. All the steps share the same geometry details and the consistency between input and output parameters. The new MAIA procedure differs from the conventional coupling approach with respect to three key aspects. First, it does not introduce homogenisations of models and loads. Second, MAIA can capture load gradients at high resolution in the three directions for all the analysis involved without requiring prohibitive computational efforts. And third, MAIA keeps the consistency between the three analyses maintaining the congruence between inputs and outputs. However, the computational effort required by the CAD-centric feature of MAIA procedure imposes the representation of BB portions and, therefore, the definition and validation of boundary conditions for each performed calculation. Regarding the nuclear analysis, it has been found that the set of reflecting and white conditions in the poloidal and toroidal directions, respectively, together with the presence of Vacuum Vessel (VV) and the definition of local neutron and photon source, produces a mismatch of -0.48 % in terms of power deposition between the DEMO and the local (e.g. slice) models. It has been demonstrated that the neutronic symmetry conditions are valid in the entire module up to the last slices nearby the caps. Furthermore, a sensitivity analysis on the angular distribution of local neutron and photon source has been performed indicating in 10 cosine bins the optimal discretisation choice in terms of compromise between the fidelity of the results obtained respect to those of the reference model and the relevant computational effort. Concerning the analysis of thermal-hydraulic boundary conditions, it has been found that the variation on mass flow rates (comprised between the ~-1.3 % and the ~0.6 %) as well as power density fluctuation (up to the ~6 % in the neighbouring domains) affect the temperature distribution for less than ±2.4 % demonstrating the applicability of poloidal symmetry conditions. As far as the thermo-mechanical analyses are concerned, it has been identified the set of boundary conditions (radial and toroidal displacements prevented to the nodes lying in the rear of the back supporting structure along the toroidal and poloidal direction, symmetry at the lower cut surface and Generalised Plane Strain to the top one) that produce a discrepancy in terms of displacement in the sub-model comprised between the -6 % and the 4 % as well as a conservative assessment of membrane and bending stresses both for primary and secondary stresses. The impact of the temperature variation has also been investigated showing that the fluctuations on total deformation are comprised between -0.3 % and the 1.7 %, on equivalent membrane stress up to 15 % while on equivalent bending stress between the -7 % and the 5 %. As a proof-of-concept, the MAIA procedure has been then used to evaluate the impact on the BB design, demonstrating that some criticalities are present in the design. In particular, the fluid-dynamic results show a violation of the temperature requirement limits that have not been solved introducing proper design solutions. Furthermore, these violations of thermal-hydraulic requirements produce very intense values of Von Mises equivalent stresses that could jeopardize the structural integrity of the segment box. This demonstrates that MAIA procedure can become the reference tool for the design of the BB. Moreover, the MAIA procedure has proven the possibility to locally map important variables such as the neutron flux and the temperature as well as the primary and secondary stress that are used for the determination of the allowable stress and applied for compering with design criteria. In order to further demonstrate the versatility and adaptability of the MAIA procedure, the water activation issue occurring within the blanket Primary Heat Transfer System (PHTS) has been studied. Using MAIA procedure, it has been possible to take into account the effects of the flow on the nitrogen concentration and to provide useful information for the development of both BB design and its PHTS.

Il Breeding Blanket (BB) del reattore DEMO rappresenta un sistema complesso in un ambiente pericoloso. Infatti, esso deve soddisfare diversi requisiti e vincoli ingegneristici sia di tipo nucleare, termo-strutturale che di sicurezza. Per questi motivi, è necessaria una progettazione omnicomprensiva che preveda l'applicazione di strumenti avanzati di simulazione basati su approcci multi-fisici. Questi strumenti devono eseguire simultaneamente diversi tipi di analisi. Tre di esse possono essere considerate prioritarie e propedeutiche per lo studio di tutti gli altri fenomeni riguardanti il BB, vale a dire l´analisi nucleare, termo-fluidodinamica e strutturale. In questa tesi, è proposto un innovativo approccio multi-fisico che copre i tre pilastri principali su cui è basato il progetto del BB (la neutronica, la termoidraulica e la termo-meccanica). Queste analisi devono essere condotte in maniera integrata, consentendo una valutazione olistica dei carichi volumetrici di potenza, delle prestazioni termiche sia del fluido di raffreddamento che delle strutture, nonché dei campi di tensione e deformazione. La strategia seguita per il conseguimento di questa sfida consiste nella creazione di una procedura “CAD-centric” e “loosely-coupled” (debolmente accoppiata) per la progettazione dei concetti di BB utilizzando una tecnica di analisi basata su sotto-modelli. Questa procedura prende il nome di Multi-physics Approach for Integrated Analysis (MAIA). Essa basa la sua architettura sull'uso di codici validati e sulla minimizzazione del loro numero. In particolare, MAIA è articolata in 10 fasi principali che vanno dalla creazione di un modello per le analisi nucleari generato dalla decomposizione in geometrie semplici di un generico CAD alla valutazione della potenza volumetrica, dal calcolo dei campi di temperatura e velocità nella struttura e nel refrigerante alla valutazione dei campi di spostamento, deformazione e stress, dalla stima dei tassi di produzione degli isotopi dell´azoto prodotti dall'attivazione dell'ossigeno presente nell'acqua al calcolo della loro distribuzione spaziale di concentrazione tenendo conto degli effetti del trasporto convettivo. Tutti i vari passaggi condividono gli stessi dettagli geometrici. In particolare, MAIA differisce dagli approcci convenzionali usati nell´accoppiamento multi-fisico su tre aspetti chiave. Innanzitutto, non introduce omogeneizzazioni dei modelli e dei carichi. In secondo luogo, MAIA permette di determinare, con un’alta risoluzione spaziale, i gradienti dei carichi per tutte le analisi coinvolte senza richiedere sforzi computazionali proibitivi. In terzo luogo, MAIA permette di mantenere la coerenza tra le tre analisi garantendo la congruenza tra gli input e gli output. Tuttavia, l´onere computazionale richiesto dall´approccio CAD-centric, su cui si basa la procedura MAIA, non permette di rappresentare il BB nel suo complesso ma solo alcune sue porzioni (una slice, per esempio). Ciò impone la definizione e, conseguentemente, la validazione di opportune condizioni al contorno per ogni sotto-modello utilizzato e per ogni analisi eseguita. A tal proposito, per quanto riguarda le analisi nucleari, le condizioni al contorno utilizzate nel modello locale della slice sono: definizione di una sorgente locale neutronica/fotonica per tener in conto l´effetto albedo dell´intero reattore, rappresentazione del Vacuum Vessel (VV) per simulare il back scattering verso il BB, e l´applicazione di condizioni di riflessione (“mirror”, specchio/simmetria, nella direzione poloidale e “white”, riflessione isotropica, in quella toroidale) per simulare la presenza delle slice adiacenti a quella analizzata. I risultati ottenuti mostrano una variazione della potenza depositata del -0.48 % tra il modello di riferimento DEMO e quello locale (slice). Inoltre, è stata eseguita un'analisi di sensibilità sulla distribuzione angolare della sorgente neutronica/fotonica locale determinando una discretizzazione ottimale in 10 suddivisioni poloidali. Questa suddivisione rappresenta un buon compromesso sia in termini di fedeltà dei risultati ottenuti, rispetto a quelli del modello di riferimento (DEMO), che di minimizzazione dell´onere computazionale. Per quanto riguarda l'analisi delle condizioni al contorno termo-idrauliche usate nel modello locale della slice, è stata applicata una condizione di simmetria termica poloidale. Assumendo una variazione delle portate comprese tra ~ -1.3% e ~ 0.6% e una fluttuazione della densità di potenza fino a ~ 6% tra slice vicine, è stata ottenuta una variazione della distribuzione delle temperature del ± 2.4% dimostrando, quindi, l'applicabilità di tali condizioni. Per quanto riguarda le analisi termo-meccaniche, le condizioni al contorno identificate per il modello locale della slice sono: simmetria sul piano inferiore della slice, Generalised Plane Strain su quello superiore e spostamenti radiali e toroidali impediti ai nodi che giacciono nella parte posteriore della back supporting structure lungo la direzione toroidale e poloidale. Queste condizioni, applicate al sotto-modello, producono una variazione compresa tra il -6% e il 4% tra gli spostamenti calcolati nella slice e quelli nel modello di riferimento DEMO, nonché una stima conservativa delle tensioni primarie e secondarie sia di membrana che di flessione. Inoltre, è stato anche studiato l'impatto della variazione (± 2.4%) di temperatura dimostrando che le fluttuazioni sulle deformazioni totale sono comprese tra il -0.3% e l’1.7%, fino a un massimo del 15% sulle tensioni equivalenti di membrana e tra il -7% e il 5% su quelle di flessione. Infine, la procedura MAIA è stata utilizzata per valutare l'impatto sul design del BB. La sua applicazione ha dimostrato la presenza di alcune criticità nel progetto. In particolare, i risultati fluidodinamici mostrano una violazione dei limiti di temperatura che non sono stati risolti introducendo soluzioni progettuali adeguate. Inoltre, queste violazioni producono, a loro volta, valori molto intensi delle tensioni equivalenti di Von Mises che potrebbero indicare un pericolo per l'integrità strutturale del BB. L´applicazione di MAIA al design del BB a permesso di dimostrare il valore aggiunto di questa procedura la quale potrebbe diventare uno strumento fondamentale e di riferimento per la progettazione del BB. Inoltre, la procedura MAIA ha permesso di mappare localmente variabili importanti come flussi neutronici e temperature, nonché le tensioni primarie e secondarie che sono utilizzate per la determinazione delle tensioni ammissibili applicate per la verifica dei criteri di progettazione. Al fine di dimostrare ulteriormente la versatilità e l'adattabilità della procedura MAIA, è stato studiato il problema di attivazione dell'acqua del sistema di trasferimento di calore primario (Primary Heat Transfer System, PHTS). Utilizzando la procedura MAIA, è stato possibile prendere in considerazione gli effetti dell´efflusso sulla concentrazione degli isotopi dell´azoto e fornire informazioni utili per lo sviluppo sia del design del BB che del suo PHTS.