The industrial application motivating this work is the fatigue computation of aircraft engines' high-pressure turbine blades. The material model involves nonlinear elastoviscoplastic behavior laws, for which the parameters depend on the temperature. For this application, the temperature loading is not accurately known and can reach values relatively close to the creep temperature: important nonlinear effects occur and the solution strongly depends on the used thermal loading. We consider a nonlinear reduced order model able to compute, in the exploitation phase, the behavior of the blade for a new temperature field loading. The sensitivity of the solution to the temperature makes the classical unenriched proper orthogonal decomposition method fail. In this work, we propose a new error indicator, quantifying the error made by the reduced order model in computational complexity independent of the size of the high-fidelity reference model. In our framework, when the error indicator becomes larger than a given tolerance, the reduced order model is updated using one time step solution of the high-fidelity reference model. The approach is illustrated on a series of academic test cases and applied on a setting of industrial complexity involving 5 million degrees of freedom, where the whole procedure is computed in parallel with distributed memory.Computing lifetime predictions for high-pressure turbine blades is a challenging task: meshes involve large numbers of degrees of freedom to account for local structures such as the internal cooling channels, the behavior laws are strongly nonlinear with many internal variables, and a large number of cycles has to be computed. Besides, the temperature loading is poorly known in the outlet section of the combustion chamber. Our team has proposed in [12] a nonintrusive reduced order model (ROM) strategy in parallel computation with distributed memory to mitigate the runtime issues: a domain decomposition method is used to compute the first cycle, and the reduced order model is used to speed up the computation of the following cycles, which can be considered as a reduced order model-based temporal extrapolation. As pointed out in [12], errors are accumulated during this temporal extrapolation. Moreover, quantifying the uncertainty on the lifetime with respect to some statistical description of the temperature loading using an already constructed reduced order model would introduce additional errors. In this context, error indicator-based enrichment of reduced order models is the topic of the present work.Error estimation for reduced model predictions is a topic that receives interest in the scientific literature. The reduced basis method [32,28] for parametrized problems is a reduced order modeling method that intrinsically relies on efficient a posteriori error bounds of the error between the reduced prediction and the reference high-fidelity (HF) solution. This method consists in a greedy enrichment of a current reduced order basis by the high-fidelity solution at the parametric...