Skip to content

g7fernandes/molecular_dynamics

Folders and files

NameName
Last commit message
Last commit date

Latest commit

Feb 13, 2020
8b3e1de · Feb 13, 2020

History

77 Commits
Aug 17, 2019
Aug 8, 2019
Feb 13, 2020
Oct 3, 2019
Sep 2, 2019
Dec 10, 2019
Dec 10, 2019
Aug 31, 2019
Sep 4, 2019
Jul 16, 2019
Sep 11, 2019
Sep 13, 2019
Dec 12, 2019
Aug 29, 2019
Jan 14, 2020
May 28, 2019
May 28, 2019
Aug 7, 2019
May 28, 2019
Sep 4, 2019
May 28, 2019
May 28, 2019
May 28, 2019
Nov 5, 2019
Nov 5, 2019
Sep 11, 2019
Sep 11, 2019
Sep 3, 2019
Jan 14, 2020
Sep 8, 2019
Dec 10, 2019
Aug 20, 2019
Sep 24, 2019
Sep 20, 2019
Sep 15, 2019
May 28, 2019
Jun 17, 2019
Sep 24, 2019
Aug 13, 2019
Jul 16, 2019
Jul 19, 2019
Aug 13, 2019
Nov 5, 2019

Repository files navigation

Molecular Dynamics

This program simulates molecular interactions of monoatomic molecules and particles with many molecules in a 2D region of calculus. The development of this code was part of a graduante reseach project of the author's Master Degree at the Faculty of Mechanical Engineering at University of Campinas.

Compilation

The Fortran code can be compiled using mpiifort (Intel Fortran Compiler) or mpifort (GCC Fortran). You may change the make file to choose the compiler that you wish. The main file is "lennard.f90"

Requirements

GNU Fortran

sudo apt-get install gfortran # Ubuntu Based distros
sudo zypper in gcc-fortran # OpenSUSE

Or Intel Fortran Compiler

If using GNU Fortran, install also the OpenMPI package

sudo apt-get install openmpi-bin libopenmpi-dev # Ubuntu 
sudo zypper in openmpi openmpi-devel # OpenSUSE

For the Python code it is necessary the pyevtk package to convert the csv files to vtk. The package progressbar2 is optional.

For the visualization, you can use Paraview or some program that reads VTK files.

Input and Output

The input needs a configuration file named settings.ini , which contains the physical properties and numerial parameters for the simulation, and files containing the positions of the molecules. Optionally a file containing the initial velocity of the particles can also be provided.

The output contains the position and velocities of the molecules. Optionally the angle of rotation of particles, and forces related to intermolecular interactions can also be extracted.

Postprocessing

The program automatically converts the .csv results to VTK, so one may use this result to postprocess on paraview. Some python scripts are provided to calculate concentration, distribution entropy, viscosity and other things. The theory behind these program is summarized at the Master Thesis that will be published soon.

References

Some of the code on this repository was based on code published by Jason R. Blevins ACM Fortran Forum 28(3), 2–7, 2009. https://github.com/jannisteunissen/config_fortran Numerical Simulation in Molecular Dynamics: Numerics, Algorithms, Parallelization, Applications, by Michael Griebel (Author), Stephan Knapek (Contributor), Gerhard Zumbusch (Contributor)

Acknowledgments

The author to acknowledge the support of this work by FAPESP-Brazil under the grant number 2017/05643-8 and FAEPEX Unicamp.

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages