kOmega::kOmega ( const volVectorField& U, const surfaceScalarField& phi, transportModel& lamTransportModel ) : RASModel(typeName, U, phi, lamTransportModel), Cmu_ ( dimensioned<scalar>::lookupOrAddToDict ( "betaStar", coeffDict_, 0.09 ) ), beta_ ( dimensioned<scalar>::lookupOrAddToDict ( "beta", coeffDict_, 0.072 ) ), alpha_ ( dimensioned<scalar>::lookupOrAddToDict ( "alpha", coeffDict_, 0.52 ) ), alphaK_ ( dimensioned<scalar>::lookupOrAddToDict ( "alphaK", coeffDict_, 0.5 ) ), alphaOmega_ ( dimensioned<scalar>::lookupOrAddToDict ( "alphaOmega", coeffDict_, 0.5 ) ), k_ ( IOobject ( "k", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateK("k", mesh_) ), omega_ ( IOobject ( "omega", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateOmega("omega", mesh_) ), nut_ ( IOobject ( "nut", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateNut("nut", mesh_) ) { nut_ = k_/(omega_ + omegaSmall_); nut_.correctBoundaryConditions(); printCoeffs(); }
kOmega::kOmega ( const volVectorField& U, const surfaceScalarField& phi, transportModel& transport, const word& turbulenceModelName, const word& modelName ) : RASModel(modelName, U, phi, transport, turbulenceModelName), Cmu_ ( dimensioned<scalar>::lookupOrAddToDict ( "betaStar", coeffDict_, 0.09 ) ), beta_ ( dimensioned<scalar>::lookupOrAddToDict ( "beta", coeffDict_, 0.072 ) ), alpha_ ( dimensioned<scalar>::lookupOrAddToDict ( "alpha", coeffDict_, 0.52 ) ), alphaK_ ( dimensioned<scalar>::lookupOrAddToDict ( "alphaK", coeffDict_, 0.5 ) ), alphaOmega_ ( dimensioned<scalar>::lookupOrAddToDict ( "alphaOmega", coeffDict_, 0.5 ) ), k_ ( IOobject ( "k", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateK("k", mesh_) ), omega_ ( IOobject ( "omega", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateOmega("omega", mesh_) ), nut_ ( IOobject ( "nut", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::AUTO_WRITE ), autoCreateNut("nut", mesh_) ) { bound(k_, kMin_); bound(omega_, omegaMin_); nut_ = k_/omega_; nut_.correctBoundaryConditions(); printCoeffs(); }