“…In our discretization of equation [ Chobanyan et al ., ], the computational domain is first geometrically tessellated using Lagrange‐type generalized curved parametric hexahedra of arbitrary geometrical orders K u , K v , and K w ( K u , K v , K w ≥ 1), shown in Figure and analytically described as [ Ilić and Notaroš , ] where r ijk = r ( u i , v j , w k ) are the position vectors of interpolation nodes and represent Lagrange interpolation polynomials in the u coordinate, with u l being the uniformly spaced interpolating nodes in u defined as u l = (2 l − K u )/ K u , l = 0, 1, …, K u , and similarly for and . Curvilinear hexahedral modeling facilitates volumetric meshes that can employ very large elements, which is consistent with the double‐higher‐order large‐domain (entire‐domain) VIE paradigm.…”