MMGUtility

From KratosWiki
Revision as of 14:51, 8 February 2017 by Vmataix (Talk | contribs)
Jump to: navigation, search

Contents

What is MMG?

Mmg is an open source software for simplicial remeshing. MMG main page[1], and you can download the code in GitHub[2]

  • It provides 3 applications and 4 libraries:
    • The mmg2d application and the libmmg2d library: adaptation and optimization of a two-dimensional triangulation and generation of a triangulation from a set of points or from given boundary edges
    • The mmgs application and the libmmgs library: adaptation and optimization of a surface triangulation and isovalue discretization
    • The mmg3d application and the libmmg3d library: adaptation and optimization of a tetrahedral mesh and implicit domain meshing
    • The libmmg library gathering the libmmg2d, libmmgs and libmmg3d libraries
  • It uses a LGPL[3] license and it has been integrated in Kratos via the mmg_utility.h in the MeshingApplication.

Important: For use it you need to download first (look in the installation section)

    • It is used like a process, using the mmg_process.py in the MeshingApplication.

How can I install this library?

The installation of the library is fortunately quite straightforward thanks to the marvels of the configure scripts. The easiest way is the following one, because requires the minimal effort, of course if you have expertise you can find an alternative way to install the library:

  • Go to MeshingApplication/custom_external_libraries/
    • Here you can see the mmg folder, but it is empty.
    • Go to the

Once it is compiled

Go to the folder tests and run:

   python3 test_MeshingApplication.py

You should get an OK, if you don't get an OK there is something wrong:

  • Check that your compilation is the correct one
    • Follow again all the steps
    • Check that the configure.sh it is correctly compiled
  • If you get an Unexpected error:
    • This problem is because your machine is not compatible with the library, unfortunately this means that you need to wait until the problem is solved

How can I used this library?

As said before, the re-meshing is based in a process structure, in this case the 'MainKratos.py' files must be modified (you can take as reference one present in the next section). The parameters that define this process are the following ones:

       ## Settings string in json format
       default_parameters = KratosMultiphysics.Parameters("""
       {
           "output_file_name"                 : "out",
           "model_part_name"                  : "MainModelPart",
           "strategy"                         : "LevelSet",
           "level_set_strategy_parameters"              :{
               "scalar_variable"                  : "DISTANCE",
               "gradient_variable"                : "DISTANCE_GRADIENT"
           },
           "hessian_strategy_parameters"              :{
               "metric_variable"                  : ["DISTANCE"],
               "interpolation_error"              : 0.04,
               "mesh_dependent_constant"          : 0.0
           },
           "enforce_current"                  : true,
           "initial_step"                     : 1,
           "step_frequency"                   : 0,
           "automatic_remesh"                 : true,
           "automatic_remesh_parameters"      :{
               "automatic_remesh_type"            : "Ratio",
               "min_size_ratio"                   : 1.0,
               "max_size_ratio"                   : 3.0,
               "refer_type"                       : "Mean",
               "min_size_current_percentage"      : 50.0,
               "max_size_current_percentage"      : 98.0
           },
           "initial_remeshing"                : true,
           "fix_contour_model_parts"          : [],
           "minimal_size"                     : 0.1,
           "maximal_size"                     : 10.0,
           "anisotropy_remeshing"             : true,
           "anisotropy_parameters":{
               "hmin_over_hmax_anisotropic_ratio" : 0.01,
               "boundary_layer_max_distance"      : 1.0,
               "boundary_layer_min_size_ratio"    : 2.0,
               "interpolation"                    : "Linear"
           },
           "save_external_files"              : false,
           "max_number_of_searchs"            : 1000,
           "echo_level"                       : 3
       }
       """)

The meaning of each of the parameters is the following one:

Can you show us a little example?

Of course, in the File:Sphere remeshed.gid.tar.gz you can find an example that you can run in your machine. This problem consists in a continuous re-meshing using the Hessian of the velocity as reference.

The resulting output should look something similar to the following[4].

The problems presented in the tests can be used as reference too.

Personal tools
Categories