shells/applications/cylinderElastic/helpers/DetonationPressure1_0.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2004 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef DETONATIONPRESSURE_H
#define DETONATIONPRESSURE_H
#include <cmath>


namespace cylExp {
    struct DetonationPressure;
    struct PressurePulse;
}


struct DetonationPressure {
    
    double operator()(double time, double x, double y, double z) const 
        {
        // Beltman - Shepherd paper eq. A.7

        const double p1 = 1.0e5;
        const double p2 = 3.33e6*0.95;
        const double p3 = 1.165e6;
        const double vcj = 2375.0;
        
        const double zbegin = 1.690;
        
        time += (zbegin/vcj);
        z += zbegin;

        double tcj = (z)/vcj;
        double T = tcj/3.0;

        const double initialPressure = p1;

        if (time>tcj) {
            double pressure = (p2-p3)*std::exp((tcj-time)/T)+p3;
            return pressure-initialPressure;
        } else {
            return p1-initialPressure;
        }

        return 0.0; 
    }   
};

#endif

---