Nektar::VariableConverter Class Reference

#include <VariableConverter.h>

Public Member Functions
	VariableConverter (const LibUtilities::SessionReaderSharedPtr &pSession, const int spaceDim)
	~VariableConverter ()
	Destructor for VariableConverter class. More...
void	GetInternalEnergy (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &energy)
	Compute the specific internal energy \( e = (E - rho*V^2/2)/rho \). More...
void	GetEnthalpy (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &enthalpy)
	Compute the specific enthalpy \( h = e + p/rho \). More...
void	GetVelocityVector (const Array< OneD, Array< OneD, NekDouble >> &physfield, Array< OneD, Array< OneD, NekDouble >> &velocity)
	Compute the velocity field \( \mathbf{v} \) given the momentum \( \rho\mathbf{v} \). More...
void	GetMach (Array< OneD, Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &soundspeed, Array< OneD, NekDouble > &mach)
	Compute the mach number \( M = \| \mathbf{v} \|^2 / c \). More...
void	GetDynamicViscosity (const Array< OneD, const NekDouble > &temperature, Array< OneD, NekDouble > &mu)
	Compute the dynamic viscosity using the Sutherland's law \( \mu = \mu_star * (T / T_star)^3/2 * (1 + C) / (T / T_star + C) \), where: = 1.7894 * 10^-5 Kg / (m * s) T_star = 288.15 K C = 110. / 288.15. More...
void	GetAbsoluteVelocity (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &Vtot)
void	GetSensor (const MultiRegions::ExpListSharedPtr &field, const Array< OneD, const Array< OneD, NekDouble >> &physarray, Array< OneD, NekDouble > &Sensor, Array< OneD, NekDouble > &SensorKappa, int offset=1)
void	GetTemperature (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &temperature)
	Compute the temperature using the equation of state. More...
void	GetPressure (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &pressure)
	Calculate the pressure using the equation of state. More...
void	GetSoundSpeed (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &soundspeed)
	Compute the sound speed using the equation of state. More...
void	GetEntropy (const Array< OneD, const Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &entropy)
	Compute the entropy using the equation of state. More...
void	GetEFromRhoP (const Array< OneD, NekDouble > &rho, const Array< OneD, NekDouble > &pressure, Array< OneD, NekDouble > &energy)
	Compute \( e(rho,p) \) using the equation of state. More...
void	GetRhoFromPT (const Array< OneD, NekDouble > &pressure, const Array< OneD, NekDouble > &temperature, Array< OneD, NekDouble > &rho)
	Compute \( rho(p,T) \) using the equation of state. More...

Protected Attributes
LibUtilities::SessionReaderSharedPtr	m_session
EquationOfStateSharedPtr	m_eos
int	m_spacedim
NekDouble	m_pInf
NekDouble	m_rhoInf
NekDouble	m_gasConstant
NekDouble	m_mu
NekDouble	m_Skappa
NekDouble	m_Kappa

Detailed Description

Definition at line 50 of file VariableConverter.h.

Constructor & Destructor Documentation

VariableConverter()

Nektar::VariableConverter::VariableConverter	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const int	spaceDim
	)

Definition at line 44 of file VariableConverter.cpp.

References Nektar::LibUtilities::NekFactory< tKey, tBase, tParam >::CreateInstance(), Nektar::GetEquationOfStateFactory(), m_eos, m_gasConstant, m_Kappa, m_mu, m_pInf, m_rhoInf, m_session, and m_Skappa.

    : m_session(pSession), m_spacedim(spaceDim)
{
    // Create equation of state object
    std::string eosType;
    m_session->LoadSolverInfo("EquationOfState", eosType, "IdealGas");
    m_eos = GetEquationOfStateFactory().CreateInstance(eosType, m_session);

    // Parameters for dynamic viscosity
    m_session->LoadParameter("pInf", m_pInf, 101325);
    m_session->LoadParameter("rhoInf", m_rhoInf, 1.225);
    m_session->LoadParameter("GasConstant", m_gasConstant, 287.058);
    m_session->LoadParameter("mu", m_mu, 1.78e-05);

    // Parameters for sensor
    m_session->LoadParameter("Skappa", m_Skappa, -1.0);
    m_session->LoadParameter("Kappa", m_Kappa, 0.25);

}

~VariableConverter()

Nektar::VariableConverter::~VariableConverter

(

)

Destructor for VariableConverter class.

Definition at line 68 of file VariableConverter.cpp.

{
}

Member Function Documentation

GetAbsoluteVelocity()

void Nektar::VariableConverter::GetAbsoluteVelocity	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	Vtot
	)

Definition at line 195 of file VariableConverter.cpp.

References GetVelocityVector(), m_spacedim, Vmath::Vsqrt(), Vmath::Vvtvp(), and Vmath::Zero().

{
    const int nPts = physfield[0].num_elements();

    // Getting the velocity vector on the 2D normal space
    Array<OneD, Array<OneD, NekDouble>> velocity(m_spacedim);

    Vmath::Zero(Vtot.num_elements(), Vtot, 1);

    for (int i = 0; i < m_spacedim; ++i)
    {
        velocity[i] = Array<OneD, NekDouble>(nPts);
    }

    GetVelocityVector(physfield, velocity);

    for (int i = 0; i < m_spacedim; ++i)
    {
        Vmath::Vvtvp(nPts, velocity[i], 1, velocity[i], 1, Vtot, 1, Vtot, 1);
    }

    Vmath::Vsqrt(nPts, Vtot, 1, Vtot, 1);
}

GetDynamicViscosity()

void Nektar::VariableConverter::GetDynamicViscosity	(	const Array< OneD, const NekDouble > &	temperature,
		Array< OneD, NekDouble > &	mu
	)

Compute the dynamic viscosity using the Sutherland's law \( \mu = \mu_star * (T / T_star)^3/2 * (1 + C) / (T / T_star + C) \), where: = 1.7894 * 10^-5 Kg / (m * s) T_star = 288.15 K C = 110. / 288.15.

WARNING, if this routine is modified the same must be done in the FieldConvert utility ProcessWSS.cpp (this class should be restructured).

Parameters

physfield	Input physical field.
mu	The resulting dynamic viscosity.

Definition at line 179 of file VariableConverter.cpp.

References m_gasConstant, m_mu, m_pInf, and m_rhoInf.

{
    const int nPts    = temperature.num_elements();
    const NekDouble C = .38175;
    NekDouble mu_star = m_mu;
    NekDouble T_star  = m_pInf / (m_rhoInf * m_gasConstant);
    NekDouble ratio;

    for (int i = 0; i < nPts; ++i)
    {
        ratio = temperature[i] / T_star;
        mu[i] = mu_star * ratio * sqrt(ratio) * (1 + C) / (ratio + C);
    }
}

GetEFromRhoP()

void Nektar::VariableConverter::GetEFromRhoP	(	const Array< OneD, NekDouble > &	rho,
		const Array< OneD, NekDouble > &	pressure,
		Array< OneD, NekDouble > &	energy
	)

Compute \( e(rho,p) \) using the equation of state.

Parameters

rho	Input density
pressure	Input pressure
energy	The resulting internal energy.

Definition at line 405 of file VariableConverter.cpp.

References m_eos.

{
    int nPts = rho.num_elements();

    for (int i = 0; i < nPts; ++i)
    {
        energy[i] = m_eos->GetEFromRhoP(rho[i], pressure[i]);
    }
}

GetEnthalpy()

void Nektar::VariableConverter::GetEnthalpy	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	enthalpy
	)

Compute the specific enthalpy \( h = e + p/rho \).

Definition at line 102 of file VariableConverter.cpp.

References GetInternalEnergy(), GetPressure(), CG_Iterations::pressure, Vmath::Vadd(), and Vmath::Vdiv().

{
    int nPts = physfield[0].num_elements();
    Array<OneD, NekDouble> energy(nPts, 0.0);
    Array<OneD, NekDouble> pressure(nPts, 0.0);

    GetInternalEnergy(physfield, energy);
    GetPressure(physfield, pressure);

    // Calculate p/rho
    Vmath::Vdiv(nPts, pressure, 1, physfield[0], 1, enthalpy, 1);
    // Calculate h = e + p/rho
    Vmath::Vadd(nPts, energy, 1, enthalpy, 1, enthalpy, 1);
}

GetEntropy()

void Nektar::VariableConverter::GetEntropy	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	entropy
	)

Compute the entropy using the equation of state.

Parameters

physfield	Input physical field
soundspeed	The resulting sound speed \( c \).

Definition at line 383 of file VariableConverter.cpp.

References GetInternalEnergy(), and m_eos.

{
    int nPts = physfield[0].num_elements();

    Array<OneD, NekDouble> energy(nPts);
    GetInternalEnergy(physfield, energy);

    for (int i = 0; i < nPts; ++i)
    {
        entropy[i] = m_eos->GetEntropy(physfield[0][i], energy[i]);
    }
}

GetInternalEnergy()

void Nektar::VariableConverter::GetInternalEnergy	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	energy
	)

Compute the specific internal energy \( e = (E - rho*V^2/2)/rho \).

Definition at line 76 of file VariableConverter.cpp.

References m_spacedim, Vmath::Smul(), Vmath::Vdiv(), Vmath::Vsub(), and Vmath::Vvtvp().

Referenced by GetEnthalpy(), GetEntropy(), GetPressure(), GetSoundSpeed(), and GetTemperature().

{
    int nPts = physfield[0].num_elements();
    Array<OneD, NekDouble> tmp(nPts, 0.0);

    // tmp = (rho * u_i)^2
    for (int i = 0; i < m_spacedim; ++i)
    {
        Vmath::Vvtvp(nPts, physfield[i + 1], 1, physfield[i + 1], 1, tmp, 1,
                     tmp, 1);
    }
    // Divide by rho and multiply by 0.5 --> tmp = 0.5 * rho * u^2
    Vmath::Vdiv(nPts, tmp, 1, physfield[0], 1, tmp, 1);
    Vmath::Smul(nPts, 0.5, tmp, 1, tmp, 1);

    // Calculate rhoe = E - rho*V^2/2
    Vmath::Vsub(nPts, physfield[m_spacedim + 1], 1, tmp, 1, energy, 1);
    // Divide by rho
    Vmath::Vdiv(nPts, energy, 1, physfield[0], 1, energy, 1);
}

GetMach()

void Nektar::VariableConverter::GetMach	(	Array< OneD, Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	soundspeed,
		Array< OneD, NekDouble > &	mach
	)

Compute the mach number \( M = \| \mathbf{v} \|^2 / c \).

Parameters

physfield	Input physical field.
soundfield	The speed of sound corresponding to physfield.
mach	The resulting mach number \( M \).

Definition at line 145 of file VariableConverter.cpp.

References m_spacedim, Vmath::Vdiv(), Vmath::Vmul(), Vmath::Vsqrt(), and Vmath::Vvtvp().

{
    const int nPts = physfield[0].num_elements();

    Vmath::Vmul(nPts, physfield[1], 1, physfield[1], 1, mach, 1);

    for (int i = 1; i < m_spacedim; ++i)
    {
        Vmath::Vvtvp(nPts, physfield[1 + i], 1, physfield[1 + i], 1, mach, 1,
                     mach, 1);
    }

    Vmath::Vdiv(nPts, mach, 1, physfield[0], 1, mach, 1);
    Vmath::Vdiv(nPts, mach, 1, physfield[0], 1, mach, 1);
    Vmath::Vsqrt(nPts, mach, 1, mach, 1);

    Vmath::Vdiv(nPts, mach, 1, soundspeed, 1, mach, 1);
}

GetPressure()

void Nektar::VariableConverter::GetPressure	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	pressure
	)

Calculate the pressure using the equation of state.

Parameters

physfield	Input momentum.
pressure	Computed pressure field.

Definition at line 320 of file VariableConverter.cpp.

References GetInternalEnergy(), and m_eos.

Referenced by GetEnthalpy().

{
    int nPts = physfield[0].num_elements();

    Array<OneD, NekDouble> energy(nPts);
    GetInternalEnergy(physfield, energy);

    for (int i = 0; i < nPts; ++i)
    {
        pressure[i] = m_eos->GetPressure(physfield[0][i], energy[i]);
    }
}

GetRhoFromPT()

void Nektar::VariableConverter::GetRhoFromPT	(	const Array< OneD, NekDouble > &	pressure,
		const Array< OneD, NekDouble > &	temperature,
		Array< OneD, NekDouble > &	rho
	)

Compute \( rho(p,T) \) using the equation of state.

Parameters

pressure	Input pressure
temperature	Input temperature
rho	The resulting density

Definition at line 424 of file VariableConverter.cpp.

References m_eos.

{
    int nPts = pressure.num_elements();

    for (int i = 0; i < nPts; ++i)
    {
        rho[i] = m_eos->GetRhoFromPT(pressure[i], temperature[i]);
    }
}

GetSensor()

void Nektar::VariableConverter::GetSensor	(	const MultiRegions::ExpListSharedPtr &	field,
		const Array< OneD, const Array< OneD, NekDouble >> &	physarray,
		Array< OneD, NekDouble > &	Sensor,
		Array< OneD, NekDouble > &	SensorKappa,
		int	offset = 1
	)

Definition at line 221 of file VariableConverter.cpp.

References Vmath::Dot(), Vmath::Fill(), Nektar::NekConstants::kNekSqrtTol, m_Kappa, m_Skappa, and Vmath::Vsub().

{
    NekDouble Skappa;
    NekDouble order;
    Array<OneD, NekDouble> tmp;
    Array<OneD, int> expOrderElement = field->EvalBasisNumModesMaxPerExp();

    for (int e = 0; e < field->GetExpSize(); e++)
    {
        int numModesElement = expOrderElement[e];
        int nElmtPoints     = field->GetExp(e)->GetTotPoints();
        int physOffset      = field->GetPhys_Offset(e);
        int nElmtCoeffs     = field->GetExp(e)->GetNcoeffs();
        int numCutOff       = numModesElement - offset;

        if (numModesElement <= offset)
        {
            Vmath::Fill(nElmtPoints, 0.0,
                    tmp = Sensor + physOffset, 1);
            Vmath::Fill(nElmtPoints, 0.0,
                    tmp = SensorKappa + physOffset, 1);
            continue;
        }

        // create vector to save the solution points per element at P = p;
        Array<OneD, NekDouble> elmtPhys(nElmtPoints,
            tmp = physarray[0] + physOffset);
        // Compute coefficients
        Array<OneD, NekDouble> elmtCoeffs(nElmtCoeffs, 0.0);
        field->GetExp(e)->FwdTrans(elmtPhys, elmtCoeffs);

        // ReduceOrderCoeffs reduces the polynomial order of the solution
        // that is represented by the coeffs given as an inarray. This is
        // done by projecting the higher order solution onto the orthogonal
        // basis and padding the higher order coefficients with zeros.
        Array<OneD, NekDouble> reducedElmtCoeffs(nElmtCoeffs, 0.0);
        field->GetExp(e)->ReduceOrderCoeffs(numCutOff, elmtCoeffs,
                                            reducedElmtCoeffs);

        Array<OneD, NekDouble> reducedElmtPhys(nElmtPoints, 0.0);
        field->GetExp(e)->BwdTrans(reducedElmtCoeffs, reducedElmtPhys);

        NekDouble numerator   = 0.0;
        NekDouble denominator = 0.0;

        // Determining the norm of the numerator of the Sensor
        Array<OneD, NekDouble> difference(nElmtPoints, 0.0);
        Vmath::Vsub(nElmtPoints, elmtPhys, 1, reducedElmtPhys, 1, difference,
                    1);

        numerator = Vmath::Dot(nElmtPoints, difference, difference);
        denominator = Vmath::Dot(nElmtPoints, elmtPhys, elmtPhys);

        NekDouble elmtSensor = sqrt(numerator / denominator);
        elmtSensor = log10(max(elmtSensor, NekConstants::kNekSqrtTol));

        Vmath::Fill(nElmtPoints, elmtSensor, tmp = Sensor + physOffset, 1);

        // Compute reference value for sensor
        order = max(numModesElement-1, 1);
        if (order > 0 )
        {
            Skappa = m_Skappa - 4.25 * log10(static_cast<NekDouble>(order));
        }
        else
        {
            Skappa = 0.0;
        }

        // Compute artificial viscosity
        NekDouble elmtSensorKappa;
        if (elmtSensor < (Skappa-m_Kappa))
        {
            elmtSensorKappa = 0;
        }
        else if (elmtSensor > (Skappa + m_Kappa))
        {
            elmtSensorKappa = 1.0;
        }
        else
        {
            elmtSensorKappa = 0.5 *
                (1 + sin(M_PI * (elmtSensor - Skappa) / (2 * m_Kappa)));
        }
        Vmath::Fill(nElmtPoints, elmtSensorKappa,
                tmp = SensorKappa + physOffset, 1);
    }
}

GetSoundSpeed()

void Nektar::VariableConverter::GetSoundSpeed	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	soundspeed
	)

Compute the sound speed using the equation of state.

Parameters

physfield	Input physical field
soundspeed	The resulting sound speed \( c \).

Definition at line 362 of file VariableConverter.cpp.

References GetInternalEnergy(), and m_eos.

{
    int nPts = physfield[0].num_elements();

    Array<OneD, NekDouble> energy(nPts);
    GetInternalEnergy(physfield, energy);

    for (int i = 0; i < nPts; ++i)
    {
        soundspeed[i] = m_eos->GetSoundSpeed(physfield[0][i], energy[i]);
    }
}

GetTemperature()

void Nektar::VariableConverter::GetTemperature	(	const Array< OneD, const Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	temperature
	)

Compute the temperature using the equation of state.

Parameters

physfield	Input physical field.
temperature	The resulting temperature \( T \).

Definition at line 341 of file VariableConverter.cpp.

References GetInternalEnergy(), and m_eos.

{
    int nPts = physfield[0].num_elements();

    Array<OneD, NekDouble> energy(nPts);
    GetInternalEnergy(physfield, energy);

    for (int i = 0; i < nPts; ++i)
    {
        temperature[i] = m_eos->GetTemperature(physfield[0][i], energy[i]);
    }
}

GetVelocityVector()

void Nektar::VariableConverter::GetVelocityVector	(	const Array< OneD, Array< OneD, NekDouble >> &	physfield,
		Array< OneD, Array< OneD, NekDouble >> &	velocity
	)

Compute the velocity field \( \mathbf{v} \) given the momentum \( \rho\mathbf{v} \).

Parameters

physfield	Momentum field.
velocity	Velocity field.

Definition at line 126 of file VariableConverter.cpp.

References m_spacedim, and Vmath::Vdiv().

Referenced by GetAbsoluteVelocity().

{
    const int nPts = physfield[0].num_elements();

    for (int i = 0; i < m_spacedim; ++i)
    {
        Vmath::Vdiv(nPts, physfield[1 + i], 1, physfield[0], 1, velocity[i], 1);
    }
}

Member Data Documentation

m_eos

EquationOfStateSharedPtr Nektar::VariableConverter::m_eos

protected

Definition at line 99 of file VariableConverter.h.

Referenced by GetEFromRhoP(), GetEntropy(), GetPressure(), GetRhoFromPT(), GetSoundSpeed(), GetTemperature(), and VariableConverter().

m_gasConstant

NekDouble Nektar::VariableConverter::m_gasConstant

protected

Definition at line 103 of file VariableConverter.h.

Referenced by GetDynamicViscosity(), and VariableConverter().

m_Kappa

NekDouble Nektar::VariableConverter::m_Kappa

protected

Definition at line 106 of file VariableConverter.h.

Referenced by GetSensor(), and VariableConverter().

m_mu

NekDouble Nektar::VariableConverter::m_mu

protected

Definition at line 104 of file VariableConverter.h.

Referenced by GetDynamicViscosity(), and VariableConverter().

m_pInf

NekDouble Nektar::VariableConverter::m_pInf

protected

Definition at line 101 of file VariableConverter.h.

Referenced by GetDynamicViscosity(), and VariableConverter().

m_rhoInf

NekDouble Nektar::VariableConverter::m_rhoInf

protected

Definition at line 102 of file VariableConverter.h.

Referenced by GetDynamicViscosity(), and VariableConverter().

m_session

LibUtilities::SessionReaderSharedPtr Nektar::VariableConverter::m_session

protected

Definition at line 98 of file VariableConverter.h.

Referenced by VariableConverter().

m_Skappa

NekDouble Nektar::VariableConverter::m_Skappa

protected

Definition at line 105 of file VariableConverter.h.

Referenced by GetSensor(), and VariableConverter().

m_spacedim

int Nektar::VariableConverter::m_spacedim

protected

Definition at line 100 of file VariableConverter.h.

Referenced by GetAbsoluteVelocity(), GetInternalEnergy(), GetMach(), and GetVelocityVector().