Thermodynamically consistent phase field model for hydrogen-assisted cracking

Questo articolo presenta un modello a campo di fase termodinamicamente consistente che simula la fessurazione assistita dall'idrogeno in materiali policristallini accoppiando la propagazione delle crepe con la segregazione dell'idrogeno e la riduzione dell'energia interfacciale, catturando con successo la transizione dal cedimento transgranulare a quello intergranulare sotto meccanismi di decoesione potenziata dall'idrogeno.

L'essenziale

Immaginate di avere una struttura metallica molto resistente, come il telaio di un'auto o un ponte. Vi aspettereste che regga sotto pressione, ma a volte atomi di idrogeno invisibili si intrufolano all'interno del metallo e causano la sua rottura improvvisa. Questo fenomeno è chiamato infragilimento da idrogeno. È come se il metallo venisse segretamente "avvelenato" dall'interno, diventando fragile e incline a spezzarsi.

Gli scienziati hanno cercato di costruire modelli informatici per prevedere esattamente come e dove questo metallo si romperà. Tuttavia, i modelli precedenti avevano un difetto fondamentale: trattavano il comportamento dell'idrogeno come una regola semplice e uniforme che si applica ovunque, anche se la struttura interna del metallo è in realtà un complesso mosaico di diversi grani e bordi.

Il Nuovo Modello "Intelligente"
Gli autori di questo articolo hanno creato una nuova, più sofisticata simulazione al computer (chiamata "modello a campo di fase") che agisce come una mappa ad alta definizione e termodinamicamente coerente. Ecco come funziona, utilizzando alcune analogie quotidiane:

• Il Metallo come una Folla: Immaginate che il metallo sia una stanza affollata piena di persone (gli atomi di metallo). I "bordi dei grani" sono le linee invisibili che separano i diversi gruppi di persone. La "crepa" è un vuoto che si allarga nella folla.
• L'Idrogeno come un Ospite Appiccicoso: Gli atomi di idrogeno sono come ospiti appiccicosi che amano nascondersi negli spazi vuoti tra le persone. Hanno una preferenza speciale: amano attaccarsi ai bordi della crepa e alle linee tra i gruppi (bordi dei grani) ancora più di quanto gli piaccia stare nel mezzo della folla.
• Il Problema della "Colla": In un metallo sano, la "colla" che tiene uniti i bordi della crepa è forte. Ma quando questi ospiti appiccicosi di idrogeno si radunano ai bordi della crepa, agiscono come un olio scivoloso, indebolendo la colla. Questo rende la crepa molto più facile da aprire.
• Il Vecchio vs il Nuovo Approccio:
  • I Vecchi Modelli: Usavano un libro di regole generico (l'isoterma di Langmuir-McLean) che assumeva l'idrogeno distribuito uniformemente e in perfetto equilibrio ovunque. Questo è come assumere che tutti nella stanza affollata siano fermi e distribuiti uniformemente, il che non è vero quando si sta formando una crepa.
  • Il Nuovo Modello: Utilizza un quadro "variazionale" flessibile (basato sul formalismo di Kim-Kim-Suzuki). Invece di imporre una regola rigida, permette all'idrogeno di "migrare" naturalmente dove vuole andare (ai bordi della crepa e ai bordi dei grani) in base alle condizioni locali. Calcola esattamente quanto la "colla" si indebolisce in tempo reale man mano che l'idrogeno si raduna.

Cosa hanno Scoperto
Il team ha testato il loro nuovo modello con due scenari principali:

1. Il Test della Crepa Singola: Hanno simulato una crepa in un singolo pezzo di metallo. Senza idrogeno, la crepa cresceva esattamente come prevede la fisica (seguendo il criterio di Griffith). Quando hanno aggiunto l'idrogeno, il modello ha mostrato che la crepa cresceva molto più facilmente perché l'idrogeno aveva indebolito l'energia superficiale. I risultati corrispondevano perfettamente alle previsioni teoriche, dimostrando che il modello funziona.

2. Il Test Policristallino (La Grande Scoperta): Hanno simulato un metallo composto da molti piccoli cristalli (grani) con bordi tra di essi.

   • Senza Idrogeno: La crepa preferiva travolgere direttamente i grani (frattura transgranulare). Era come una palla da demolizione che abbatte i muri di una casa perché i muri sono più deboli della malta tra di essi.
   • Con l'Idrogeno: L'idrogeno si è radunato pesantemente ai bordi tra i grani, indebolendo la "malta" significativamente più delle "pareti". Improvvisamente, la crepa ha cambiato percorso. Invece di travolgere i grani, ha iniziato a serpeggiare lungo i bordi (frattura intergranulare). È stato come se l'idrogeno avesse trasformato la malta in sabbia bagnata, facendo crollare la casa lungo le giunture piuttosto che attraverso i mattoni.

Perché questo è Importante
Questo nuovo modello è un importante aggiornamento perché non si limita a indovinare dove va l'idrogeno; lo calcola in base alla termodinamica effettiva del sistema. Cattura con successo la transizione da un tipo di frattura all'altro, il che è fondamentale per comprendere perché i materiali falliscono in presenza di idrogeno.

Gli autori osservano che, sebbene questo modello sia un grande passo avanti, attualmente si concentra su un meccanismo specifico (l'indebolimento della colla da parte dell'idrogeno). Il lavoro futuro dovrà aggiungere altri fattori complessi, come il modo in cui il metallo si piega e si torce (plasticità) e come altri tipi di difetti interagiscono con l'idrogeno. Ma per ora, questo modello fornisce un modo chiaro, coerente e accurato per vedere come l'idrogeno trasformi un metallo forte in uno fragile.

Sommario tecnico

Dichiarazione del Problema
L'infragilimento da idrogeno (HE) rappresenta una sfida critica per le tecnologie dell'idrogeno, causando una degradazione catastrofica dei materiali e fratture impreviste nei sistemi energetici decarbonizzati. Sebbene siano stati proposti meccanismi quali la decoesione indotta dall'idrogeno (HEDE), la plasticità localizzata indotta dall'idrogeno (HELP) e le vacanze indotte dallo sforzo indotte dall'idrogeno (HESIV), la previsione dell'HE rimane difficile a causa della sua natura multifisica che spazia dalla scala atomica a quella macroscopica. I modelli esistenti di phase field (PF) per la frattura assistita dall'idrogeno tipicamente accoppiano la meccanica della frattura fragile con il meccanismo HEDE tramite la diffusione dell'idrogeno. Tuttavia, questi modelli soffrono di due limitazioni primarie: non riproducono esplicitamente la segregazione dell'idrogeno (il fenomeno responsabile della riduzione dell'energia superficiale) e si affidano all'isoterma di Langmuir-McLean per calcolare la copertura di idrogeno. Quest'ultima è una relazione di equilibrio macroscopica tra fasi bulk e superficiali omogenee, il che la rende fisicamente ingiustificata come legge locale per la modellazione di microstrutture eterogenee che coinvolgono interazioni elastiche.

Metodologia
Gli autori propongono un modello phase field termodinamicamente consistente basato sul formalismo Kim-Kim-Suzuki (KKS) per simulare la frattura assistita dall'idrogeno in materiali policristallini. Il modello opera all'interno di un quadro variazionale in cui l'energia libera totale include i contributi dell'energia superficiale della cricca, dell'energia elastica, dell'accoppiamento dei bordi di grano (GB) e dell'energia chimica.

Le caratteristiche metodologiche chiave includono:

• Variabili di Campo: Il sistema è descritto da un campo di danno ($\eta$), una probabilità di occupazione dell'idrogeno ($c$) e un insieme di parametri d'ordine ($\phi_i$) che rappresentano le orientazioni cristallografiche.
• Formalismo KKS: Invece di imporre una legge locale di Langmuir-Mclean, il modello definisce fasi distinte per il bulk ($c_b$) e la superficie della cricca ($c_{ck}$). La densità di energia libera chimica è costruita come una somma ponderata dei contributi bulk e della cricca, vincolata dall'uguaglianza dei potenziali chimici ($\mu_b = \mu_{ck}$). Ciò consente di disaccoppiare le proprietà interfacciali dai campi di concentrazione.
• Consistenza Termodinamica: Il modello assume soluzioni solide ideali per l'energia chimica e deriva le equazioni di equilibrio governanti (equilibrio meccanico e bilancio del potenziale chimico) dal funzionale dell'energia libera totale.
• Implementazione Numerica: Le equazioni di equilibrio sono risolte utilizzando metodi Fourier-spettrali e un solutore FFT discreto per l'equilibrio meccanico.

Contributi Chiave

1. Modellazione Esplicita della Segregazione: Il modello tiene conto naturalmente della segregazione dell'idrogeno sulle superfici delle cricche e sui bordi di grano senza assumere un isomero macroscopico come legge locale.
2. Derivazione Variazionale della Degradazione dell'Energia Superficiale: La degradazione effettiva dell'energia superficiale della cricca indotta dalla segregazione dell'idrogeno è calcolata analiticamente all'interno del quadro termodinamico.
3. Cattura della Transizione: Il modello dimostra la capacità di catturare la transizione dalla frattura transgranulare alla frattura intergranulare assistita dall'idrogeno tramite il meccanismo HEDE.

Risultati

• Validazione in Monocristalli: Le simulazioni numeriche di una cricca preesistente in un monocristallo sotto carico di Modo I hanno mostrato che lo sforzo critico per la crescita della cricca ($\sigma_c$) segue il criterio di Griffith. Quando il meccanismo HEDE è stato attivato, $\sigma_c$ è diminuito a causa di una riduzione dell'energia superficiale effettiva ($\Delta \gamma_{ch} \approx -4,29$ J/m$^2$), in linea con le previsioni analitiche derivate dalle equazioni del modello.
• Applicazione Policristallina: Le simulazioni in un materiale policristallino a base di nichel hanno rivelato una distinta transizione nei percorsi di propagazione della cricca. In un ambiente privo di idrogeno, la cricca si è propagata in modo transgranulare a causa dell'elevata tenacità dei bordi di grano (GB). Con l'introduzione dell'idrogeno, il modello ha previsto uno spostamento verso la frattura intergranulare. Ciò è stato attribuito a una drastica riduzione del rapporto tra l'energia di frattura GB-bulk effettiva (da 0,5 a 0,07), causata da energie di segregazione dell'idrogeno più elevate sulle superfici di tipo GB rispetto alle superfici di tipo bulk. La mappa della concentrazione di idrogeno ha confermato una significativa segregazione ai bordi di grano e alle superfici di bordo di grano fratturate.

Significatività
L'articolo sostiene che questo lavoro fornisce un metodo rigoroso e generico per modellare il meccanismo HEDE alla scala microstrutturale. Utilizzando il formalismo KKS, il modello riproduce la segregazione dell'idrogeno e la conseguente riduzione dell'energia superficiale senza fare affidamento sull'isoterma di Langmuir-McLean, termodinamicamente limitata per situazioni eterogenee. Lo studio convalida che l'isoterma emerge naturalmente solo in specifici limiti omogenei e statici, giustificando così il nuovo approccio per microstrutture complesse. Il modello cattura con successo la transizione osservata sperimentalmente dalla frattura transgranulare a quella intergranulare utilizzando dati di input realistici per il nichel, offrendo uno strumento robusto per comprendere il cedimento assistito dall'idrogeno nei metalli policristallini. Gli autori osservano che il lavoro futuro affronterà l'anisotropia dell'energia superficiale e l'incorporazione di altri meccanismi di HE, come il modello HELP, attraverso la modellazione delle dislocazioni.