Skip to content

Latest commit

 

History

History
100 lines (79 loc) · 3.87 KB

README.md

File metadata and controls

100 lines (79 loc) · 3.87 KB

DynEarthSol3D_energy : A numerical tool for tectonic modeling with thermomechanics

Basic Information:

DynEarthSol3D is a finite element code that solves the momentum balance and the heat transfer in Lagrangian form using unstructured meshes. It can be used to study the long-term deformation of Earth's lithosphere and problems alike.

The purpose of this fork of DES3D (https://bitbucket.org/tan2/dynearthsol3d) is to show that libadaptivity from Fluidity (https://github.com/FluidityProject/fluidity), a self-contained library for anisotropic adaptive mesh refinement, works well as a mesh optimizer for DES3D's remeshing.

Thermomechanics module:

In this code i have added the full energy balance equation with mass conservation. For anything about these implementations please feel free to contact :

Sabber Ahamed [email protected] Senior Data Scientist Bridgestone America Nashville, TN, USA

Installation

  • You will need a recent C++ compiler that supports C++11 standard. (GNU g++ 4.4 or newer version will suffice.)
  • You will need a recent version of Boost::Program_options library (1.42 or newer version). Instructions for building the library: -- Download the source code from www.boost.org -- In the untarred source directory, run "./bootstrap.sh" -- In the same directory, run "./b2 --with-program_options -q" to build the library.
  • You will need Python 2.6+ or 3.2+ and the Numpy package.

    For libadaptiviy, you further need

  • VTK (v.5.10 tested) built from source or development packages.
  • MPI (openmpi-1.6.1 tested).

    Build procedure:

    • libadaptivity
  • Run "LDFLAGS=-L${VTK_LIBDIR} ./configure" in libadaptivity
  • Run "make".
  • Run "make" in tests/. Depending on the flavor of MPI, one might need to add '-lmpi_f77' to LIBS in Makefile.
    • DES3D
  • Edit 'Makefile', 1) modify BOOST_ROOT_DIR if you manually built or installed boost library. If you followed the instructions above to build Boost::Program_options library, set BOOST_ROOT_DIR to the untarred boost directory. 2) turn on 'useadapt' to use libadaptivity. 3) Set 'ndims = 3'. libadaptivity works only for 3D. 4) Check if VTK_INC path is correct. 5) Copy LIBS in libadaptivity/tests/Makefile to LIBADAPTIVITY_LIBS. 6) Make sure that the path to libadaptivity.a is correctly set in the beginning of LIBADAPTIVITY_LIBS. $(LIBADAPTIVITY_LIB)/libadaptivity.a will do.
  • Run "make opt=0" to build a debugging executable.
  • Run "make openmp=0" to build the executable without OpenMP. This is necessary to debug the code under valgrind.

    How to run

  • Execute "dynearthsol2d inputfile".
  • Several example input files are provided under 'examples/' directory. The format of the input file is described in 'examples/defaults.cfg'.
  • Benchmark cases with analytical solution can be found under 'benchmarks/' directory.
  • Execute the executable with '-h' flag to see the available input parameters and their descriptions.

    How to plot

  • Run "2vtk.py modelname" to convert the binary output to VTK files.
  • Execute 2vtk.py with '-h' flag to see more usage information.
  • Some of the simulation outputs might be disabled. Edit 2vtk.py and output.cxx to disable/enable them.
  • Plot the VTK files with Paraview or LLNL's Visit program.

    Bug reports

    Bug reports, comments, and suggestions are always welcome. The best channel is to create an issue on the Issue Tracker here at the original code: http://bitbucket.org/tan2/dynearthsol3d or this version with thermomechanics : https://github.com/msahamed/dynearthsol3d_energy

    License

    This program is free software: you can redistribute it and/or modify it under the terms of the MIT / X Windows System license (see the file LICENSE for the full text).

    The files under the subdirectories 3x3-C/, ann/, tetgen/, and triangles/ are distributed by their own license(s).

    •