Dear MRST Developers,
I'm working with the code, initially constructed for a two-phase (oil-water) simulation, we are numerically simulating hydraulic fracturing (with one water injecting well) - coupled flow-geomechanical simulation. I am trying to rebuild it to work on a three-phase (blackoil model); Some parts of the code are attached, I added third phase (gas) into the simulation, previous version of code was done for two phase. However, I am encountering some errors, after added gas to the simulation:
Error using cellfun
Non-scalar in Uniform output, at index 1, output 1.
Set 'UniformOutput' to false.
Error in PhysicalModel/stepFunction (line 489)
if any(cellfun(@(d) ~all(isfinite(d)), dx))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveMinistep (line 356)
model.stepFunction(state, state0, dt, drivingForces, ...
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveTimestep (line 220)
solveMinistep(solver, model, state, state0_inner, dt, drivingForces);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in MechFluidFixedStressSplitModel/stepFunction (line 107)
[mstate, mreport] = mechsolver.solveTimestep(mstate0, dt, mechModel, ...
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveMinistep (line 356)
model.stepFunction(state, state0, dt, drivingForces, ...
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveTimestep (line 220)
solveMinistep(solver, model, state, state0_inner, dt, drivingForces);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in simulateScheduleAD (line 250)
[state, report] = solver.solveTimestep(state0, dt(i), model,...
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in THREE_phase_Base (line 755)
[w, s, report] = simulateScheduleAD(initState, model, schedule, ...
Thank you in advance!
Timur Sabirov