An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials
The purpose of the present letter is to propose an efficient, accurate and robust FFT-based solver for transient diffusion in heterogeneous materials with “realistic” BC, taking advantage of two recent advances in terms of boundary conditions and finite difference schemes to overcome their actual li...
Saved in:
| Main Author: | |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Académie des sciences
2025-01-01
|
| Series: | Comptes Rendus. Mécanique |
| Subjects: | |
| Online Access: | https://comptes-rendus.academie-sciences.fr/mecanique/articles/10.5802/crmeca.281/ |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1825205960856043520 |
|---|---|
| author | Gélébart, Lionel |
| author_facet | Gélébart, Lionel |
| author_sort | Gélébart, Lionel |
| collection | DOAJ |
| description | The purpose of the present letter is to propose an efficient, accurate and robust FFT-based solver for transient diffusion in heterogeneous materials with “realistic” BC, taking advantage of two recent advances in terms of boundary conditions and finite difference schemes to overcome their actual limitations (periodic BC and spurious oscillations). It is an essential step towards couplings between mechanics and other physics (such as the diffusion of species) through FFT-based solvers. Discrete Trigonometric Transform are used to implement non-periodic boundary conditions, and a finite difference (FD) scheme recently proposed by Finel is advantageously compared to the common hexahedral FD scheme. “Accurate” refers to two properties: accurate in term of locality with a small size of Finite Difference pencil to capture fluctuations around material heterogeneities, and accurate in term of precision with the absence of spurious spatial oscillations (at least in the reported cases with well-separated inclusions). The “robustness” is here associated to the stability of the solver, especially associated to the implicit time integration method. The description of the method focuses on thermal diffusion but applies to any other similar diffusion process (with the same type of parabolic equation). As a by-product, the FD scheme proposed by Finel is introduced in a more general (for mixing finite different schemes) and simple way (no introduction of FCC subgrids), extending its usage to any type of grid parity (not only even grids). |
| format | Article |
| id | doaj-art-d1d27daded754d0084d168d94eff413e |
| institution | Kabale University |
| issn | 1873-7234 |
| language | English |
| publishDate | 2025-01-01 |
| publisher | Académie des sciences |
| record_format | Article |
| series | Comptes Rendus. Mécanique |
| spelling | doaj-art-d1d27daded754d0084d168d94eff413e2025-02-07T13:49:01ZengAcadémie des sciencesComptes Rendus. Mécanique1873-72342025-01-01353G111312510.5802/crmeca.28110.5802/crmeca.281An accurate and robust FFT-based solver for transient diffusion in heterogeneous materialsGélébart, Lionel0https://orcid.org/0000-0002-1387-1978Université Paris Saclay, CEA, SRMA, 91191, Gif/Yvette, FranceThe purpose of the present letter is to propose an efficient, accurate and robust FFT-based solver for transient diffusion in heterogeneous materials with “realistic” BC, taking advantage of two recent advances in terms of boundary conditions and finite difference schemes to overcome their actual limitations (periodic BC and spurious oscillations). It is an essential step towards couplings between mechanics and other physics (such as the diffusion of species) through FFT-based solvers. Discrete Trigonometric Transform are used to implement non-periodic boundary conditions, and a finite difference (FD) scheme recently proposed by Finel is advantageously compared to the common hexahedral FD scheme. “Accurate” refers to two properties: accurate in term of locality with a small size of Finite Difference pencil to capture fluctuations around material heterogeneities, and accurate in term of precision with the absence of spurious spatial oscillations (at least in the reported cases with well-separated inclusions). The “robustness” is here associated to the stability of the solver, especially associated to the implicit time integration method. The description of the method focuses on thermal diffusion but applies to any other similar diffusion process (with the same type of parabolic equation). As a by-product, the FD scheme proposed by Finel is introduced in a more general (for mixing finite different schemes) and simple way (no introduction of FCC subgrids), extending its usage to any type of grid parity (not only even grids).https://comptes-rendus.academie-sciences.fr/mecanique/articles/10.5802/crmeca.281/FFTHeterogeneous materialsTransient diffusionFinite differenceTetrahedra |
| spellingShingle | Gélébart, Lionel An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials Comptes Rendus. Mécanique FFT Heterogeneous materials Transient diffusion Finite difference Tetrahedra |
| title | An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials |
| title_full | An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials |
| title_fullStr | An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials |
| title_full_unstemmed | An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials |
| title_short | An accurate and robust FFT-based solver for transient diffusion in heterogeneous materials |
| title_sort | accurate and robust fft based solver for transient diffusion in heterogeneous materials |
| topic | FFT Heterogeneous materials Transient diffusion Finite difference Tetrahedra |
| url | https://comptes-rendus.academie-sciences.fr/mecanique/articles/10.5802/crmeca.281/ |
| work_keys_str_mv | AT gelebartlionel anaccurateandrobustfftbasedsolverfortransientdiffusioninheterogeneousmaterials AT gelebartlionel accurateandrobustfftbasedsolverfortransientdiffusioninheterogeneousmaterials |