Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo ( const fvMesh& mesh ) : psiThermo(mesh, word::null), twoPhaseMixture(mesh, *this), thermo1_(NULL), thermo2_(NULL) { { volScalarField T1(IOobject::groupName("T", phase1Name()), T_); T1.write(); } { volScalarField T2(IOobject::groupName("T", phase2Name()), T_); T2.write(); } thermo1_ = rhoThermo::New(mesh, phase1Name()); thermo2_ = rhoThermo::New(mesh, phase2Name()); thermo1_->validate(phase1Name(), "e"); thermo2_->validate(phase2Name(), "e"); correct(); }
Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo ( const volVectorField& U, const surfaceScalarField& phi ) : psiThermo(U.mesh(), word::null), twoPhaseMixture(U.mesh(), *this), interfaceProperties(alpha1(), U, *this), thermo1_(nullptr), thermo2_(nullptr) { { volScalarField T1(IOobject::groupName("T", phase1Name()), T_); T1.write(); } { volScalarField T2(IOobject::groupName("T", phase2Name()), T_); T2.write(); } thermo1_ = rhoThermo::New(U.mesh(), phase1Name()); thermo2_ = rhoThermo::New(U.mesh(), phase2Name()); // thermo1_->validate(phase1Name(), "e"); // thermo2_->validate(phase2Name(), "e"); correct(); }