Figure (E) presents the estimates of the linearization error for 𝜁, as a function of 𝑘 9CD , for case study a. The results were obtained according to the methodology in Section 6.5.2, considering 𝑃 = 10 s . For 𝑘 9CD > 10 M the linearization error decreases as 𝑘 9CD increases, as expected, with a linear behavior in a log-log plot. For 𝑘 9CD < 10 M , however, the linearization error behavior is erratic, as observed in Fig. (B)-(D), i.e., the physics of the problem is not effectively captured in this discretization level. Only after a certain onset, in this case identified as 𝑘 9CD > 10 M , the actual profiles for 𝜁 are adequately captured, and successive refinements are obtained as 𝑘 9CD increases, allowing for proper interpretation of the solution performance -in this case in terms of linearization error.