MRSTlogo_web_banner_light_thinner_Epi.png

Function omega_l.m
tutorial main page computeComposition.m computeWaterComp.m equationCompositional.m flash_calculation.m getHenryCoef.m
getR.m initStateBravo.m omega_l.m quadraticRelPerm.m runBravo.m setNonlinearSolverParameters.m setupControls.m
setupGeometry.m setupSystem.m solvefi.m vaporPressure.m

Returns compressibility coefficients

function omega = omega_l(p, system)

   [nComp, mv, cl, p_ref] = deal(system.nComp, system.fluid.mv, system.fluid.cl, ...
                                 system.fluid.p_ref);
   omega = cell(nComp + 1, 1);
   for ic = 1 : nComp + 1
      omega{ic} = mv*exp(-cl*(p - p_ref));
   end

end

Published October 8, 2014