Simulator should take model as an argument Harmonic Mean for transmissibility should take cell sizes into account Saturation impes solver asssembly One base class for multiphase and geomechanics Inheritance for pressure value and saturation values normalize wells by segment in phase controlled injectors Upon reading mesh node coordinate needs to be multiplied by units Field report geometry Units PressureSolver and Wells work with refined mesh https://www.dealii.org/8.4.0/doxygen/deal.II/step_46.html
Refinement at saturation fronts Initial refinement at water-oil-contact Test 3D heterogeneous bitmap Use trilinos parallel classes
Properly extract fevalues
Properly distribute into the stiffness matrix Locate wellbore on cell interfaces between two processes Use cell->point_inside() to check if wellbore is in a cell
Well segment function is either inaccurate or slow Inclined well calculation is wrong How to properly distribute flux between wells I am calculating face transmissibilities twice! Pieceman formula valid only for r << r_pieceman add check as a runtime parameter!!!!