Vai al contenuto principale
it

Gradienti e limitatori su una griglia scadente

In un solutore a celle centrate ogni flusso poggia su una stima di gradiente. È qui che la qualità della griglia entra per prima nella soluzione.

Nel metodo ai volumi finiti le incognite sono immagazzinate nei centri cella, ma i flussi servono sulle facce delle celle. Il ponte fra i due è la stima del gradiente, ed è lì che la qualità della griglia entra per prima nella soluzione — su celle oblique un errore di gradiente diventa direttamente un errore di flusso.

Due metodi comuni

Minimi quadrati. Si adatta un piano ai centri delle celle vicine. Più robusto su griglie irregolari; si indebolisce quando i vicini sono addensati in una sola direzione.

Green-Gauss. Si integra sui valori di faccia. Economico e accurato su griglia regolare; su celle oblique l'interpolazione del valore di faccia degrada e l'errore cresce.

La scelta dipende dalla qualità della griglia ed è una delle voci da controllare su una soluzione che non converge. Si veda la nota sulle impostazioni del solutore.

Che cosa proteggono i limitatori, e quanto costano

La ricostruzione del secondo ordine produce sovraelongazioni non fisiche ai gradienti ripidi: energia cinetica turbolenta negativa, frazioni volumetriche oltre il 100 %, temperature sotto lo zero assoluto. Un limitatore le taglia.

Il costo è questo: ovunque il limitatore intervenga, la soluzione scende localmente al primo ordine e si genera diffusione numerica. Più il limitatore è aggressivo, più la soluzione è stabile e sfocata.

Il punto pratico

Guardate dove interviene il limitatore. Se interviene a una transizione davvero ripida — un'onda d'urto o una superficie libera — è normale. Se interviene in mezzo a un flusso regolare, è un sintomo di qualità della griglia: in quella zona ci sono celle oblique o salti bruschi di volume.

La maggior parte dei casi in cui la qualità della griglia sembra buona in media mentre la soluzione degrada viene proprio da qui — ciò che va esaminato sono le celle peggiori. Si veda la nota sul tipo di cella.

Riferimenti

  1. Ferziger, J. H., Perić, M. & Street, R. L. (2020). Computational Methods for Fluid Dynamics, 4th ed. Springer. DOI: 10.1007/978-3-319-99693-6
  2. Patankar, S. V. (1980). Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing. ISBN 978-0-89116-522-4.

Parliamo del tuo progetto

Definiamo insieme il sistema da analizzare e la domanda a cui rispondere.

Parliamo del tuo progetto