createFields.H File Reference
Include dependency graph for createFields.H:

Go to the source code of this file.

Functions

word continuousPhaseName (IOdictionary(IOobject("physicalProperties", runTime.constant(), mesh, IOobject::MUST_READ)).lookup("continuousPhaseName"))
 
Info<< "Reading field U\"<< endl;volVectorField Uc(IOobject(IOobject::groupName("U", continuousPhaseName), runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field p\"<< endl;volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading/calculating continuous-phase face flux field phic\"<< endl;surfaceScalarField phic(IOobject(IOobject::groupName("phi", continuousPhaseName), runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), linearInterpolate(Uc) &mesh.Sf());pressureReference pressureReference(p, pimple.dict());mesh.schemes().setFluxRequired(p.name());Info<< "Creating turbulence model\"<< endl;autoPtr< viscosityModel > continuousPhaseViscosity (viscosityModel::New(mesh))
 
dimensionedScalar rhocValue (IOobject::groupName("rho", continuousPhaseName), dimDensity, continuousPhaseViscosity->lookup(IOobject::groupName("rho", continuousPhaseName)))
 
volScalarField rhoc (IOobject(rhocValue.name(), runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, rhocValue)
 
volScalarField muc (IOobject(IOobject::groupName("mu", continuousPhaseName), runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), rhoc *continuousPhaseViscosity->nu())
 
Info<< "Creating field alphac\"<< endl;volScalarField alphac(IOobject(IOobject::groupName("alpha", continuousPhaseName), runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimless, 0));Info<< "Constructing clouds"<< endl;parcelCloudList clouds(rhoc, Uc, muc, g);scalar alphacMin(1 - mesh.solution().solverDict(alphac.name()).lookup< scalar >"max"));alphac=max(1.0 - clouds.theta(), alphacMin);alphac.correctBoundaryConditions();surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));surfaceScalarField alphaPhic(IOobject::groupName("alphaPhi", continuousPhaseName), alphacf *phic);autoPtr< phaseIncompressible::momentumTransportModel > continuousPhaseTurbulence (phaseIncompressible::momentumTransportModel::New(alphac, Uc, alphaPhic, phic, continuousPhaseViscosity))
 

Function Documentation

◆ continuousPhaseName()

word continuousPhaseName ( IOdictionary(IOobject("physicalProperties", runTime.constant(), mesh, IOobject::MUST_READ)).lookup("continuousPhaseName")  )

◆ continuousPhaseViscosity()

Info<< "Reading field U\n" << endl;volVectorField Uc( IOobject ( IOobject::groupName("U", continuousPhaseName), runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Reading field p\n" << endl;volScalarField p( IOobject ( "p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Reading/calculating continuous-phase face flux field phic\n" << endl;surfaceScalarField phic( IOobject ( IOobject::groupName("phi", continuousPhaseName), runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), linearInterpolate(Uc) & mesh.Sf());pressureReference pressureReference(p, pimple.dict());mesh.schemes().setFluxRequired(p.name());Info<< "Creating turbulence model\n" << endl;autoPtr<viscosityModel> continuousPhaseViscosity ( viscosityModel::New(mesh )

◆ rhocValue()

dimensionedScalar rhocValue ( IOobject::groupName("rho", continuousPhaseName)  ,
dimDensity  ,
continuousPhaseViscosity->  lookupIOobject::groupName("rho", continuousPhaseName) 
)

◆ rhoc()

◆ muc()

◆ continuousPhaseTurbulence()

Info<< "Creating field alphac\n" << endl;volScalarField alphac( IOobject ( IOobject::groupName("alpha", continuousPhaseName), runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), mesh, dimensionedScalar(dimless, 0));Info<< "Constructing clouds" << endl;parcelCloudList clouds(rhoc, Uc, muc, g);scalar alphacMin( 1 - mesh.solution().solverDict(alphac.name()).lookup<scalar>"max"));alphac = max(1.0 - clouds.theta(), alphacMin);alphac.correctBoundaryConditions();surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));surfaceScalarField alphaPhic( IOobject::groupName ( "alphaPhi", continuousPhaseName ), alphacf*phic);autoPtr<phaseIncompressible::momentumTransportModel> continuousPhaseTurbulence ( phaseIncompressible::momentumTransportModel::New(alphac, Uc, alphaPhic, phic, continuousPhaseViscosity )