OxideModel.h

Go to the documentation of this file.

#ifndef included_AMP_OxideModel
#define included_AMP_OxideModel
 
 
#include <string>
 
 
namespace AMP::TimeIntegrator {
 
 
class OxideModel
{
public:
    static void get_equilibrium_concetration( double T, double *C );
 
 
    static void get_diffusion_coefficients( double T, double *D );
 
 
    static int integrateOxide( double dT,
                               int N,
                               const int *N2,
                               const double *x0,
                               const double *Cb,
                               double **C0,
                               const double *D,
                               double **C1,
                               double *x1,
                               double *v1 );
 
 
private:
    static void computeVelocity( const int N,
                                 const double *x,
                                 const double *Cb,
                                 const int *Nl,
                                 double const *const *C,
                                 const double *D,
                                 double *v );
 
 
    /* This function computes the maximum timestep based on diffusion
     * @param[in] N     Number of zones
     * @param[in] x     Boundary position
     * @param[in] Cb    Concentration at boundaries
     * @param[in] C     Concentration of each layer (N)
     * @param[in] D     Diffusion coefficient
     */
    static double computeDiffustionTimestep(
        const int N, const double x[2], const double Cb[2], const double *C, const double D );
 
 
    /* This function performs a linear solve to get the future oxygen concentration
     * given the current oxygen concentration, the current and future boundary
     * velocities, and the current boundary positions.
     * @param[in] N     Number of zones
     * @param[in] dt    Timestep
     * @param[in] x0    Initial boundary position
     * @param[in] v0    Velocity of boundaries
     * @param[in] Cb    Concentration at boundaries
     * @param[in] C0    Initial concentration (N)
     * @param[in] D     Diffusion coefficient
     * @param[out] C1   Final concentration (N)
     * @param[out] x1   Final boundary position
     */
    static void solveLinearDiffusionLayer( const int N,
                                           const double dt,
                                           const double x0[2],
                                           const double v[2],
                                           const double Cb[2],
                                           const double *C0,
                                           const double D,
                                           double *C1,
                                           double *x1 );
};
} // namespace AMP::TimeIntegrator
 
#endif