Homepage Books Releases Docs Discourse

Is it bug?

Hi All,
To obtain the initial state condition for stress, the SUB_DOMAIN is used as follows. However, when I plot them using VisIt, it looks like STRESS_ZZ swapped with STRESS_YY.

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_XX

$GEO_TYPE

SUB_DOMAIN

3

0 -6005191.5+7798.95*z

1 -6005191.5+7798.95*z

2 -6005191.5+7798.95*z

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_YY

$GEO_TYPE

SUB_DOMAIN

3

0 -6005191.5+7798.95*z

1 -6005191.5+7798.95*z

2 -6005191.5+7798.95*z

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_ZZ

$GEO_TYPE

SUB_DOMAIN

3

0 -20017305+25996.5*z

1 -20017305+25996.5*z

2 -20017305+25996.5*z

Thanks in advance.

Regards,

Mohammad

KETZIN_DEFORMATION1.vtu (533 KB)

Hi Mohammad,

Just debugged the source code with your input, and found there is no bug for this sort of input. The data are distinguished by the name, e.g. STRESS_XX, STRESS_YY, STRESS_ZZ, etc by the following piece of code in ogs5

···

================================================
for (j = 0; j < (long)ic_vector.size(); j++)
{
m_ic = ic_vector[j];
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XX)
stress_ic[0] = m_ic;
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_YY)
stress_ic[1] = m_ic;
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_ZZ)
stress_ic[2] = m_ic;
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XY)
stress_ic[3] = m_ic;
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XZ)
stress_ic[4] = m_ic;
if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_YZ)
stress_ic[5] = m_ic;
}

On 06/04/2018 03:41 PM, Mohammad Islam wrote:

#INITIAL_CONDITION
$PCS_TYPE
DEFORMATION
$PRIMARY_VARIABLE
STRESS_XX
$GEO_TYPE
SUB_DOMAIN
3
0 -6005191.5+7798.95*z
1 -6005191.5+7798.95*z
2 -6005191.5+7798.95*z
#INITIAL_CONDITION
$PCS_TYPE
DEFORMATION
$PRIMARY_VARIABLE
STRESS_YY
$GEO_TYPE
SUB_DOMAIN
3
0 -6005191.5+7798.95*z
1 -6005191.5+7798.95*z
2 -6005191.5+7798.95*z
#INITIAL_CONDITION
$PCS_TYPE
DEFORMATION
$PRIMARY_VARIABLE
STRESS_ZZ
$GEO_TYPE
SUB_DOMAIN
3
0 -20017305+25996.5*z
1 -20017305+25996.5*z
2 -20017305+25996.5*z

Best regards,

Wenqing

Dear Wenqing,
To understand the initial state condition, I used following two references:

3D : https://github.com/ufz/ogs5-benchmarks/tree/master/M/excavation/3D_Time_Controlled

2D(XZ Plane) : https://github.com/ufz/ogs5-benchmarks/tree/master/THM/init

In the ViSit plot, for the 3D case, STRESS_XX, STRESS_YY & STRESS_ZZ are matching with the initial condition input file. However, the initial state conditions for the 2D(XZ Plane) are not matching. In this regard, I have attached ViSit plot files and initial condition file.

thm-decov.ic (1.03 KB)

···

On Wed, Jun 6, 2018 at 5:07 AM, Wenqing Wang wenqing.wang@ufz.de wrote:

On 06/04/2018 03:41 PM, Mohammad Islam wrote:

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_XX

$GEO_TYPE

SUB_DOMAIN

3

0 -6005191.5+7798.95*z

1 -6005191.5+7798.95*z

2 -6005191.5+7798.95*z

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_YY

$GEO_TYPE

SUB_DOMAIN

3

0 -6005191.5+7798.95*z

1 -6005191.5+7798.95*z

2 -6005191.5+7798.95*z

#INITIAL_CONDITION

$PCS_TYPE

DEFORMATION

$PRIMARY_VARIABLE

STRESS_ZZ

$GEO_TYPE

SUB_DOMAIN

3

0 -20017305+25996.5*z

1 -20017305+25996.5*z

2 -20017305+25996.5*z

You received this message because you are subscribed to a topic in the Google Groups “ogs-users” group.

To unsubscribe from this topic, visit https://groups.google.com/d/topic/ogs-users/Vfw_hD4IOCg/unsubscribe.

To unsubscribe from this group and all its topics, send an email to ogs-users+unsubscribe@googlegroups.com.

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

Hi Mohammad,

Just debugged the source code with your input, and found there is no bug for this sort of input. The data are distinguished by the name, e.g. STRESS_XX, STRESS_YY, STRESS_ZZ, etc by the following piece of code in ogs5

================================================

for (j = 0; j < (long)ic_vector.size(); j++)

{

    m_ic = ic_vector[j];

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XX)

        stress_ic[0] = m_ic;

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_YY)

        stress_ic[1] = m_ic;

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_ZZ)

        stress_ic[2] = m_ic;

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XY)

        stress_ic[3] = m_ic;

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_XZ)

        stress_ic[4] = m_ic;

    if (m_ic->getProcessPrimaryVariable() == FiniteElement::STRESS_YZ)

        stress_ic[5] = m_ic;

}

================================================Best regards,

Wenqing

-----------------------------------------------------------------

Dr. Mohammad Nurul Islam

ORISE Postdoctoral Researcher

The National Energy Technology Laboratory (NETL)

US Department of Energy (DoE)

NETL-PGH, 84/309

Email: Mohammad.Islam@netl.doe.gov

https://orcid.org/0000-0002-1798-1162

For 2D, whatever XY or XZ plane with coordinates in input, it is always converted into XY system for computation and for output.

···

On 06/06/2018 02:55 PM, Mohammad Islam wrote:

Dear Wenqing,
To understand the initial state condition, I used following two references:

3D : https://github.com/ufz/ogs5-benchmarks/tree/master/M/excavation/3D_Time_Controlled
2D(XZ Plane) : https://github.com/ufz/ogs5-benchmarks/tree/master/THM/init

In the ViSit plot, for the 3D case, STRESS_XX, STRESS_YY & STRESS_ZZ are matching with the initial condition input file. However, the initial state conditions for the 2D(XZ Plane) are not matching. In this regard, I have attached ViSit plot filesandinitial condition file.

---------------

In the *.ic file, for STRESS_YY

#INITIAL_CONDITION
$PCS_TYPE
DEFORMATION
$PRIMARY_VARIABLE
STRESS_YY
$GEO_TYPE
SUB_DOMAIN
3
0 -13243500+26487*y
1 -13243500+26487*y
2 -13243500+26487*y

It looks like for the 2D(XZ Plane) case, STRESS_ZZ swapped with STRESS_YY.

Regards,
Mohammad

*/-----------------------------------------------------------------/*

*/Dr. Mohammad Nurul Islam/*

ORISE Postdoctoral Researcher

The National Energy Technology Laboratory (NETL)

US Department of Energy (DoE)

NETL-PGH, 84/309

Email: Mohammad.Islam@netl.doe.gov <mailto:Mohammad.Islam@netl.doe.gov>

https://orcid.org/0000-0002-1798-1162

On Wed, Jun 6, 2018 at 5:07 AM, Wenqing Wang <wenqing.wang@ufz.de > <mailto:wenqing.wang@ufz.de>> wrote:

    Hi Mohammad,

    Just debugged the source code with your input, and found there is
    no bug for this sort of input. The data are distinguished by the
    name, e.g. STRESS_XX, STRESS_YY, STRESS_ZZ, etc by the following
    piece of code in ogs5
    ================================================
     for (j = 0; j < (long)ic_vector.size(); j++)
     {
     m_ic = ic_vector[j];
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_XX)
     stress_ic[0] = m_ic;
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_YY)
     stress_ic[1] = m_ic;
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_ZZ)
     stress_ic[2] = m_ic;
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_XY)
     stress_ic[3] = m_ic;
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_XZ)
     stress_ic[4] = m_ic;
     if (m_ic->getProcessPrimaryVariable() ==
    FiniteElement::STRESS_YZ)
     stress_ic[5] = m_ic;
     }
    ================================================

    On 06/04/2018 03:41 PM, Mohammad Islam wrote:

        #INITIAL_CONDITION
         $PCS_TYPE
         DEFORMATION
         $PRIMARY_VARIABLE
         STRESS_XX
         $GEO_TYPE
         SUB_DOMAIN
         3
         0 -6005191.5+7798.95*z
         1 -6005191.5+7798.95*z
         2 -6005191.5+7798.95*z
        #INITIAL_CONDITION
         $PCS_TYPE
         DEFORMATION
         $PRIMARY_VARIABLE
         STRESS_YY
         $GEO_TYPE
         SUB_DOMAIN
         3
         0 -6005191.5+7798.95*z
         1 -6005191.5+7798.95*z
         2 -6005191.5+7798.95*z
        #INITIAL_CONDITION
         $PCS_TYPE
         DEFORMATION
         $PRIMARY_VARIABLE
         STRESS_ZZ
         $GEO_TYPE
         SUB_DOMAIN
         3
         0 -20017305+25996.5*z
         1 -20017305+25996.5*z
         2 -20017305+25996.5*z

    Best regards,

    Wenqing

    -- You received this message because you are subscribed to a topic in
    the Google Groups "ogs-users" group.
    To unsubscribe from this topic, visit
    https://groups.google.com/d/topic/ogs-users/Vfw_hD4IOCg/unsubscribe
    <https://groups.google.com/d/topic/ogs-users/Vfw_hD4IOCg/unsubscribe>.
    To unsubscribe from this group and all its topics, send an email
    to ogs-users+unsubscribe@googlegroups.com
    <mailto:ogs-users%2Bunsubscribe@googlegroups.com>.
    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 ogs-users+unsubscribe@googlegroups.com <mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.

Dear Wenqing,
Thank you very much.

Regards,

Mohammad

···