View Issue Details

IDProjectCategoryView StatusLast Update
0001899OpenFOAMBugpublic2015-11-05 19:57
Reporterprojectionist Assigned Tohenry  
Status resolvedResolutionfixed 
PlatformGNU/LinuxOSUbuntuOS Version15.04
Product Versiondev 
Summary0001899: [reactingTwoPhaseEulerFoam] Call to aspect ratio model from Tomiyama lift model fails
DescriptionWhen I select the Tomiyama lift model in reactingTwoPhaseEulerFoam, the lift model calls the aspect ratio model via phasePair, as you can see from the stacktrace.

From the solver output, we can see, that there is an aspect ratio model selected. However, the lift model's call for an aspect ratio results in a fatal error.

Other lift models are not affected, as they do not call for the aspect ratio.
Steps To ReproduceChoose in the phaseProperties file of, say the bubbleColumn tutorial case in $FOAM_TUTORIALS/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumn, the Tomiyama lift model:

    (air in water)
        type Tomiyama;

All other settings are left as they are.
Additional InformationSelecting default blending method: linear
Selecting surfaceTensionModel for (air and water): constant
Selecting aspectRatioModel for (air in water): constant
Selecting aspectRatioModel for (water in air): constant
Selecting dragModel for (air and water): segregated
Selecting swarmCorrection for (air and water): none
Selecting dragModel for (air in water): SchillerNaumann
Selecting swarmCorrection for (air in water): none
Selecting dragModel for (water in air): SchillerNaumann
Selecting swarmCorrection for (water in air): none
Selecting virtualMassModel for (air in water): constantCoefficient
Selecting virtualMassModel for (water in air): constantCoefficient
Selecting liftModel for (air in water): Tomiyama
Selecting heatTransferModel for (air in water): RanzMarshall
Selecting heatTransferModel for (water in air): RanzMarshall
Calculating field g.h

/* other out omitted */


    request for aspectRatioModel aspectRatioModel.airInWater from objectRegistry region0 failed
    available objects of type aspectRatioModel are

    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/gerhard/OpenFOAM/OpenFOAM-dev/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::aspectRatioModel const& Foam::objectRegistry::lookupObject<Foam::aspectRatioModel>(Foam::word const&) const at ??:?
#3 Foam::orderedPhasePair::E() const at ??:?
#4 Foam::phasePair::EoH2() const at ??:?
#5 Foam::liftModels::TomiyamaLift::Cl() const at ??:?
#6 Foam::liftModel::Fi() const at ??:?
#7 Foam::liftModel::F() const at ??:?
#8 Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModel<Foam::liftModel>::F<Foam::Vector<double> >() const at ??:?
#9 Foam::MomentumTransferPhaseSystem<Foam::twoPhaseSystem>::F(Foam::phasePairKey const&) const at ??:?
#10 Foam::twoPhaseSystem::F() const at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/"
#13 ? at ??:?
Abgebrochen (Speicherabzug geschrieben)
TagsNo tags attached.



2015-11-05 14:18


aspectRatioModel.H (3,105 bytes)


2015-11-05 14:18


aspectRatioModel.C (2,293 bytes)


2015-11-05 14:20

reporter   ~0005564

The attached aspectRatioModel.H and aspectRatioModel.C should fix the issue?


2015-11-05 19:57

manager   ~0005566

@Juho: Thanks for the patch but I think it is better that the aspectRatio model is obtained from the table in phaseSystem rather than requiring it to be registered to the database; see commit 76a910ab9fb3295dee28688e1df4fd4f74c2e00f

Issue History

Date Modified Username Field Change
2015-11-05 08:43 projectionist New Issue
2015-11-05 14:18 Juho File Added: aspectRatioModel.H
2015-11-05 14:18 Juho File Added: aspectRatioModel.C
2015-11-05 14:20 Juho Note Added: 0005564
2015-11-05 19:57 henry Note Added: 0005566
2015-11-05 19:57 henry Status new => resolved
2015-11-05 19:57 henry Resolution open => fixed
2015-11-05 19:57 henry Assigned To => henry