In modeling and numerically implementing a follower pressure in a geometrically nonlinear setting, one needs to compute the volume enclosed by a surface and its variation. For closed surfaces, the volume can be expressed as a surface integral invoking the divergence theorem. For periodic systems, widely used in computational physics and materials science, the enclosed volume calculation and its variation is more delicate and has not been examined before. Here, we develop simple expressions involving integrals on the surface, on its boundary lines, and point contributions. We consider two specific situations, a periodic tubular surface and a doubly periodic surface enclosing a volume with a nearby planar substrate, which are useful to model systems such as pressurized carbon nanotubes, supported lipid bilayers or graphene. We provide a set of numerical examples, which show that the familiar surface integral term alone leads to an incorrect volume evaluation and spurious forces at the periodic boundaries.
CitationRahimi, M.; Zhang, K.; Arroyo, M. Computing the volume enclosed by a periodic surface and its variation to model a follower pressure. "Computational Mechanics", 01 Març 2015, vol. 55, núm. 3, p. 519-525.
All rights reserved. This work is protected by the corresponding intellectual and industrial property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public communication or transformation of this work are prohibited without permission of the copyright holder. If you wish to make any use of the work not provided for in the law, please contact: email@example.com