_fluides[0]= new StiffenedGas(594.,_Pref,_Tref,1.6e6, 621.,3100.); //stiffened gas law for water at pressure 155 bar, and temperature 345°C
}
}
+
+ _fileName = "SolverlabSinglePhase";
+ PetscPrintf(PETSC_COMM_WORLD,"\n Navier-Stokes equations for single phase flow\n");
}
void SinglePhase::initialize(){
- cout<<"Initialising the Navier-Stokes model"<<endl;
- *_runLogFile<<"Initialising the Navier-Stokes model"<<endl;
+ cout<<"\n Initialising the Navier-Stokes model\n"<<endl;
+ *_runLogFile<<"\n Initialising the Navier-Stokes model\n"<<endl;
_Uroe = new double[_nVar];
_gravite = vector<double>(_nVar,0);//Not to be confused with _GravityField3d (size _Ndim). _gravite (size _Nvar) is usefull for dealing with source term and implicitation of gravity vector
}
}
void SinglePhase::save(){
+ PetscPrintf(PETSC_COMM_WORLD,"Saving numerical results at time step number %d \n\n", _nbTimeStep);
+ *_runLogFile<< "Saving numerical results at time step number "<< _nbTimeStep << endl<<endl;
+
string prim(_path+"/SinglePhasePrim_");///Results
string cons(_path+"/SinglePhaseCons_");
string allFields(_path+"/");
}
}
void SinglePhaseStaggered::initialize(){
- cout<<"Initialising the Navier-Stokes model"<<endl;
- *_runLogFile<<"Initialising the Navier-Stokes model"<<endl;
+ cout<<"\n Initialising the Navier-Stokes model\n"<<endl;
+ *_runLogFile<<"\n Initialising the Navier-Stokes model\n"<<endl;
_Uroe = new double[_nVar];
_gravite = vector<double>(_nVar,0);//Not to be confused with _GravityField3d (size _Ndim). _gravite (size _Nvar) is usefull for dealing with source term and implicitation of gravity vector
}
}
void SinglePhaseStaggered::save(){
+ PetscPrintf(PETSC_COMM_WORLD,"Saving numerical results at time step number %d \n\n", _nbTimeStep);
+ *_runLogFile<< "Saving numerical results at time step number "<< _nbTimeStep << endl<<endl;
+
string prim(_path+"/SinglePhaseStaggeredPrim_");///Results
string cons(_path+"/SinglePhaseStaggeredCons_");
prim+=_fileName;