Any way to print out material or fluid properties of elements in possibly .vtk file?

Chan-Hee_Park · March 29, 2013, 12:29am

Hi guys,

Just have a quick sit to print out material(i.e., permeability) or fluid properties (i.e., viscosity, density) of elements(or nodes) in the output file. Anybody has an example configured for this or should I code to do it?

Cheers,

Chan-Hee

Marc_Walther · March 29, 2013, 8:43am

Hej Chan-Hee,

I think, you can use Tecplot output and give Viscosity and fluid density in element out. Norbert implemented some of this output some time ago. Maybe you can use these parts of the code for permeability out…

Cheers

Marc

···

On 29 Mar 2013 01:29, “Chan-Hee Park” nbeyond@gmail.com wrote:

Hi guys,

Just have a quick sit to print out material(i.e., permeability) or fluid properties (i.e., viscosity, density) of elements(or nodes) in the output file. Anybody has an example configured for this or should I code to do it?

Cheers,

Chan-Hee

–

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+unsubscribe@googlegroups.com.

For more options, visit https://groups.google.com/groups/opt_out.

Georg_Kosakowski · March 29, 2013, 9:01am

Hoi,

permeability output in vtk legacy files is supported:

$ELE_VALUES

PERMEABILITY

Other properties need probably modifications of the source code.

Greetings, Georg

···

Am Freitag, 29. März 2013 01:29:35 UTC+1 schrieb Chan-Hee Park:

Hi guys,

Just have a quick sit to print out material(i.e., permeability) or fluid properties (i.e., viscosity, density) of elements(or nodes) in the output file. Anybody has an example configured for this or should I code to do it?

Cheers,

Chan-Hee

Chan-Hee_Park · April 1, 2013, 3:28am

Thanks Marc and Georg,

I have tried permeability output in vtk with the latest trunk. It works with VTK output as Georg mentioned. I didn’t work with VISOCITY as an ELE_VALUES though. Meanwhile, “PVD” output did not work with PERMEABILITY, either. Currently, VTK works only for PERMEABILITY by the following lines.

// PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR? // JTARON 2010

  	else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)

  		vtk_file << "VECTORS permeability double " << "\n";

  		for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)

  			MeshLib::CElem* ele = _mesh->ele_vector[j];

  			CMediumProperties* MediaProp =

  					mmp_vector[ele->GetPatchIndex()];

  			for (size_t i = 0; i < 3; i++)

  				vtk_file <<

  				MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";

  			vtk_file << "\n";

As for VISCOSITY, I guess this is a good starting point how I revise the code. First, I have to choose whether I put it as node or element value. Thanks.

···

On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:

Hoi,
permeability output in vtk legacy files is supported:

$ELE_VALUES

PERMEABILITY

Other properties need probably modifications of the source code.

Greetings, Georg

Am Freitag, 29. März 2013 01:29:35 UTC+1 schrieb Chan-Hee Park:

Hi guys,

Just have a quick sit to print out material(i.e., permeability) or fluid properties (i.e., viscosity, density) of elements(or nodes) in the output file. Anybody has an example configured for this or should I code to do it?

Cheers,

Chan-Hee

Norihiro_Watanabe2 · April 1, 2013, 10:52am

Chan-Hee,

Attached is my ad-hoc version of PVD output supporting some MMP and MFP properties as element data (evaluated at the 1st gauss point). you cannot compile it as it is because some MMP functions exist only in my branch. but you can simply remove lines to fix errors. Viscosity is not supported but you can easily add it.

best,
nori

vtk.cpp (39.4 KB)

···

On 04/01/2013 05:28 AM, Chan-Hee Park wrote:

Thanks Marc and Georg,

I have tried permeability output in vtk with the latest trunk. It works
with VTK output as Georg mentioned. I didn't work with VISOCITY as an
ELE_VALUES though. Meanwhile, "PVD" output did not work with
PERMEABILITY, either. Currently, VTK works only for PERMEABILITY by the
following lines.

    // PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR? // JTARON 2010
    else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)
    {
    vtk_file << "VECTORS permeability double " << "\n";
    for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)
    {
    MeshLib::CElem* ele = _mesh->ele_vector[j];
    CMediumProperties* MediaProp =
    mmp_vector[ele->GetPatchIndex()];
    for (size_t i = 0; i < 3; i++)
    vtk_file <<
    MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";
    vtk_file << "\n";
    }

As for VISCOSITY, I guess this is a good starting point how I revise the
code. First, I have to choose whether I put it as node or element value.
Thanks.

On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:

    Hoi,
    permeability output in vtk legacy files is supported:
    $ELE_VALUES
    PERMEABILITY

    Other properties need probably modifications of the source code.
    Greetings, Georg

    Am Freitag, 29. Mï¿½rz 2013 01:29:35 UTC+1 schrieb Chan-Hee Park:

        Hi guys,

        Just have a quick sit to print out material(i.e., permeability)
        or fluid properties (i.e., viscosity, density) of elements(or
        nodes) in the output file. Anybody has an example configured for
        this or should I code to do it?

        Cheers,

        Chan-Hee

--
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+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.

--
Dr.-Ing. Norihiro Watanabe
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum fï¿½r Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraï¿½e 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschï¿½ftsfï¿½hrer: Prof. Dr. Georg Teutsch
Administrativer Geschï¿½ftsfï¿½hrer: Dr. Heike Graï¿½mann

Chan-Hee_Park · April 1, 2013, 11:37am

Thanks Nori for your customized file. I could make the compiler happy after revising the following.

in fem_ele_std.h

public:

void Config();

in vtk.cpp

/* PCH Comment out just for compilation

case 3:

mat_value = mmp->FractureAperture(i_e);

break;

case 4:

mat_value = mmp->FractureContactAreaRatio(i_e);

break;

*/

And with Georg suggestion of permeability and PVD output option as you suggested.

in .out file

$ELE_VALUES

PERMEABILITY

That does not make it work. It prints all zero instead of 1.e-11 I have assigned. Just for your info and thanks. In fact, I prefer the PVD option as you may notice it from our conversation here in the google group.

Cheers,

Chan-Hee

···

On Monday, April 1, 2013 7:52:08 PM UTC+9, Norihiro Watanabe wrote:

Chan-Hee,

Attached is my ad-hoc version of PVD output supporting some MMP and MFP
properties as element data (evaluated at the 1st gauss point). you
cannot compile it as it is because some MMP functions exist only in my
branch. but you can simply remove lines to fix errors. Viscosity is not
supported but you can easily add it.

best,
nori

On 04/01/2013 05:28 AM, Chan-Hee Park wrote:
Thanks Marc and Georg,

I have tried permeability output in vtk with the latest trunk. It works
with VTK output as Georg mentioned. I didn’t work with VISOCITY as an
ELE_VALUES though. Meanwhile, “PVD” output did not work with
PERMEABILITY, either. Currently, VTK works only for PERMEABILITY by the
following lines.
// PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR?   // JTARON 2010
else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)
{
vtk_file << "VECTORS permeability double " << "\n";
for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)
{
MeshLib::CElem* ele = _mesh->ele_vector[j];
CMediumProperties* MediaProp =
mmp_vector[ele->GetPatchIndex()];
for (size_t i = 0; i < 3; i++)
vtk_file <<
MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";
vtk_file << "\n";
}
}
As for VISCOSITY, I guess this is a good starting point how I revise the
code. First, I have to choose whether I put it as node or element value.
Thanks.

On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:
Hoi,
permeability output in vtk legacy files is supported:
$ELE_VALUES
PERMEABILITY

Other properties need  probably modifications of the source code.
Greetings, Georg

Am Freitag, 29. Mï¿½rz 2013 01:29:35 UTC+1 schrieb Chan-Hee Park:

    Hi guys,

    Just have a quick sit to print out material(i.e., permeability)
    or fluid properties (i.e., viscosity, density) of elements(or
    nodes) in the output file. Anybody has an example configured for
    this or should I code to do it?

    Cheers,

    Chan-Hee
–
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/groups/opt_out.
–
Dr.-Ing. Norihiro Watanabe
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum fï¿½r Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraï¿½e 15 / 04318 Leipzig / Germany

norihiro...@ufz.de / http://www.ufz.de

Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschï¿½ftsfï¿½hrer: Prof. Dr. Georg Teutsch
Administrativer Geschï¿½ftsfï¿½hrer: Dr. Heike Graï¿½mann

Norihiro_Watanabe2 · April 1, 2013, 11:41am

The PVD output works differently from the VTK output. PVD checks $MMP_VALUES and $MFP_VALUES. For example, you need to add the followings lines in OUT files:

$MMP_VALUES
PERMEABILITY
$MFP_VALUES
DENSITY

best,
nori

···

On 04/01/2013 01:37 PM, Chan-Hee Park wrote:

Thanks Nori for your customized file. I could make the compiler happy
after revising the following.

in fem_ele_std.h

public:
void Config();

in vtk.cpp

/* PCH Comment out just for compilation
                 case 3:
                     mat_value = mmp->FractureAperture(i_e);
                     break;
                 case 4:
                     mat_value = mmp->FractureContactAreaRatio(i_e);
                     break;
                     */

And with Georg suggestion of permeability and PVD output option as you
suggested.

in .out file

$ELE_VALUES
PERMEABILITY

That does not make it work. It prints all zero instead of 1.e-11 I have
assigned. Just for your info and thanks. In fact, I prefer the PVD
option as you may notice it from our conversation here in the google group.

Cheers,

Chan-Hee

On Monday, April 1, 2013 7:52:08 PM UTC+9, Norihiro Watanabe wrote:

    Chan-Hee,

    Attached is my ad-hoc version of PVD output supporting some MMP and MFP
    properties as element data (evaluated at the 1st gauss point). you
    cannot compile it as it is because some MMP functions exist only in my
    branch. but you can simply remove lines to fix errors. Viscosity is not
    supported but you can easily add it.

    best,
    nori

    On 04/01/2013 05:28 AM, Chan-Hee Park wrote:
     > Thanks Marc and Georg,
     >
     > I have tried permeability output in vtk with the latest trunk. It
    works
     > with VTK output as Georg mentioned. I didn't work with VISOCITY
    as an
     > ELE_VALUES though. Meanwhile, "PVD" output did not work with
     > PERMEABILITY, either. Currently, VTK works only for PERMEABILITY
    by the
     > following lines.
     >
     > // PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR? //
    JTARON 2010
     > else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)
     > {
     > vtk_file << "VECTORS permeability double " << "\n";
     > for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)
     > {
     > MeshLib::CElem* ele = _mesh->ele_vector[j];
     > CMediumProperties* MediaProp =
     > mmp_vector[ele->GetPatchIndex()];
     > for (size_t i = 0; i < 3; i++)
     > vtk_file <<
     > MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";
     > vtk_file << "\n";
     > }
     >
     > As for VISCOSITY, I guess this is a good starting point how I
    revise the
     > code. First, I have to choose whether I put it as node or element
    value.
     > Thanks.
     >
     > On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:
     >
     > Hoi,
     > permeability output in vtk legacy files is supported:
     > $ELE_VALUES
     > PERMEABILITY
     >
     > Other properties need probably modifications of the source
    code.
     > Greetings, Georg
     >
     > Am Freitag, 29. Mï¿½rz 2013 01:29:35 UTC+1 schrieb Chan-Hee > Park:
     >
     > Hi guys,
     >
     > Just have a quick sit to print out material(i.e.,
    permeability)
     > or fluid properties (i.e., viscosity, density) of
    elements(or
     > nodes) in the output file. Anybody has an example
    configured for
     > this or should I code to do it?
     >
     > Cheers,
     >
     > Chan-Hee
     >
     > --
     > 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/groups/opt_out
    <https://groups.google.com/groups/opt_out>.
     >

    --
    Dr.-Ing. Norihiro Watanabe
    Department of Environmental Informatics (ENVINF)
    Wissenschaftler

    Helmholtz-Zentrum fï¿½r Umweltforschung GmbH - UFZ
    Helmholtz Centre for Environmental Research GmbH - UFZ
    Permoserstraï¿½e 15 / 04318 Leipzig / Germany

    norihiro...@ufz.de <javascript:> / http://www.ufz.de
    Telefon +49 341 235 1806

    Sitz der Gesellschaft: Leipzig
    Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
    Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
    Wissenschaftlicher Geschï¿½ftsfï¿½hrer: Prof. Dr. Georg Teutsch
    Administrativer Geschï¿½ftsfï¿½hrer: Dr. Heike Graï¿½mann

--
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+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.

--
Dr.-Ing. Norihiro Watanabe
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum fï¿½r Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraï¿½e 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschï¿½ftsfï¿½hrer: Prof. Dr. Georg Teutsch
Administrativer Geschï¿½ftsfï¿½hrer: Dr. Heike Graï¿½mann

Chan-Hee_Park · April 2, 2013, 4:05am

Thanks Nori. It works beautiful with the following revision.

In the .out file,

$MMP_VALUES

PERMEABILITY

$MFP_VALUES

DENSITY

VISCOSITY

$GEO_TYPE

DOMAIN

$DAT_TYPE

PVD

In vtk.cpp that Nori provided in this conversation.

struct ELEMENT_MFP_VALUES

{

static double getValue(CFluidProperties* mfp, int mfp_id) {

double mat_value = .0;

switch (mfp_id)

{

case 0:

mat_value = mfp->Density();

break;

case 1:

mat_value = mfp->Viscosity();

break;

default:

cout << “CVTK::WriteElementValues: no MFP values specified” << endl;

break;

}

return mat_value;

}

static int getMFPIndex(const std::string &mfp_name) {

int mfp_id = -1;

if (mfp_name.compare(“DENSITY”) == 0) {

mfp_id = 0;

} else if (mfp_name.compare(“VISCOSITY”) == 0) {

** mfp_id = 1;**

} else {

cout << "CVTK::WriteElementValues: no valid MFP values specified. " << mfp_name << endl;

}

return mfp_id;

and some more revision that I have made to have viscosity as a function of temperature, now I can print out viscosity change due to temperature effect. Once again, thanks everybody.

Cheers,

Chan-Hee

···

On Monday, April 1, 2013 8:41:33 PM UTC+9, Norihiro Watanabe wrote:

The PVD output works differently from the VTK output. PVD checks
$MMP_VALUES and $MFP_VALUES. For example, you need to add the followings
lines in OUT files:

$MMP_VALUES
PERMEABILITY
$MFP_VALUES
DENSITY

best,
nori

On 04/01/2013 01:37 PM, Chan-Hee Park wrote:

Thanks Nori for your customized file. I could make the compiler happy
after revising the following.

in fem_ele_std.h

public:
void Config();

in vtk.cpp

/* PCH Comment out just for compilation
case 3:
mat_value = mmp->FractureAperture(i_e);
break;
case 4:
mat_value = mmp->FractureContactAreaRatio(i_e);
break;
*/

And with Georg suggestion of permeability and PVD output option as you
suggested.

in .out file

$ELE_VALUES
PERMEABILITY

That does not make it work. It prints all zero instead of 1.e-11 I have
assigned. Just for your info and thanks. In fact, I prefer the PVD
option as you may notice it from our conversation here in the google group.

Cheers,

Chan-Hee

On Monday, April 1, 2013 7:52:08 PM UTC+9, Norihiro Watanabe wrote:

Chan-Hee,

Attached is my ad-hoc version of PVD output supporting some MMP and MFP
properties as element data (evaluated at the 1st gauss point). you
cannot compile it as it is because some MMP functions exist only in my
branch. but you can simply remove lines to fix errors. Viscosity is not
supported but you can easily add it.

best,
nori

On 04/01/2013 05:28 AM, Chan-Hee Park wrote:
 > Thanks Marc and Georg,
 >
 > I have tried permeability output in vtk with the latest trunk. It
works
 > with VTK output as Georg mentioned. I didn't work with VISOCITY
as an
 > ELE_VALUES though. Meanwhile, "PVD" output did not work with
 > PERMEABILITY, either. Currently, VTK works only for PERMEABILITY
by the
 > following lines.
 >
 >     // PRINT CHANGING (OR CONSTANT) PERMEABILITY TENSOR?   //
JTARON 2010
 >     else if (_cellArrayNames[k].compare("PERMEABILITY") == 0)
 >     {
 >     vtk_file << "VECTORS permeability double " << "\n";
 >     for (int j = 0l; j < (long) _mesh->ele_vector.size(); j++)
 >     {
 >     MeshLib::CElem* ele = _mesh->ele_vector[j];
 >     CMediumProperties* MediaProp =
 >     mmp_vector[ele->GetPatchIndex()];
 >     for (size_t i = 0; i < 3; i++)
 >     vtk_file <<
 >     MediaProp->PermeabilityTensor(j)[i * 3 + i] << " ";
 >     vtk_file << "\n";
 >     }
 >     }
 >
 >
 > As for VISCOSITY, I guess this is a good starting point how I
revise the
 > code. First, I have to choose whether I put it as node or element
value.
 > Thanks.
 >
 >
 > On Friday, March 29, 2013 6:01:29 PM UTC+9, Georg Kosakowski wrote:
 >
 >     Hoi,
 >     permeability output in vtk legacy files is supported:
 >     $ELE_VALUES
 >     PERMEABILITY
 >
 >     Other properties need  probably modifications of the source
code.
 >     Greetings, Georg
 >
 >     Am Freitag, 29. Mï¿½rz 2013 01:29:35 UTC+1 schrieb Chan-Hee > >     Park:
 >
 >         Hi guys,
 >
 >         Just have a quick sit to print out material(i.e.,
permeability)
 >         or fluid properties (i.e., viscosity, density) of
elements(or
 >         nodes) in the output file. Anybody has an example
configured for
 >         this or should I code to do it?
 >
 >         Cheers,
 >
 >         Chan-Hee
 >
 > --
 > 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/groups/opt_out](https://groups.google.com/groups/opt_out)

<[https://groups.google.com/groups/opt_out](https://groups.google.com/groups/opt_out)>.
 >
 >

--