Atomistic-Continuum Coupling by Homogenization

Sintesi Tecnica: Accoppiamento Atomistico-Continuo tramite Omogeneizzazione

Enunciato del Problema
Le simulazioni atomistiche classiche basate su potenziali interatomici offrono una risoluzione ad alta fedeltà delle instabilità reticolari, della nucleazione dei difetti e dell'evoluzione della microstruttura. Tuttavia, le dimensioni dei sistemi accessibili sono limitate alla scala nanometrica, rendendole inadatte all'analisi di componenti strutturali su scala micrometrica e delle loro condizioni al contorno. Sebbene i metodi multiscala concorrenti (ad esempio, il quasicontinuo) accoppino domini atomistici e continui nello spazio, spesso soffrono di sfide algoritmiche come forze spurie all'interfaccia ("forze fantasma") e difficoltà nel mantenere la coerenza energetica globale. I metodi gerarchici, che si basano sull'omogeneizzazione, evitano interfacce spaziali dirette ma tradizionalmente richiedono una legge costitutiva continua alla microscala. La sfida affrontata in questo lavoro è stabilire un accoppiamento gerarchico in cui la microscala è governata strettamente dalla fisica atomistica discreta (potenziali interatomici) piuttosto che da una legge continua fenomenologica, mantenendo al contempo coerenza energetica e robustezza per l'analisi strutturale non lineare.

Metodologia
Gli autori propongono un framework atomistico-continuo a due scale basato sull'omogeneizzazione computazionale e sul principio di equivalenza energetica di Hill–Mandel. L'approccio accoppia un problema ai valori al contorno non lineare agli elementi finiti (FE) alla microscala (livello strutturale) con problemi di cella di statica molecolare (MS) periodici nei punti di quadratura della mesh FE.

• Transizione di Scala: Il gradiente di deformazione continuo ($F$) guida il volume elementare rappresentativo (RVE) atomistico. Vengono applicate condizioni al contorno periodiche (PBC) alla cella atomistica. Il solver atomistico esegue una minimizzazione dell'energia a temperatura zero per trovare la configurazione atomica di equilibrio per la deformazione prescritta.
• Trasferimento di Sforzi e Tangente: A differenza dei metodi FE$^2$ standard che si basano su una funzione di energia immagazzinata, la cella atomistica restituisce gli sforzi di Cauchy mediati per volume derivati dal teorema del viriale. Questi vengono trasformati negli sforzi di Piola–Kirchhoff del primo ordine ($P$) e del secondo ordine ($S$) coniugati al lavoro, richiesti dalla formulazione FE lagrangiana.
• Tangente Algoritmica: Per garantire iterazioni di Newton robuste alla scala grossolana, viene calcolato un modulo tangente efficace. Ciò viene ottenuto mediante differenziazione numerica (differenza in avanti) della risposta agli sforzi rispetto al gradiente di deformazione. Ciò comporta la perturbazione del gradiente di deformazione, l'esecuzione di rilassamenti MS indipendenti per ogni perturbazione e il calcolo degli incrementi di sforzo risultanti.
• Implementazione: Il framework integra LAMMPS come solver MS all'interno di un workflow scalabile basato su MPI. Il parallelismo è sfruttato attraverso i punti di quadratura, dove ogni punto di Gauss attiva un insieme indipendente di simulazioni atomistiche (una per lo stato di base e sei per le perturbazioni tangenti).

Contributi Chiave

1. Framework Energeticamente Coerente: Gli autori stabiliscono una transizione dalla scala nanometrica a quella micrometrica basata strettamente sul postulato di Hill–Mandel, garantendo l'equivalenza energetica tra la potenza virtuale macroscopica e la potenza virtuale microscopica mediata per volume.
2. Interfaccia Discreto-Continuo: Viene definita un'interfaccia specifica in cui la deformazione continua guida un RVE atomistico periodico, e il solver atomistico restituisce sforzi mediati per volume in forma coniugata al lavoro insieme a moduli tangenti efficaci.
3. Robustezza in Regimi Non Lisci: Il lavoro dimostra che un solver Newton statico-implicito può rimanere robusto anche quando la risposta atomistica non è liscia a causa della nucleazione di difetti e dell'instabilità. Il metodo recupera con successo una convergenza quasi quadratica nelle iterazioni finali dei passi di carico critici.
4. Implementazione Pratica: Lo studio delinea una strategia per integrare blocchi di simulazione consolidati (LAMMPS) in un workflow scalabile, gestendo il costo computazionale di molteplici risoluzioni atomistiche concorrenti per iterazione macroscopica.

Risultati Numerici
Il framework è stato testato su rame monocristallino (potenziale EAM) in diversi scenari:

• Trazione e Compressione Uniaxiale: Il modello ha catturato una marcata asimmetria tra trazione e compressione. In trazione, l'instabilità ha portato alla formazione di difetti di impilamento e lamelle correlate ai geminati. In compressione, la microstruttura è evoluta in fasci ampi e stratificati di difetti di impilamento e regioni a impacchettamento esagonale compatto (hcp). I risultati hanno mostrato che l'insorgenza dell'instabilità dipende dalle dimensioni della scatola di simulazione, mentre la risposta elastica pre-instabilità è indipendente dalle dimensioni.
• Caricamento Ciclico: Sotto caricamento ciclico inverso, il sistema ha mostrato una transizione da un regime elastico iniziale ad alta forza a uno stato governato da un'architettura di difetti metastabili "addestrata". Ciò ha portato a cicli di isteresi ampi e a una rapida stabilizzazione ciclica dopo il primo evento di instabilità.
• Flessione di una Mensola: È stata simulata una trave su scala micrometrica (lunghezza 10 µm). Il framework ha risolto una transizione localizzata elastico-plastica in cui la nucleazione di difetti è avvenuta bruscamente in nanodomini vicino all'estremità incassata (fibre inferiori), mentre il resto della trave è rimasto elastico. Il solver Newton ha gestito efficacemente i bruschi cambiamenti di rigidità associati a queste instabilità localizzate.

Significato e Affermazioni
Il documento afferma di fornire il primo upscaling da atomistico a continuo che sia simultaneamente: (i) basato sul principio energetico rigoroso di Hill–Mandel, (ii) compatibile con una formulazione agli elementi finiti statica-implicita, e (iii) dotato di una tangente efficace adatta ai solver Newton in regimi fortemente non lineari.

Gli autori sottolineano che il loro approccio estende la modellazione atomistica basata su potenziali a scale di lunghezza strutturali (micrometri) inaccessibili alla simulazione atomistica diretta in un contesto quasi-statico e atermico. Una scoperta critica è che il solver Newton alla scala grossolana rimane robusto e recupera una convergenza quasi quadratica anche quando la risposta alla scala fine subisce bruschi cambiamenti di rigidità a causa della nucleazione di difetti. Ciò suggerisce che le informazioni sugli sforzi e sulla tangente coerenti dal punto di vista energetico provenienti dalle celle atomistiche rimangono strumenti utili per l'analisi strutturale anche in scenari altamente non lineari e instabili.

Il lavoro è esplicitamente limitato alle risposte quasi-statiche e atermiche a temperatura zero. Gli autori notano che future estensioni potrebbero coinvolgere modelli a temperatura finita, condizioni al contorno alternative e l'integrazione di potenziali interatomici basati sull'apprendimento automatico.