OxideTimeIntegrator.h

Go to the documentation of this file.

#ifndef included_AMP_OxideTimeIntegrator
#define included_AMP_OxideTimeIntegrator
 
 
#include "AMP/time_integrators/TimeIntegrator.h"
#include "AMP/time_integrators/oxide/OxideTimeIntegratorParameters.h"
#include "AMP/vectors/Vector.h"
 
#include <string>
 
 
namespace AMP::TimeIntegrator {
 
class OxideTimeIntegrator : public TimeIntegrator
{
public:
    explicit OxideTimeIntegrator( std::shared_ptr<TimeIntegratorParameters> parameters );
 
    virtual ~OxideTimeIntegrator();
 
    void initialize( std::shared_ptr<TimeIntegratorParameters> parameters ) override;
 
    void reset( std::shared_ptr<const TimeIntegratorParameters> parameters ) override;
 
    int advanceSolution( const double dt,
                         const bool first_step,
                         std::shared_ptr<AMP::LinearAlgebra::Vector> in,
                         std::shared_ptr<AMP::LinearAlgebra::Vector> out ) override;
 
    bool checkNewSolution( void ) override;
 
    void updateSolution( void ) override;
 
    double getNextDt( const bool good_solution ) override;
 
private:
    // declare the default constructor to be private
    OxideTimeIntegrator() = delete;
 
    // The mesh over which we define the oxide
    std::shared_ptr<AMP::Mesh::Mesh> d_mesh;
 
    // Some internal vectors
    AMP::LinearAlgebra::Vector::shared_ptr d_oxide; // Oxide depth of each point (m)
    AMP::LinearAlgebra::Vector::shared_ptr d_alpha; // Alpha depth of each point (m)
    AMP::LinearAlgebra::Vector::shared_ptr d_temp;  // Temperature of each point (K)
 
    // Internal data for calculating the oxide
    std::vector<int> N_layer;
    AMP::LinearAlgebra::Vector::shared_ptr depth; // Depth of each layer (cm)
    AMP::LinearAlgebra::Vector::shared_ptr conc;  // Oxygen concentration of each layer (g/cm^3)
};
} // namespace AMP::TimeIntegrator
 
#endif