//! Constructeur par défaut
ProblemCoreFlows();
virtual ~ProblemCoreFlows();
- // -*-*-*- Gestion du calcul -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
-
+
+ // -*-*-*- Gestion du calcul (interface ICoCo -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
/** \fn initialize
* \brief Alloue la mémoire et vérifie que le maillage et les conditions limites/initiales sont bien définis
virtual Field& getOutputField(const string& nameField )=0;//Renvoie un champs pour le postraitement
*/
+ Field getUnknownField() const;
+
//paramètres du calcul -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
/** \fn setPresentTime
* @param void
* @return the number of phases considered in the model
* */
- int getNumberOfPhases(){
+ int getNumberOfPhases() const {
return _nbPhases;
};
* \brief returns the value of the latent heat
* @param double L, the value of the latent heat
* */
- double getLatentHeat(){
+ double getLatentHeat() const{
return _latentHeat;
}
* \brief sets the saturation temperature
* @param Tsat double corresponds to saturation temperature
* */
- double getSatTemp(){
+ double getSatTemp() const {
return _Tsat;
}
* \brief returns the porosity field;
* @param
* */
- Field getPorosityField(){
+ Field getPorosityField() const {
return _porosityField;
}
* \param [in] void
* \param [out] enum NonLinearFormulation
* */
- NonLinearFormulation getNonLinearFormulation(){
+ NonLinearFormulation getNonLinearFormulation() const{
return _nonLinearFormulation;
}
* \param [in] void
* \param [out] enum SpaceScheme(upwind, centred, pressureCorrection, pressureCorrection, staggered)
* */
- SpaceScheme getSpaceScheme();
+ SpaceScheme getSpaceScheme() const;
/** \fn setNumericalScheme
* \brief sets the numerical method (upwind vs centered and explicit vs implicit)
virtual Field& getOutputField(const string& nameField )=0;//Renvoie un champs pour le postraitement
*/
+ Field getConservativeField() const ;
+ Field getPrimitiveField() const;
+
protected :
/** Number of phases in the fluid **/
int _nbPhases;
_spaceScheme=upwind;
}
-SpaceScheme ProblemFluid::getSpaceScheme()
+SpaceScheme ProblemFluid::getSpaceScheme() const
{
return _spaceScheme;
}
}
}
+Field ProblemFluid::getConservativeField() const
+{
+ if(!_initializedMemory)
+ {
+ _runLogFile->close();
+ throw CdmathException("ProblemFluid::getConservativeField call initialize() first");
+ }
+ return _UU;
+}
+
+Field ProblemFluid::getPrimitiveField() const
+{
+ if(!_initializedMemory)
+ {
+ _runLogFile->close();
+ throw CdmathException("ProblemFluid::getPrimitiveField call initialize() first");
+ }
+ return _VV;
+}
+
void ProblemFluid::terminate(){
delete[] _AroePlus;