Dear Wenqing,
To understand the initial state condition, I used following two references:
3D : ogs5-benchmarks/M/excavation/3D_Time_Controlled at master · ufz/ogs5-benchmarks · GitHub
2D(XZ Plane) : ogs5-benchmarks/THM/init at master · ufz/ogs5-benchmarks · GitHub
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>
ORCID
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\.
