OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

store fracStress as a VoigtArray in MechContainer

Browse Source
This commit is contained in:
Arne Morten Kvarving 2025-02-05 11:30:20 +01:00
parent 7136f97d6e
commit be0fbac68b
2 changed files with 4 additions and 31 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
opm/simulators/flow/MechContainer.cpp
View File

@ -67,23 +67,11 @@ allocate(const std::size_t bufferSize,
}; };
resizeAndRegister(delstress_, "DELSTR"); resizeAndRegister(delstress_, "DELSTR");
resizeAndRegister(fracstress_, "FRCSTR");
resizeAndRegister(linstress_, "LINSTR"); resizeAndRegister(linstress_, "LINSTR");
resizeAndRegister(strain_, "STRAIN"); resizeAndRegister(strain_, "STRAIN");
resizeAndRegister(stress_, "STRESS"); resizeAndRegister(stress_, "STRESS");
this->fracstressXX_.resize(bufferSize,0.0);
rstKeywords["FRCSTRXX"] = 0;
this->fracstressYY_.resize(bufferSize,0.0);
rstKeywords["FRCSTRYY"] = 0;
this->fracstressZZ_.resize(bufferSize,0.0);
rstKeywords["FRCSTRZZ"] = 0;
this->fracstressXY_.resize(bufferSize,0.0);
rstKeywords["FRCSTRXY"] = 0;
this->fracstressXZ_.resize(bufferSize,0.0);
rstKeywords["FRCSTRXZ"] = 0;
this->fracstressYZ_.resize(bufferSize,0.0);
rstKeywords["FRCSTRYZ"] = 0;
allocated_ = true; allocated_ = true;
} }
@ -110,12 +98,7 @@ void MechContainer<Scalar>::
assignFracStress(const unsigned globalDofIdx, assignFracStress(const unsigned globalDofIdx,
const Dune::FieldVector<Scalar,6>& fracStress) const Dune::FieldVector<Scalar,6>& fracStress)
{ {
this->fracstressXX_[globalDofIdx] = fracStress[0]; this->fracstress_.assign(globalDofIdx, VoigtContainer<Scalar>(fracStress));
this->fracstressYY_[globalDofIdx] = fracStress[1];
this->fracstressZZ_[globalDofIdx] = fracStress[2];
this->fracstressYZ_[globalDofIdx] = fracStress[3];
this->fracstressXZ_[globalDofIdx] = fracStress[4];
this->fracstressXY_[globalDofIdx] = fracStress[5];
} }
template<class Scalar> template<class Scalar>
@ -179,12 +162,7 @@ outputRestart(data::Solution& sol) const
const auto solutionVectors = std::array{ const auto solutionVectors = std::array{
DataEntry{"DELSTR", UnitSystem::measure::pressure, &delstress_}, DataEntry{"DELSTR", UnitSystem::measure::pressure, &delstress_},
DataEntry{"DISP", UnitSystem::measure::length, &disp_}, DataEntry{"DISP", UnitSystem::measure::length, &disp_},
DataEntry{"FRCSTRXX", UnitSystem::measure::pressure, &fracstressXX_}, DataEntry{"FRCSTR", UnitSystem::measure::pressure, &fracstress_},
DataEntry{"FRCSTRYY", UnitSystem::measure::pressure, &fracstressYY_},
DataEntry{"FRCSTRZZ", UnitSystem::measure::pressure, &fracstressZZ_},
DataEntry{"FRCSTRXY", UnitSystem::measure::pressure, &fracstressXY_},
DataEntry{"FRCSTRXZ", UnitSystem::measure::pressure, &fracstressXZ_},
DataEntry{"FRCSTRYZ", UnitSystem::measure::pressure, &fracstressYZ_},
DataEntry{"LINSTR", UnitSystem::measure::pressure, &linstress_}, DataEntry{"LINSTR", UnitSystem::measure::pressure, &linstress_},
DataEntry{"MECHPOTF", UnitSystem::measure::pressure, &potentialForce_}, DataEntry{"MECHPOTF", UnitSystem::measure::pressure, &potentialForce_},
DataEntry{"PRESPOTF", UnitSystem::measure::pressure, &potentialPressForce_}, DataEntry{"PRESPOTF", UnitSystem::measure::pressure, &potentialPressForce_},

7
opm/simulators/flow/MechContainer.hpp
View File

@ -84,15 +84,10 @@ public:
std::array<ScalarBuffer,3> disp_; std::array<ScalarBuffer,3> disp_;
VoigtArray<Scalar> delstress_; VoigtArray<Scalar> delstress_;
VoigtArray<Scalar> fracstress_;
VoigtArray<Scalar> linstress_; VoigtArray<Scalar> linstress_;
VoigtArray<Scalar> strain_; VoigtArray<Scalar> strain_;
VoigtArray<Scalar> stress_; VoigtArray<Scalar> stress_;
ScalarBuffer fracstressXX_;
ScalarBuffer fracstressYY_;
ScalarBuffer fracstressZZ_;
ScalarBuffer fracstressXY_;
ScalarBuffer fracstressXZ_;
ScalarBuffer fracstressYZ_;
}; };
} // namespace Opm } // namespace Opm