Finite Element & Finite Volume Methods

an isoparametric element

Real domains have curved boundaries — circular pipes, turbine blades, the surface of the Earth — but straight-sided triangles tile a curved boundary only crudely, leaving little chordal gaps. The isoparametric idea bends the elements themselves to follow the curve. The trick: define everything on one simple 'reference' element (a perfect unit triangle or square), then MAP it onto the actual, possibly curved, physical element — and use the SAME shape functions both to interpolate the solution AND to describe the geometry of the map. 'Iso-parametric' means 'same parameters' for both jobs.

Concretely, you work on a master element with coordinates (xi, eta) where the shape functions N_i(xi, eta) are simple polynomials. The physical coordinates of a point come from x = sum N_i(xi, eta) x_i, where x_i are the actual node positions — the same N_i that interpolate the solution u = sum N_i u_i. If you add mid-edge nodes and use quadratic N_i, the edges of the mapped element become parabolas, hugging a curved boundary far better than straight lines. All the element integrals (stiffness, mass, load) are then computed on the tidy reference element, pulling factors through the map: an integral over the physical element becomes an integral over the master element weighted by the determinant of the Jacobian of the map, evaluated by Gaussian quadrature at a few points. This is why numerical integration is woven into finite element codes — the mapped integrands are rarely polynomials you can integrate by hand.

Isoparametric elements are how production finite element software handles real geometry and high accuracy together, but two honesties matter. First, the mapping's Jacobian must stay positive everywhere inside the element — if an element is too distorted the map folds over, the Jacobian determinant goes negative, and the element becomes invalid (garbage results, not just inaccuracy). Second, the quadrature must be accurate enough: too few Gauss points UNDER-integrates and can introduce spurious zero-energy 'hourglass' modes, while too many wastes time; choosing the right rule for each element type is part of the craft.

A quadratic triangle has nodes at the 3 corners plus 3 edge midpoints. Placing the midpoint nodes ON a circular boundary and mapping with quadratic shape functions bends the element's edges into parabolas that follow the circle closely — capturing curvature that straight-sided linear triangles could only approximate by piling on tiny flat elements.

Same shape functions map geometry and interpolate solution; element integrals are done on a reference element by quadrature.

The mapping's Jacobian determinant must stay strictly positive throughout the element — an over-distorted element folds the map, turns the Jacobian negative, and produces nonsense, not merely inaccuracy. Too-low quadrature can also seed spurious hourglass modes.

Also called
isoparametric mappingcurved elementreference element mapping等參映射參考元素