Zanurzam się w fascynujący świat analizy elementów skończonych i chciałbym rozwiązać duży problem termomechaniczny (tylko mechaniczny efekt termiczny , brak sprzężenia zwrotnego).
W przypadku problemu mechanicznego zrozumiałem już odpowiedź Geoffa , że muszę użyć iteracyjnego solvera ze względu na rozmiar mojej siatki. W odpowiedzi Matta przeczytałem dalej , że wybór właściwego algorytmu iteracyjnego jest trudnym zadaniem.
Pytam tutaj, czy istnieją jakieś doświadczenia dotyczące dużych trójwymiarowych problemów liniowo-elastycznych, które pomogłyby mi zawęzić moje poszukiwania najlepszej wydajności? W moim przypadku jest to struktura z cienkimi, wzorzystymi filmami i nieregularnie rozmieszczonymi materiałami (zarówno o wysokim CTE, jak i niskim CTE). W tej analizie termomechanicznej nie ma dużych odkształceń. Mogę korzystać z HPC mojego uniwersytetu [1.314 węzłów, z 2 procesorami AMD Opteron (każdy 2,2 GHz / 8 rdzeni)].
Myślę, że PETScmoże zawierać coś interesującego, zwłaszcza algorytmy, które dokonują pewnego rodzaju dekompozycji domen (FETI, multigrid), ale opcje mnie trochę przytłaczają i nie mam doświadczenia. Podoba mi się również wyrażenie „geometrycznie poinformowane warunki wstępne” , ale nie jestem pewien, czy to mi pomaga. Nie znalazłem jeszcze czegoś, co skupiałoby się na mechanice ciągłości liniowej.
Silne skalowanie (Amdahl) jest bardzo ważne w mojej aplikacji, ponieważ moi partnerzy przemysłowi nie mogą długo czekać na wyniki symulacji. Zdecydowanie doceniam nie tylko odpowiedzi, ale także rekomendacje do dalszego czytania w komentarzach.