Monte-Carlo study of scaling exponents of rough surfaces and correlated percolation

We calculate the scaling exponents of the two-dimensional correlated percolation cluster's hull and unscreened perimeter. Correlations are introduced through an underlying correlated random potential, which is used to define the state of bonds of a two-dimensional bond percolation model. Monte-Carlo simulations are run and the values of the scaling exponents are determined as functions of the Hurst exponent H in the range -0.75 <= H <= 1. The results confirm the conjectures of earlier studies

mdFoam+: Advanced molecular dynamics in OpenFOAM

This paper introduces mdFoam+, which is an MPI parallelised molecular dynamics (MD) solver implemented entirely within the OpenFOAM software framework. It is open-source and released under the same GNU General Public License (GPL) as OpenFOAM. The source code is released as a publicly open software repository that includes detailed documentation and tutorial cases. Since mdFoam+ is designed entirely within the OpenFOAM C++ object-oriented framework, it inherits a number of key features. The code is designed for extensibility and flexibility, so it is aimed first and foremost as an MD research tool, in which new models and test cases can be developed and tested rapidly. Implementing mdFoam+ in OpenFOAM also enables easier development of hybrid methods that couple MD with continuum-based solvers. Setting up MD cases follows the standard OpenFOAM format, as mdFoam+ also relies upon the OpenFOAM dictionary-based directory structure. This ensures that useful pre- and post-processing capabilities provided by OpenFOAM remain available even though the fully Lagrangian nature of an MD simulation is not typical of most OpenFOAM applications. Results show that mdFoam+ compares well to another well-known MD code (e.g. LAMMPS) in terms of benchmark problems, although it also has additional functionality that does not exist in other open-source MD codes