Imported meshes

A mesh read from a file must behave exactly like a generated one, including unstructured topology, where neighbouring elements use different local axes. (i) Every Gmsh fixture of the test suite is read; the files are genuine Gmsh output (tools/make_mesh_fixtures.py), and on the straight-sided ones every high-order node must be the multilinear image of the element's vertices, which tests the node-ordering tables for quadrilaterals of degree 1 to 5 and hexahedra of degree 1 to 4. (ii) On curved unstructured meshes from Gmsh in 2-D and 3-D, and on a 3-D mesh whose elements were randomly re-oriented and passed through a file, a uniform flow must remain steady, the residuals of a perturbed state must sum to zero, and the inviscid residual of a continuous state must not depend on the Riemann solver away from the boundary (the jump across every face is zero only if its nodes are matched with the right ones of the neighbour; the viscous terms cancel in this comparison, so it is made for the Euler equations only). The coupling of the viscous terms across faces (gradient traces of the neighbour, lifted jumps) is tested separately, for both sets of equations: a rough state with varying velocity and temperature and a jump on every face must give the same residual, node by node, when every element of the mesh is stored in other, randomly chosen local axes. (iii) The coarse annulus of the cylinder case, written to a Gmsh file after randomly re-orienting its elements and read back through a case specification, must give the force coefficients of the generated mesh to solver tolerance. (iv) Laminar flow past a cylinder in a channel (Re = 40, M = 0.2, slip walls) is solved on the unstructured example mesh at degree 2, and the adjoint derivative of the drag coefficient with respect to the Mach number is compared with a central finite difference of re-converged flows.

• usable Gmsh fixtures read15 of 15expected allpass
• unusable Gmsh fixtures (triangles, binary, untagged) refused with a message3 of 3expected allpass
• straight-sided fixtures cover quadrilaterals of degree 1–5 and hexahedra of degree 1–49expected = 9 (dimension, degree) pairspass
• node ordering: largest distance of a node of a straight-sided mesh from the multilinear map of its vertices4.44e-15expected < 1e-13 (round-off; coordinates of order one)pass
• straight-sided meshes: relative error of the area / volume3.06e-16expected < 1e-12 (the quadrature is exact for multilinear elements)pass
• every straight-sided fixture is unstructured: fewest face pairs with a non-identity orientation44expected > 0pass
• curved quadrilateral meshes (circular hole): area error at geometry degree pg + 1 relative to pg, largest0.1473expected < 1 (the boundary interpolant converges with the degree)pass
• MSH 2.2 and MSH 4.1 files of one periodic mesh: largest node difference0expected < 1e-12pass
• periodic node pairs of the file became topology (slave nodes)13expected > 0, and no boundary group on the periodic sidespass
• hole_quad_p3.msh (Gmsh, curved, 2-D), inviscid: perturbed state is not trivial (max |Rw|)0.4158expected > 1e-3pass
• hole_quad_p3.msh (Gmsh, curved, 2-D), viscous: perturbed state is not trivial (max |Rw|)0.4359expected > 1e-3pass
• hole_quad_p3.msh (Gmsh, curved, 2-D): face pairs with a non-identity orientation48expected > 0pass
• hole_hex_p2.msh (Gmsh, curved, 3-D), inviscid: perturbed state is not trivial (max |Rw|)0.6203expected > 1e-3pass
• hole_hex_p2.msh (Gmsh, curved, 3-D), viscous: perturbed state is not trivial (max |Rw|)0.3062expected > 1e-3pass
• hole_hex_p2.msh (Gmsh, curved, 3-D): face pairs with a non-identity orientation72expected > 0pass
• periodic warped 3x3x3 box, re-oriented, written and read back, inviscid: perturbed state is not trivial (max |Rw|)0.1151expected > 1e-3pass
• periodic warped 3x3x3 box, re-oriented, written and read back, viscous: perturbed state is not trivial (max |Rw|)0.1204expected > 1e-3pass
• periodic warped 3x3x3 box, re-oriented, written and read back: face pairs with a non-identity orientation74expected > 0pass
• re-oriented 3-D box read from a file: relative face orientations present8expected = 8 (all)pass
• free stream on imported meshes: largest max |Rw|4.44e-15expected < 1e-12pass
• conservation on imported meshes: largest |sum of Rw| of a perturbed state1.38e-14expected < 1e-12 (residuals of order 0.1 to 1)pass
• neighbour matching on imported meshes (inviscid): largest |Rw(Roe) − Rw(LLF)| off the boundary for a continuous state9.99e-16expected < 1e-12 (zero jump across every correctly matched face; residuals of order 0.1 to 1)pass
• neighbour coupling on imported meshes, Euler and Navier–Stokes: largest difference of the residual of a rough state between two descriptions of the mesh (elements re-oriented)2.48e-14expected < 1e-12 (round-off; residuals of order 0.1 to 1)pass
• the rough state exercises the viscous terms: smallest max |Rw(Navier–Stokes) − Rw(Euler)|0.02271expected > 1e-3 (velocity and temperature vary; far above the bound of the coupling check)pass
• generated and imported annulus: both solves convergedyesexpected truepass
• imported annulus: face pairs with a non-identity orientation42expected > 0 (elements were re-oriented before writing)pass
• force and moment coefficients, |imported − generated| (largest)2.13e-14expected < 1e-10 (both solves reduce the residual by 1e-12)pass
• cylinder in a channel (examples/meshes/cylinder_channel_p3.msh), Re = 40, p = 2: convergedyesexpected truepass
• cylinder in a channel: dcd/dmach, adjoint vs central finite difference (rel.)1.18e-07expected < 1e-6pass

Run notes Measured in this run: the largest free-stream residual on the imported meshes was 4.4e-15, the largest conservation defect 1.4e-14 and the largest dependence of the inviscid residual of a continuous state on the Riemann solver 1.0e-15. The residual of a rough state changed by at most 2.5e-14 when the elements were re-oriented, with and without viscous terms; the viscous terms contributed between 2.3e-02 and 1.4e-01 to that residual. Conservation alone would not notice a wrong index map between two straight faces (their normals are constant); the Riemann-solver comparison (trace matching, an inviscid test: the viscous terms cancel in it), the comparison of two descriptions (which includes the gradient traces and lifted jumps of the viscous discretization) and the solves of parts (iii) and (iv) do. The annulus written to a file with re-oriented elements (42 of 84 face pairs with a non-identity orientation) gave force coefficients within 2.1e-14 of those on the generated mesh. On the example mesh the adjoint derivative of the drag coefficient with respect to the Mach number differs from the central difference by a relative 1.2e-07. The Gmsh fixtures contain only two of the eight relative face orientations of hexahedra (see the table); all eight occur in the re-oriented box, which is written by Gyre's own writer, not by Gmsh.

gyre 0.4.0 · 39979fcf90 · 2026-11-12 11:46 · cpu · quick suite · record (JSON)Same code as this build