EOS.cpp

Go to the documentation of this file.

#include "EOS.h"
#include "EOSio.h"
#include "helm.h"
#include <string>
 
namespace serif::eos {
    EOS::EOS(const EOSio& reader) : m_reader(reader) {}
    EOS::EOS(const std::string& filename, const EOSFormat format) : m_reader(EOSio(filename, format)) {}
 
    EOSOutput EOS::get(const EOSInput& in) {
        EOSOutput output;
        if (getFormat() == EOSFormat::HELM) {
            helmholtz::HELMEOSInput q;
            q.T = in.temperature; // Temperature in K
            q.rho = in.density; // Density in g/cm^3
 
            q.abar = in.composition.getMeanParticleMass(); // Mean atomic mass in g
            q.zbar = in.composition.getMeanAtomicNumber(); // Mean atomic number (dimensionless)
 
            helmholtz::HELMEOSOutput tempOutput;
            tempOutput = helmholtz::get_helm_EOS(q, *std::get<std::unique_ptr<helmholtz::HELMTable>>(m_reader.getTable()));
 
            output.electronFraction = tempOutput.ye;
            output.electronChemicalPotential = tempOutput.etaele;
            output.neutronExcessFraction = tempOutput.xnefer;
 
            // --- Pressure Variables ---
            output.pressure.total = tempOutput.ptot;
            output.pressure.gas = tempOutput.pgas;
            output.pressure.radiation = tempOutput.ptot;
            output.pressure.dDensity = tempOutput.dpresdd;
            output.pressure.dTemperature = tempOutput.dpresdt;
            output.pressure.dMeanAtomicMassNumber = tempOutput.dpresda;
            output.pressure.dMeanAtomicNumber = tempOutput.dpresdz;
 
            // --- Energy Variables ---
            output.energy.total = tempOutput.etot;
            output.energy.gas = tempOutput.egas;
            output.energy.radiation = tempOutput.erad;
            output.energy.dDensity = tempOutput.denerdd;
            output.energy.dTemperature = tempOutput.denerdt;
            output.energy.dMeanAtomicMassNumber = tempOutput.denerda;
            output.energy.dMeanAtomicNumber = tempOutput.denerdz;
 
            // --- Entropy Variables ---
            output.entropy.total = tempOutput.stot;
            output.entropy.gas = tempOutput.sgas;
            output.entropy.radiation = tempOutput.srad;
            output.entropy.dDensity = tempOutput.dentrdd;
            output.entropy.dTemperature = tempOutput.dentrdt;
            output.entropy.dMeanAtomicMassNumber = tempOutput.dentrda;
            output.entropy.dMeanAtomicNumber = tempOutput.dentrdz;
        }
        return output;
    }
 
    EOSFormat EOS::getFormat() const {
        return m_reader.getFormat();
    }
 
    const EOSio& EOS::getReader() const {
        return m_reader;
    }
}