Sphere definition and material groups association

Dear All,
I am a new user of OGS and I have a few problems getting started.

I have a 3D rectangular model with two rock layers and within the lower layer is embedded a sphere representing a magma chamber. I am doing an only mechanical simulation for the beginning. I want to simulate pressure increase inside the magma chamber and see the deformation on the overlying rocks.

So I request your assistance in the following questions if possible:

  1. How to create/define a sphere in the .GLI file? I managed to define it with GMSH and create the .MSH file but I don’t know how to define it in the .GLI file. Source terms will be applied to it so I think I need to define it in the .GLI file.

  2. How to associate different parts of the model to the material groups created during the meshing procedure?

  3. in the .MFP and .MSP files, what is the keyword/subkeyword that allows an association between the properties being defined and a particular material group (or area in the model)?

Thanks in advance!

Guillaume Doudou Faye

Department of Earth Resources Engineering

Graduate School of Engineering

Kyushu University

Japan

1 Like

1. How to create/define a sphere in the .GLI file? I managed to define it with GMSH and create the .MSH file but I don't know how to define it in the .GLI file. Source terms will be applied to it so I think I need to define it in the .GLI file.
>You can generate a triangular mesh of the sphere, and convert the mesh data to tin data, which takes the format of
    triangle index (start from 0), coordinate of node 1, coordinate of node 2, coordinate of node 3
    Here is an example (say foo.tin ):
   0 24.0424 -0.13103 -0.0393981 25 -0.308142 -0.41 25 0 0.74
   1 23.5294 0 0.718235 24.0424 -0.13103 -0.0393981 25 0 0.74
   2 23.9901 -0.128073 1.48185 23.5294 0 0.718235 25 0 0.74
   3 25 0 0.74 25 -0.308142 1.89 23.9901 -0.128073 1.48185
   4 24.0424 -0.13103 -0.0393981 23.9901 -0.128073 1.48185 25 0 0.74
   5 23.9901 -0.128073 1.48185 23.5294 0 0.718235 24.0424 -0.13103 -0.0393981
(you have to prepare a script for the conversion by yourself. )

Then define the sphere in gli as
  #SURFACE
  $NAME
   shpere
  $EPSILON
   1.0e-4
  $TIN
   fool.tin

2. How to associate different parts of the model to the material groups created during the meshing procedure?
You can define physics groups in gmsh

3. in the .MFP and .MSP files, what is the keyword/subkeyword that allows an association between the properties being defined and a particular material group (or area in the model)?
The the order of appearance of #MEDIUM_PROPERTIES (mmp) or #SOLID_PROPERTIES (msp) gives the index of material group, and the index starts from 1.
For example
GeoSys-MMP: Material Medium Properties ------------------------------------
#MEDIUM_PROPERTIES // material ID: 1
  $GEOMETRY_DIMENSION
   3
  $GEOMETRY_AREA
   1.000000e+000
  $POROSITY
   1 0.15
  $STORAGE
   1 0.
  $PERMEABILITY_TENSOR
   ANISOTROPIC 6.e-20 0 0 0 6.e-20 0 0. 0 3e-20
#MEDIUM_PROPERTIES // material ID: 2
  $GEOMETRY_DIMENSION
   3
....

···

On 03/29/2017 11:18 AM, Guillaume Doudou FAYE wrote:

Dear All,
I am a new user of OGS and I have a few problems getting started.
I have a 3D rectangular model with two rock layers and within the lower layer is embedded a sphere representing a magma chamber. I am doing an only mechanical simulation for the beginning. I want to simulate pressure increase inside the magma chamber and see the deformation on the overlying rocks.
So I request your assistance in the following questions if possible:

1. How to create/define a sphere in the .GLI file? I managed to define it with GMSH and create the .MSH file but I don't know how to define it in the .GLI file. Source terms will be applied to it so I think I need to define it in the .GLI file.
2. How to associate different parts of the model to the material groups created during the meshing procedure?
3. in the .MFP and .MSP files, what is the keyword/subkeyword that allows an association between the properties being defined and a particular material group (or area in the model)?

Thanks in advance!

Guillaume Doudou Faye

Department of Earth Resources Engineering

Graduate School of Engineering

Kyushu University

Japan

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

Hi there,

if you have a triangulated sphere as a mesh you can convert this into geometry using the Data Explorer. Just right-click the mesh and select "Convert to geometry". Save and you're done.

(The functionality for conversion is not dependent on the GUI but I don't think you can directly call this from anywhere else at the moment. If you wanted to you could write a a small programme to call it from the command line.)

Also, you can easily create all kinds of Meshes using ParaView. Among other filters there's a SphereSource-filter for creating, well, spheres. You can save it as a vtu-file and import this into the DataExplorer, convert to geometry, etc.

Cheers,
Karsten

···

Am 29.03.2017 um 11:18 schrieb Guillaume Doudou FAYE:

Dear All,
I am a new user of OGS and I have a few problems getting started.
I have a 3D rectangular model with two rock layers and within the lower
layer is embedded a sphere representing a magma chamber. I am doing an
only mechanical simulation for the beginning. I want to simulate
pressure increase inside the magma chamber and see the deformation on
the overlying rocks.
So I request your assistance in the following questions if possible:

1. How to create/define a sphere in the .GLI file? I managed to define
it with GMSH and create the .MSH file but I don't know how to define it
in the .GLI file. Source terms will be applied to it so I think I need
to define it in the .GLI file.
2. How to associate different parts of the model to the material groups
created during the meshing procedure?
3. in the .MFP and .MSP files, what is the keyword/subkeyword that
allows an association between the properties being defined and a
particular material group (or area in the model)?

Thanks in advance!

Guillaume Doudou Faye

Department of Earth Resources Engineering

Graduate School of Engineering

Kyushu University

Japan

--
You received this message because you are subscribed to the Google
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an email to [email protected]
<mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

--

Dr. Karsten Rink
Department of Environmental Informatics
Helmholtz Centre for Environmental Research - UFZ
Permoserstrasse 15, 04318 Leipzig, Germany
Tel: +49(0)-341-235-1067
e-mail: [email protected]

Thanks for the responses. I try the recommended solutions and will give feedback ASAP.
Regards,

Guillaume

···

Le mercredi 29 mars 2017 18:18:45 UTC+9, Guillaume Doudou FAYE a écrit :

Dear All,
I am a new user of OGS and I have a few problems getting started.

I have a 3D rectangular model with two rock layers and within the lower layer is embedded a sphere representing a magma chamber. I am doing an only mechanical simulation for the beginning. I want to simulate pressure increase inside the magma chamber and see the deformation on the overlying rocks.

So I request your assistance in the following questions if possible:

  1. How to create/define a sphere in the .GLI file? I managed to define it with GMSH and create the .MSH file but I don’t know how to define it in the .GLI file. Source terms will be applied to it so I think I need to define it in the .GLI file.
  1. How to associate different parts of the model to the material groups created during the meshing procedure?
  1. in the .MFP and .MSP files, what is the keyword/subkeyword that allows an association between the properties being defined and a particular material group (or area in the model)?

Thanks in advance!

Guillaume Doudou Faye

Department of Earth Resources Engineering

Graduate School of Engineering

Kyushu University

Japan

Hi Everyone,
I finished writing my input files and proceeded to run the simulation.

I am having an error message: ’ error: Error while parsing the XML file ‘file name’ at line 0: cannot open file. at ConfigTreeUtil.ccp, line 63’

I tried to run the model example thm_decov as I thought that I made a mistake in my input files, but I get the exact same error message.

Someone probably knows why this is happening, please let me know how I could solve the problem.

PS: I am using ogs 6.0.8 .

Thanks in advance!

Guillaume

Hello Guillaume,

in your earlier mails you mentioned GLI, MSP and MFP files. They are not
used with OGS version 6 but with OGS version 5. So probably you should
use OGS 5 for your work. Or am I getting something wrong?

Best regards,
Christoph

···

Am Thu, 30 Mar 2017 07:48:48 -0700 (PDT) schrieb Guillaume Doudou FAYE <[email protected]>:

Hi Everyone,
I finished writing my input files and proceeded to run the simulation.
I am having an error message: ' error: Error while parsing the XML
file 'file name' at line 0: cannot open file. at ConfigTreeUtil.ccp,
line 63' I tried to run the model example thm_decov as I thought that
I made a mistake in my input files, but I get the exact same error
message. Someone probably knows why this is happening, please let me
know how I could solve the problem.

PS: I am using ogs 6.0.8 .
Thanks in advance!

Guillaume

--
Christoph Lehmann
Computational Energy Systems
Department of Environmental Informatics

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 | 04318 Leipzig | Germany
Tel. +49 341 235 1805
E-Mail: [email protected]

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Dr. h.c. Georg Teutsch
Administrative Geschäftsführerin: Prof. Dr. Heike Graßmann

Hi Christoph,
Thank you for the hint. I used OGS 5.7.0 . The simulation could start. But it takes forever to read the .MFP file (see screenshot). When I delete the .MFP file from the directory, it works fine but stops a few seconds later without any error message in the command prompt ( see attached screenshots 1 and 2).

I have also attached a zip file containing all my input files.

Any idea of why this happens?

Just to make my idea clear, I want to model the magma chamber as an incompressible fluid and apply pressure to it as a source term. This should increase the volume and lead to the uplift and deformation of overlying rocks. Maybe I am getting something wrong in the approach?

Thanks in advance!

Best Regards,

Guillaume

Sorry I could not attach the input files. They can be seen via this link:

https://www.dropbox.com/sh/4thyw3nd6gsvmid/AABlVBDKerT226635SeJnd1wa?dl=0

Best Regards,

Guillaume

Hello Guillaume,

on first glance it does not look like there is an obvious error in the
MFP file. But since I don't use OGS 5, I cannot help very much here.
Maybe someone else on this list can?

Best regards,
Christoph

···

Am Fri, 31 Mar 2017 02:12:35 -0700 (PDT) schrieb Guillaume Doudou FAYE <[email protected]>:

Hi Christoph,
Thank you for the hint. I used OGS 5.7.0 . The simulation could
start. But it takes forever to read the .MFP file (see screenshot).
When I delete the .MFP file from the directory, it works fine but
stops a few seconds later without any error message in the command
prompt ( see attached screenshots 1 and 2).
I have also attached a zip file containing all my input files.
Any idea of why this happens?
Just to make my idea clear, I want to model the magma chamber as an
incompressible fluid and apply pressure to it as a source term. This
should increase the volume and lead to the uplift and deformation of
overlying rocks. Maybe I am getting something wrong in the approach?
Thanks in advance!

Best Regards,

Guillaume

--
Christoph Lehmann
Computational Energy Systems
Department of Environmental Informatics

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 | 04318 Leipzig | Germany
Tel. +49 341 235 1805
E-Mail: [email protected]

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Dr. h.c. Georg Teutsch
Administrative Geschäftsführerin: Prof. Dr. Heike Graßmann

I will take a look later in today.

···

On 03/31/2017 11:30 AM, Guillaume Doudou FAYE wrote:

    Sorry I could not attach the input files. They can be seen via
    this link:

Dropbox

Best Regards,

Guillaume
--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

Thanks a lot Christoph and Wenqing,

I will wait for the feedback.

Best Regards,

Guillaume

Four corrections:
1) changed all files extensions to small letters (not important).
2) removed one domain ddc file (sequential computation with DDC increase computation time)
3) added one solid property. The original file has only two solid properties, however the mesh has three material domains.
4) corrected subdomain initial condition input under keyword FUNCTION, which takes the form:
     a+b*x+c*y+d*z
   where any of the term can be ommited, e.g. just d*z

Correction 3) avoid the the program crashing.

The corrected files are in the the attached zip file. Before extracting the files, please rename the zip file extension zipp to zip.

Basement_uplift.zzip (411 KB)

···

On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

    Thanks a lot Christoph and Wenqing,

    I will wait for the feedback.
    Best Regards,
    Guillaume

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

Thank you very much for your time and help Wenqing!
Best Regards,

Guillaume

In Correction 3), I just copied a section of solid property from the existing one for testing, and you have to change it according to the definition of your example.

···

On 03/31/2017 03:22 PM, Wenqing Wang wrote:

Four corrections:
1) changed all files extensions to small letters (not important).
2) removed one domain ddc file (sequential computation with DDC increase computation time)
3) added one solid property. The original file has only two solid properties, however the mesh has three material domains.
4) corrected subdomain initial condition input under keyword FUNCTION, which takes the form:
    a+b*x+c*y+d*z
  where any of the term can be ommited, e.g. just d*z

Correction 3) avoid the the program crashing.

The corrected files are in the the attached zip file. Before extracting the files, please rename the zip file extension zipp to zip.

On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

    Thanks a lot Christoph and Wenqing,

    I will wait for the feedback.
    Best Regards,
    Guillaume

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

I have a problem there. the 3rd material group (material ID 2 in the mesh file) (the sphere representing a magma chamber) is an incompressible and viscous liquid. How should I define its properties in the solid properties files?
The Fluid properties file contains the properties for this material group.

I tried to run the simulation. It was reading the .MFP file for 30 minutes. I stopped the computation and deleted the . MFP file from the folder. Now it is running properly. But I guess the results will not be good until I properly define the material groups and properties.

Thanks in advance!

Guillaume

···

Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :

  In Correction 3), I just copied a section of solid property from the existing one for testing,  and you have to  change it according to the definition of your example.



  On 03/31/2017 03:22 PM, Wenqing Wang wrote:

Four corrections:

    1) changed all files extensions to small letters (not important).

    2) removed one domain ddc file (sequential computation with DDC increase computation time)

    3) added one solid property. The original file has only two solid properties, however the mesh has three material domains.

    4) corrected subdomain initial condition input under keyword FUNCTION, which takes the form:

        a+b*x+c*y+d*z

      where any of the term can be ommited, e.g. just d*z



    Correction 3) avoid the the program crashing.



    The corrected files are in the the attached zip file. Before extracting the files, please rename the zip file extension zipp to zip.



    On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

Thanks a lot Christoph and Wenqing,

I will wait for the feedback.

Best Regards,

Guillaume

    You received this message because you are subscribed to the Google Groups "ogs-users" group.

    To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

    For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

Did you use the files in my last email? I have tested them with the master in GitHub - ufz/ogs5: OpenGeoSys THM/C simulator version 5, and the program ran without any problem in reading files.

For the 3rd material group of an incompressible and viscous liquid, you can use treated as hard solid with a relatively large Youngs' modulus.

···

On 03/31/2017 04:28 PM, Guillaume Doudou FAYE wrote:

I have a problem there. the 3rd material group (material ID 2 in the mesh file) (the sphere representing a magma chamber) is an incompressible and viscous liquid. How should I define its properties in the solid properties files?
The Fluid properties file contains the properties for this material group.
I tried to run the simulation. It was reading the .MFP file for 30 minutes. I stopped the computation and deleted the . MFP file from the folder. Now it is running properly. But I guess the results will not be good until I properly define the material groups and properties.

Thanks in advance!

Guillaume

Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :

    In Correction 3), I just copied a section of solid property from
    the existing one for testing, and you have to change it
    according to the definition of your example.

    On 03/31/2017 03:22 PM, Wenqing Wang wrote:

    Four corrections:
    1) changed all files extensions to small letters (not important).
    2) removed one domain ddc file (sequential computation with DDC
    increase computation time)
    3) added one solid property. The original file has only two solid
    properties, however the mesh has three material domains.
    4) corrected subdomain initial condition input under keyword
    FUNCTION, which takes the form:
        a+b*x+c*y+d*z
      where any of the term can be ommited, e.g. just d*z

    Correction 3) avoid the the program crashing.

    The corrected files are in the the attached zip file. Before
    extracting the files, please rename the zip file extension zipp
    to zip.

    On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

        Thanks a lot Christoph and Wenqing,

        I will wait for the feedback.
        Best Regards,
        Guillaume

    -- You received this message because you are subscribed to the
    Google Groups "ogs-users" group.
    To unsubscribe from this group and stop receiving emails from
    it, send an email to ogs-users+...@googlegroups.com <javascript:>.
    For more options, visit https://groups.google.com/d/optout
    <https://groups.google.com/d/optout&gt;\.

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

Yes I used the files in your last email. WIthout the Fluid Properties file the simulation goes through in few minutes. But ehn the .MFP file was included it stuck at ‘MFPRead’ for nearly 30 minutes.

···

Le samedi 1 avril 2017 00:24:14 UTC+9, Wenqing Wang a écrit :

  Did you use the files in my last email? I have tested them with the master in [https://github.com/ufz/ogs5](https://github.com/ufz/ogs5)      , and the program ran  without any problem in reading files.



  For the 3rd material group of an incompressible and viscous liquid, you can use treated as hard solid with a relatively large Youngs' modulus.

  On 03/31/2017 04:28 PM, Guillaume Doudou FAYE wrote:
    I have a problem there. the 3rd material group (material ID 2 in the mesh file) (the sphere representing a magma chamber) is an incompressible and viscous liquid. How should I define its properties in the solid properties files?
      The Fluid properties file contains the properties for this material group.
      I tried to run the simulation. It was reading the .MFP file for 30 minutes. I stopped the computation and deleted the . MFP file from the folder. Now it is running properly. But I guess the results will not be good until I properly define the material groups and properties.

Thanks in advance!

Guillaume

      Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :
            In Correction 3), I just copied a section of solid property from the existing one for testing,  and you have to  change it according to the definition of your example.



            On 03/31/2017 03:22 PM, Wenqing Wang wrote:

Four corrections:

              1) changed all files extensions to small letters (not important).

              2) removed one domain ddc file (sequential computation with DDC increase computation time)

              3) added one solid property. The original file has only two solid properties, however the mesh has three material domains.

              4) corrected subdomain initial condition input under keyword FUNCTION, which takes the form:

                  a+b*x+c*y+d*z

                where any of the term can be ommited, e.g. just d*z



              Correction 3) avoid the the program crashing.



              The corrected files are in the the attached zip file. Before extracting the files, please rename the zip file extension zipp to zip.



              On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

Thanks a lot Christoph and Wenqing,

I will wait for the feedback.

Best Regards,

Guillaume

              You received this message because you are subscribed to the Google Groups "ogs-users" group.

              To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

              For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

  You received this message because you are subscribed to the Google Groups "ogs-users" group.

  To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

  For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

It's weired. Just tested on my machines again. Even with the debug model executable file, running to the start of time stepping took a couple of seconds. If you use Linux, please run dos2unix to all input files.

···

On 03/31/2017 05:47 PM, Guillaume Doudou FAYE wrote:

Yes I used the files in your last email. WIthout the Fluid Properties file the simulation goes through in few minutes. But ehn the .MFP file was included it stuck at 'MFPRead' for nearly 30 minutes.

Le samedi 1 avril 2017 00:24:14 UTC+9, Wenqing Wang a écrit :

    Did you use the files in my last email? I have tested them with
    the master in GitHub - ufz/ogs5: OpenGeoSys THM/C simulator version 5, and the program ran without any problem in reading files.

    For the 3rd material group of an incompressible and viscous
    liquid, you can use treated as hard solid with a relatively large
    Youngs' modulus.
    On 03/31/2017 04:28 PM, Guillaume Doudou FAYE wrote:

    I have a problem there. the 3rd material group (material ID 2 in
    the mesh file) (the sphere representing a magma chamber) is an
    incompressible and viscous liquid. How should I define its
    properties in the solid properties files?
    The Fluid properties file contains the properties for this
    material group.
    I tried to run the simulation. It was reading the .MFP file for
    30 minutes. I stopped the computation and deleted the . MFP file
    from the folder. Now it is running properly. But I guess the
    results will not be good until I properly define the material
    groups and properties.

    Thanks in advance!

    Guillaume

    Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :

        In Correction 3), I just copied a section of solid property
        from the existing one for testing, and you have to change
        it according to the definition of your example.

        On 03/31/2017 03:22 PM, Wenqing Wang wrote:

        Four corrections:
        1) changed all files extensions to small letters (not
        important).
        2) removed one domain ddc file (sequential computation with
        DDC increase computation time)
        3) added one solid property. The original file has only two
        solid properties, however the mesh has three material domains.
        4) corrected subdomain initial condition input under keyword
        FUNCTION, which takes the form:
            a+b*x+c*y+d*z
          where any of the term can be ommited, e.g. just d*z

        Correction 3) avoid the the program crashing.

        The corrected files are in the the attached zip file. Before
        extracting the files, please rename the zip file extension
        zipp to zip.

        On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

            Thanks a lot Christoph and Wenqing,

            I will wait for the feedback.
            Best Regards,
            Guillaume

        -- You received this message because you are subscribed to the
        Google Groups "ogs-users" group.
        To unsubscribe from this group and stop receiving emails
        from it, send an email to ogs-users+...@googlegroups.com.
        For more options, visit https://groups.google.com/d/optout
        <https://groups.google.com/d/optout&gt;\.

    -- You received this message because you are subscribed to the
    Google Groups "ogs-users" group.
    To unsubscribe from this group and stop receiving emails from it,
    send an email to ogs-users+...@googlegroups.com <javascript:>.
    For more options, visit https://groups.google.com/d/optout
    <https://groups.google.com/d/optout&gt;\.

--
You received this message because you are subscribed to the Google Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [email protected] <mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout\.

Unfortunately I am using windows. I have just retried with all input files but it is stuck again at the same step.

···

Le samedi 1 avril 2017 01:15:04 UTC+9, Wenqing Wang a écrit :

  It's weired. Just tested on my machines again. Even with the debug model executable file,  running to the start of time stepping took  a couple of seconds. If you use Linux, please run dos2unix to all input files.



  On 03/31/2017 05:47 PM, Guillaume Doudou FAYE wrote:
    Yes I used the files in your last email. WIthout the Fluid Properties file the simulation goes through in few minutes. But ehn the .MFP file was included it stuck at 'MFPRead' for nearly 30 minutes.


      Le samedi 1 avril 2017 00:24:14 UTC+9, Wenqing Wang a écrit :
            Did you use the files in my last email? I have tested them with the master in [https://github.com/ufz/ogs5](https://github.com/ufz/ogs5)                , and the program ran  without any problem in reading files.



            For the 3rd material group of an incompressible and viscous liquid, you can use treated as hard solid with a relatively large Youngs' modulus.

            On 03/31/2017 04:28 PM, Guillaume Doudou FAYE wrote:
              I have a problem there. the 3rd material group (material ID 2 in the mesh file) (the sphere representing a magma chamber) is an incompressible and viscous liquid. How should I define its properties in the solid properties files?
                The Fluid properties file contains the properties for this material group.
                I tried to run the simulation. It was reading the .MFP file for 30 minutes. I stopped the computation and deleted the . MFP file from the folder. Now it is running properly. But I guess the results will not be good until I properly define the material groups and properties.

Thanks in advance!

Guillaume

                Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :
                      In Correction 3), I just copied a section of solid property from the existing one for testing,  and you have to  change it according to the definition of your example.



                      On 03/31/2017 03:22 PM, Wenqing Wang wrote:

Four corrections:

                        1) changed all files extensions to small letters (not important).

                        2) removed one domain ddc file (sequential computation with DDC increase computation time)

                        3) added one solid property. The original file has only two solid properties, however the mesh has three material domains.

                        4) corrected subdomain initial condition input under keyword FUNCTION, which takes the form:

                            a+b*x+c*y+d*z

                          where any of the term can be ommited, e.g. just d*z



                        Correction 3) avoid the the program crashing.



                        The corrected files are in the the attached zip file. Before extracting the files, please rename the zip file extension zipp to zip.



                        On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

Thanks a lot Christoph and Wenqing,

I will wait for the feedback.

Best Regards,

Guillaume

                        You received this message because you are subscribed to the Google Groups "ogs-users" group.

                        To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

                        For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

            You received this message because you are subscribed to the Google Groups "ogs-users" group.

            To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

            For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

  You received this message because you are subscribed to the Google Groups "ogs-users" group.

  To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

  For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

maybe you should try without $COMPRESSIBILITY part in the MFP file. I
think the keyword is used only for gas phase.

···

On Sat, Apr 1, 2017 at 1:20 AM, Guillaume Doudou FAYE <[email protected]> wrote:

Unfortunately I am using windows. I have just retried with all input files
but it is stuck again at the same step.

Le samedi 1 avril 2017 01:15:04 UTC+9, Wenqing Wang a écrit :

It's weired. Just tested on my machines again. Even with the debug model
executable file, running to the start of time stepping took a couple of
seconds. If you use Linux, please run dos2unix to all input files.

On 03/31/2017 05:47 PM, Guillaume Doudou FAYE wrote:

Yes I used the files in your last email. WIthout the Fluid Properties file
the simulation goes through in few minutes. But ehn the .MFP file was
included it stuck at 'MFPRead' for nearly 30 minutes.

Le samedi 1 avril 2017 00:24:14 UTC+9, Wenqing Wang a écrit :

Did you use the files in my last email? I have tested them with the
master in GitHub - ufz/ogs5: OpenGeoSys THM/C simulator version 5, and the program ran without any
problem in reading files.

For the 3rd material group of an incompressible and viscous liquid, you
can use treated as hard solid with a relatively large Youngs' modulus.
On 03/31/2017 04:28 PM, Guillaume Doudou FAYE wrote:

I have a problem there. the 3rd material group (material ID 2 in the mesh
file) (the sphere representing a magma chamber) is an incompressible and
viscous liquid. How should I define its properties in the solid properties
files?
The Fluid properties file contains the properties for this material
group.
I tried to run the simulation. It was reading the .MFP file for 30
minutes. I stopped the computation and deleted the . MFP file from the
folder. Now it is running properly. But I guess the results will not be good
until I properly define the material groups and properties.

Thanks in advance!

Guillaume

Le vendredi 31 mars 2017 22:38:50 UTC+9, Wenqing Wang a écrit :

In Correction 3), I just copied a section of solid property from the
existing one for testing, and you have to change it according to the
definition of your example.

On 03/31/2017 03:22 PM, Wenqing Wang wrote:

Four corrections:
1) changed all files extensions to small letters (not important).
2) removed one domain ddc file (sequential computation with DDC increase
computation time)
3) added one solid property. The original file has only two solid
properties, however the mesh has three material domains.
4) corrected subdomain initial condition input under keyword FUNCTION,
which takes the form:
    a+b*x+c*y+d*z
  where any of the term can be ommited, e.g. just d*z

Correction 3) avoid the the program crashing.

The corrected files are in the the attached zip file. Before extracting
the files, please rename the zip file extension zipp to zip.

On 03/31/2017 03:02 PM, Guillaume Doudou FAYE wrote:

Thanks a lot Christoph and Wenqing,

    I will wait for the feedback.
    Best Regards,
    Guillaume

--
You received this message because you are subscribed to the Google
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an email to ogs-users+...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout\.

--
You received this message because you are subscribed to the Google Groups
"ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an
email to ogs-users+...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout\.

--
You received this message because you are subscribed to the Google Groups
"ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an
email to ogs-users+...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout\.

--
You received this message because you are subscribed to the Google Groups
"ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send an
email to [email protected].
For more options, visit https://groups.google.com/d/optout\.

--
Norihiro Watanabe