#6004 Pore Compressibility: incorrect sign first term in computation.

This commit is contained in:
Kristian Bendiksen 2020-05-29 15:04:36 +02:00 committed by Magne Sjaastad
parent 73b02d59e6
commit ee20ec4107

View File

@ -230,7 +230,7 @@ RigFemScalarResultFrames*
if ( deltaPorePressure != 0.0 && porosity != 0.0 ) if ( deltaPorePressure != 0.0 && porosity != 0.0 )
{ {
poreCompressibility = poreCompressibility =
( biotCoefficient * evData[elmNodResIdx] ) / ( deltaPorePressure * porosity ); -( biotCoefficient * evData[elmNodResIdx] ) / ( deltaPorePressure * porosity );
// Guard against divide by zero: second term can be ignored when bulk modulus is zero, // Guard against divide by zero: second term can be ignored when bulk modulus is zero,
// which can happens when biot coefficient is 1.0 // which can happens when biot coefficient is 1.0
if ( biotCoefficient != 1.0 && porosity != 1.0 ) if ( biotCoefficient != 1.0 && porosity != 1.0 )