diff --git a/ApplicationCode/Application/Tools/RiaEclipseUnitTools.h b/ApplicationCode/Application/Tools/RiaEclipseUnitTools.h index f18da30dc5..28f34c0cae 100644 --- a/ApplicationCode/Application/Tools/RiaEclipseUnitTools.h +++ b/ApplicationCode/Application/Tools/RiaEclipseUnitTools.h @@ -45,8 +45,10 @@ public: static double mmToMeter( double mm ) { return mm / 1000.0; } static double meterToMm( double meter ) { return 1000.0 * meter; } + static double barToPascal( double bar ) { return bar * 100000.0; } static double barToPsi( double bar ) { return bar * 14.5038; } static double barPerMeterToPsiPerFeet( double barPerMeter ) { return barPerMeter * 4.42075; } + static double gigaPascalToPascal( double gigaPascal ) { return gigaPascal * 1.0e9; } static double darcysConstant( UnitSystem unitSystem ); diff --git a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp index af03909c1b..3b00f116c8 100644 --- a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp +++ b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp @@ -18,6 +18,7 @@ #include "RigFemPartResultCalculatorPoreCompressibility.h" +#include "RiaEclipseUnitTools.h" #include "RiaLogging.h" #include "RigFemPart.h" @@ -216,8 +217,8 @@ RigFemScalarResultFrames* // Calculate bulk modulus for solids (grains). // Incoming unit for Young's Modulus is GPa: convert to Pa. - double poissonRatio = poissonRatioData[elmIdx]; - double youngsModuli = youngsModuliData[elmIdx] * 1.0e9; + double poissonRatio = poissonRatioData[elmIdx]; + double youngsModuli = RiaEclipseUnitTools::gigaPascalToPascal( youngsModuliData[elmIdx] ); double bulkModulusFrame = youngsModuli / ( 3.0 * ( 1.0 - 2.0 * poissonRatio ) ); double bulkModulus = bulkModulusFrame / ( 1.0 - biotCoefficient ); @@ -229,7 +230,8 @@ RigFemScalarResultFrames* // and convert unit from Bar to Pascal. double referencePorePressure = referencePorFrameData[nodeIdx]; double framePorePressure = porFrameData[nodeIdx]; - double deltaPorePressure = ( framePorePressure - referencePorePressure ) * 100000.0; + double deltaPorePressure = + RiaEclipseUnitTools::barToPascal( framePorePressure - referencePorePressure ); // Calculate pore compressibility double poreCompressibility = inf;