Automated Parallel and Body-Fitted Mesh Generation in Finite Element Simulation of Macromolecular Systems

Automated Parallel and Body-Fitted Mesh Generation in Finite Element Simulation of Macromolecular Systems

Year:    2016

Communications in Computational Physics, Vol. 19 (2016), Iss. 3 : pp. 582–602

Abstract

Mesh generation is a bottleneck for finite element simulations of biomolecules. A robust and efficient approach, based on the immersed boundary method proposed in [8], has been developed and implemented to generate large-scale mesh body-fitted to molecular shape for general parallel finite element simulations. The molecular Gaussian surface is adopted to represent the molecular surface, and is finally approximated by piecewise planes via the tool phgSurfaceCut in PHG [43], which is improved and can reliably handle complicated molecular surfaces, through mesh refinement steps. A coarse background mesh is imported first and then is distributed into each process using a mesh partitioning algorithm such as space filling curve [5] or METIS [22]. A bisection method is used for the mesh refinements according to the molecular PDB or PQR file which describes the biomolecular region. After mesh refinements, the mesh is optionally repartitioned and redistributed for load balancing. For finite element simulations, the modification of region mark and boundary types is done in parallel. Our parallel mesh generation method has been successfully applied to a sphere cavity model, a DNA fragment, a gramicidin A channel and a huge Dengue virus system. The results of numerical experiments show good parallel efficiency. Computations of electrostatic potential and solvation energy also validate the method. Moreover, the meshing process and adaptive finite element computation can be integrated as one PHG project to avoid the mesh importing and exporting costs, and improve the convenience of application as well.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.161114.021015a

Communications in Computational Physics, Vol. 19 (2016), Iss. 3 : pp. 582–602

Published online:    2016-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    21

Keywords:   

  1. High-order Numerical Quadratures in a Tetrahedron with an Implicitly Defined Curved Interface

    Cui, Tao | Leng, Wei | Liu, Huaqing | Zhang, Linbo | Zheng, Weiying

    ACM Transactions on Mathematical Software, Vol. 46 (2020), Iss. 1 P.1

    https://doi.org/10.1145/3372144 [Citations: 13]
  2. Frontiers in biomolecular mesh generation and molecular visualization systems

    Gui, Sheng | Khan, Dawar | Wang, Qin | Yan, Dong-Ming | Lu, Ben-Zhuo

    Visual Computing for Industry, Biomedicine, and Art, Vol. 1 (2018), Iss. 1

    https://doi.org/10.1186/s42492-018-0007-0 [Citations: 6]
  3. A flux-based moving mesh method applied to solving the Poisson–Nernst–Planck equations

    Lv, Minrui | Lu, Benzhuo

    Journal of Computational Physics, Vol. 513 (2024), Iss. P.113169

    https://doi.org/10.1016/j.jcp.2024.113169 [Citations: 1]