One-shot Parareal Approach for Topology Optimisation of Transient Heat Flow

TL;DR

This paper introduces a one-shot Parareal-based method for topology optimization of transient heat flow, achieving significant speedup and comparable design quality, with modifications to handle adjoint analysis and iterative refinement.

Contribution

The paper proposes a novel one-shot Parareal approach for topology optimization that reduces iteration counts and improves computational efficiency in transient heat problems.

Findings

Achieved up to 5x speedup with 16 threads.

Produced designs similar to reference methods in quality.

Identified coarse propagator computation as a bottleneck.

Abstract

This paper presents a method of performing topology optimisation of transient heat conduction problems using the parallel-in-time method Parareal. To accommodate the adjoint analysis, the Parareal method was modified to store intermediate time steps. Preliminary tests revealed that Parareal requires many iterations to achieve accurate results and, thus, achieves no appreciable speedup. To mitigate this, a one-shot approach was used, where the time history is iteratively refined over the optimisation process. The method estimates objectives and sensitivities by introducing cumulative objectives and sensitivities and solving for these using a single iteration of Parareal, after which it updates the design using the Method of Moving Asymptotes. The resulting method was applied to a test problem where a power mean of the temperature was minimised. It achieved a peak speedup relative to a sequential reference method of $5\times$ using 16 threads. The resulting designs were similar to the one found by the reference method, both in terms of objective values and qualitative appearance. The one-shot Parareal method was compared to the Parallel Local-in-Time method of topology optimisation. This revealed that the Parallel Local-in-Time method was unstable for the considered test problem, but it achieved a peak speedup of $12\times$ using 32 threads. It was determined that the dominant bottleneck in the one-shot Parareal method was the time spent on computing coarse propagators.