Solver.cpp

Go to the documentation of this file.

 
// OpenCAEPoro header files
#include "Solver.hpp"
 
void Solver::Setup(Reservoir& rs, const OCPControl& ctrl)
{
    OCPModel = ctrl.GetModel();
    switch (OCPModel) {
        case ISOTHERMALMODEL:
            SetupIsoT(rs, ctrl);
            break;
        case THERMALMODEL:
            SetupT(rs, ctrl);
            break;
        default:
            OCP_ABORT("Wrong model type specified!");
            break;
    }
}
 
void Solver::SetupIsoT(Reservoir& rs, const OCPControl& ctrl)
{
    // Setup static infomation for reservoir
    rs.SetupIsoT();
    IsoTSolver.SetupMethod(rs, ctrl);
}
 
void Solver::SetupT(Reservoir& rs, const OCPControl& ctrl)
{
    rs.SetupT();
    TSolver.SetupMethod(rs, ctrl);
}
 
void Solver::InitReservoir(Reservoir& rs) const
{
    switch (OCPModel) {
        case ISOTHERMALMODEL:
            IsoTSolver.InitReservoir(rs);
            break;
        case THERMALMODEL:
            TSolver.InitReservoir(rs);
            break;
        default:
            OCP_ABORT("Wrong model type specified!");
            break;
    }
}
 
void Solver::RunSimulation(Reservoir& rs, OCPControl& ctrl, OCPOutput& output)
{
    GetWallTime timer;
    timer.Start();
    output.PrintInfoSched(rs, ctrl, timer.Stop());
    USI numTSteps = ctrl.GetNumTSteps();
    for (USI d = 0; d < numTSteps - 1; d++) {
        rs.ApplyControl(d);
        ctrl.ApplyControl(d, rs);
        while (!ctrl.IsCriticalTime(d + 1)) {
            GoOneStep(rs, ctrl);
            output.SetVal(rs, ctrl);
            if (ctrl.printLevel >= PRINT_ALL) {
                // Print Summary and critical information at every time step
                output.PrintInfo();
            }
        }
        output.PrintInfoSched(rs, ctrl, timer.Stop());
        // rs.allWells.ShowWellStatus(rs.bulk);
    }
    rs.OutInfoFinal();
    ctrl.RecordTotalTime(timer.Stop() / 1000);
}
 
void Solver::GoOneStep(Reservoir& rs, OCPControl& ctrl)
{
 
    if (ctrl.printLevel >= PRINT_SOME) {
        cout << "### DEBUG: " << setprecision(3) << fixed << ctrl.GetCurTime()
             << " Days";
        cout << ",  NR: " << ctrl.GetNRiterT() << ",  LS: " << ctrl.GetLSiterT()
             << ",  Last dt: " << ctrl.last_dt << " Days" << endl;
    }
 
    switch (OCPModel) {
        case ISOTHERMALMODEL:
            GoOneStepIsoT(rs, ctrl);
            break;
        case THERMALMODEL:
            GoOneStepT(rs, ctrl);
            break;
        default:
            OCP_ABORT("Wrong model type specified!");
            break;
    }
}
 
void Solver::GoOneStepIsoT(Reservoir& rs, OCPControl& ctrl)
{
    // Prepare for time marching
    IsoTSolver.Prepare(rs, ctrl);
 
    // Time marching with adaptive time stepsize
    while (OCP_TRUE) {
        if (ctrl.GetCurDt() < MIN_TIME_CURSTEP)
            OCP_ABORT("Time stepsize is too small!");
        // Assemble linear system
        IsoTSolver.AssembleMat(rs, ctrl);
        // Solve linear system
        IsoTSolver.SolveLinearSystem(rs, ctrl);
        if (!IsoTSolver.UpdateProperty(rs, ctrl)) {
            continue;
        }
        if (IsoTSolver.FinishNR(rs, ctrl)) break;
    }
 
    // Finish current time step
    IsoTSolver.FinishStep(rs, ctrl);
}
 
void Solver::GoOneStepT(Reservoir& rs, OCPControl& ctrl)
{
    // Prepare for time marching
    TSolver.Prepare(rs, ctrl);
 
    // Time marching with adaptive time stepsize
    while (OCP_TRUE) {
        if (ctrl.GetCurDt() < MIN_TIME_CURSTEP)
            OCP_ABORT("Time stepsize is too small!");
        // Assemble linear system
        TSolver.AssembleMat(rs, ctrl);
        // Solve linear system
        TSolver.SolveLinearSystem(rs, ctrl);
        if (!TSolver.UpdateProperty(rs, ctrl)) {
            ctrl.ResetIterNRLS();
            continue;
        }
        if (TSolver.FinishNR(rs, ctrl)) break;
    }
 
    // Finish current time step
    TSolver.FinishStep(rs, ctrl);
}
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/21/2021      Create file                          */
/*----------------------------------------------------------------------------*/