  An interior penalty discontinuous Galerkin method is devised to approximate
minimizers of a linear folding model by discontinuous isoparametric finite
element functions that account for an approximation of a folding arc. The
numerical analysis of the discrete model includes an a priori error estimate in
case of an accurate representation of the folding curve by the isoparametric
mesh. Additional estimates show that geometric consistency errors may be
controlled separately if the folding arc is approximated by piecewise
polynomial curves. Various numerical experiments are carried out to validate
the a priori error estimate for the folding model.
