const Foam::volSymmTensorField& Foam::constSolidThermo::directionalK() const { dimensionedSymmTensor t ( constK_.name(), constK_.dimensions(), symmTensor ( constK_.value(), 0.0, 0.0, constK_.value(), 0.0, constK_.value() ) ); return tmp<volSymmTensorField> ( new volSymmTensorField ( IOobject ( "K", mesh_.time().timeName(), mesh_, IOobject::NO_READ, IOobject::NO_WRITE ), mesh_, t ) ); }
XPP_SE::XPP_SE ( const word& name, const volScalarField& alpha, const volVectorField& U, const surfaceScalarField& phi, const dictionary& dict ) : viscoelasticLaw(name, alpha, U, phi), tau_ ( IOobject ( "tau" + name, U.time().timeName(), U.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), U.mesh() ), I_ ( dimensionedSymmTensor ( "I", dimensionSet(0, 0, 0, 0, 0, 0, 0), symmTensor ( 1, 0, 0, 1, 0, 1 ) ) ), etaS1_(dict.subDict("phase1").lookup("etaS")), etaS2_(dict.subDict("phase2").lookup("etaS")), etaP1_(dict.subDict("phase1").lookup("etaP")), etaP2_(dict.subDict("phase2").lookup("etaP")), alpha1_(dict.subDict("phase1").lookup("alpha")), alpha2_(dict.subDict("phase2").lookup("alpha")), lambdaOb1_(dict.subDict("phase1").lookup("lambdaOb")), lambdaOb2_(dict.subDict("phase2").lookup("lambdaOb")), lambdaOs1_(dict.subDict("phase1").lookup("lambdaOs")), lambdaOs2_(dict.subDict("phase2").lookup("lambdaOs")), q1_(dict.subDict("phase1").lookup("q")), q2_(dict.subDict("phase2").lookup("q")), pt_ ( dimensionedScalar ( "zero", lambdaOb1_.dimensions(), scalar( SMALL ) ) ) {}
Foam::S_MDCPP::S_MDCPP ( const word& name, const volVectorField& U, const surfaceScalarField& phi, const dictionary& dict ) : viscoelasticLaw(name, U, phi), tau_ ( IOobject ( "tau" + name, U.time().timeName(), U.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), U.mesh() ), I_ ( dimensionedSymmTensor ( "I", dimensionSet(0, 0, 0, 0, 0, 0, 0), symmTensor ( 1, 0, 0, 1, 0, 1 ) ) ), rho_(dict.lookup("rho")), etaS_(dict.lookup("etaS")), etaP_(dict.lookup("etaP")), zeta_(dict.lookup("zeta")), lambdaOb_(dict.lookup("lambdaOb")), lambdaOs_(dict.lookup("lambdaOs")), q_(dict.lookup("q")) {}
Foam::XPP_DE::XPP_DE ( const word& name, const volVectorField& U, const surfaceScalarField& phi, const dictionary& dict ) : viscoelasticLaw(name, U, phi), S_ ( IOobject ( "S" + name, U.time().timeName(), U.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), U.mesh() ), Lambda_ ( IOobject ( "Lambda" + name, U.time().timeName(), U.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), U.mesh() ), tau_ ( IOobject ( "tau" + name, U.time().timeName(), U.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), U.mesh(), dimensionedSymmTensor ( "zero", dimensionSet(1, -1, -2, 0, 0, 0, 0), symmTensor::zero ) ), I_ ( dimensionedSymmTensor ( "I", dimensionSet(0, 0, 0, 0, 0, 0, 0), symmTensor ( 1, 0, 0, 1, 0, 1 ) ) ), rho_(dict.lookup("rho")), etaS_(dict.lookup("etaS")), etaP_(dict.lookup("etaP")), alpha_(dict.lookup("alpha")), lambdaOb_(dict.lookup("lambdaOb")), lambdaOs_(dict.lookup("lambdaOs")), q_(dict.lookup("q")) {}
Foam::kineticTheoryModel::kineticTheoryModel ( const Foam::phaseModel& phasea, const Foam::volVectorField& Ub, const Foam::volScalarField& alpha, const Foam::dragModel& draga ) : phasea_(phasea), Ua_(phasea.U()), Ub_(Ub), alpha_(alpha), phia_(phasea.phi()), draga_(draga), rhoa_(phasea.rho()), da_(phasea.d()), nua_(phasea.nu()), kineticTheoryProperties_ ( IOobject ( "kineticTheoryProperties", Ua_.time().constant(), Ua_.mesh(), IOobject::MUST_READ, IOobject::NO_WRITE ) ), kineticTheory_(kineticTheoryProperties_.lookup("kineticTheory")), equilibrium_(kineticTheoryProperties_.lookup("equilibrium")), viscosityModel_ ( kineticTheoryModels::viscosityModel::New ( kineticTheoryProperties_ ) ), conductivityModel_ ( conductivityModel::New ( kineticTheoryProperties_ ) ), radialModel_ ( radialModel::New ( kineticTheoryProperties_ ) ), granularPressureModel_ ( granularPressureModel::New ( kineticTheoryProperties_ ) ), frictionalStressModel_ ( frictionalStressModel::New ( kineticTheoryProperties_ ) ), e_(kineticTheoryProperties_.lookup("e")), alphaMax_(kineticTheoryProperties_.lookup("alphaMax")), alphaMinFriction_(kineticTheoryProperties_.lookup("alphaMinFriction")), Fr_(kineticTheoryProperties_.lookup("Fr")), eta_(kineticTheoryProperties_.lookup("eta")), p_(kineticTheoryProperties_.lookup("p")), phi_(dimensionedScalar(kineticTheoryProperties_.lookup("phi"))*M_PI/180.0), Theta_ ( IOobject ( "Theta", Ua_.time().timeName(), Ua_.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), Ua_.mesh() ), mua_ ( IOobject ( "mua", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0) ), lambda_ ( IOobject ( "lambda", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::NO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0) ), pa_ ( IOobject ( "pa", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0) ), ppMagf_ ( IOobject ( "ppMagf", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0) ), kappa_ ( IOobject ( "kappa", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::NO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0) ), gs0_ ( IOobject ( "gs0", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::NO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 1.0) ), gs0Prime_ ( IOobject ( "gs0prime", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::NO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 0.0) ), //_AO_09/01/2014 // Shear stress ratio eta (here, we call it as upsilon) upsilon_ ( IOobject ( "upsilon", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 0.0) ), //_AO_09/01/2014 // Shear stress tau_ (here, we call it as upsilon) tau_ ( IOobject ( "tau", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedSymmTensor("zero", dimensionSet(1, -1, -2, 0, 0), symmTensor(0,0,0,0,0,0)) ), //_AO_09/01/2014 // Effective restitution coefficient (function of e and mu (friction coefficient) ) eEff_(0), // Friction coefficient muFric_(0), // Dilute inertial regime 0 < alpha < alphaf (p. 2) alphaf ~ 0.49 alphaf_(0), // Dense intertial regime alphac_(0), // YG 12/27/2014 alphad_(0), // Yield stress ratio upsilons_(0), // Modified kinetic theory by Chialvo-Sundaresan on/off mofidiedKineticTheoryPU_(kineticTheoryProperties_.lookup("modifiedKineticTheoryPU")), // //-AO, YG - Decompose particle pressure, Sundar's idea decomposePp_(false), //The option to use Berzi's Model in the SKT implementation Berzi_(false), paStar_ ( IOobject ( "paStar", Ua_.time().timeName(), Ua_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), Ua_.mesh(), dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0) ) { if(kineticTheoryProperties_.found("Berzi")) { Berzi_ = true; } if(kineticTheoryProperties_.found("decomposePp")) { decomposePp_ = true; } }
Foam::Leonov::Leonov ( const word& name, const volScalarField& alpha, const volVectorField& U, const surfaceScalarField& phi, const dictionary& dict ) : viscoelasticLaw(name, alpha, U, phi), sigma_ ( IOobject ( "sigma" + name, U.time().timeName(), U.mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), U.mesh() ), tau_ ( IOobject ( "tau" + name, U.time().timeName(), U.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), U.mesh(), dimensionedSymmTensor ( "zero", dimensionSet(1, -1, -2, 0, 0, 0, 0), symmTensor::zero ) ), I_ ( dimensionedSymmTensor ( "I", dimensionSet(0, 0, 0, 0, 0, 0, 0), symmTensor ( 1, 0, 0, 1, 0, 1 ) ) ), rho1_(dict.subDict("phase1").lookup("rho")), rho2_(dict.subDict("phase2").lookup("rho")), etaS1_(dict.subDict("phase1").lookup("etaS")), etaS2_(dict.subDict("phase2").lookup("etaS")), etaP1_(dict.subDict("phase1").lookup("etaP")), etaP2_(dict.subDict("phase2").lookup("etaP")), lambda1_(dict.subDict("phase1").lookup("lambda")), lambda2_(dict.subDict("phase2").lookup("lambda")) {}