View Issue Details

IDProjectCategoryView StatusLast Update
0004047OpenFOAMBugpublic2024-01-17 16:06
ReporterFederico Municchi Assigned Towill  
PrioritynormalSeveritycrashReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version22.04.03
Product Version11 
Fixed in Versiondev 
Summary0004047: MPPIC "moment" averaging method causes Goldschmidt to crash
DescriptionThis is a common issue with the previous versions of OpenFOAM.

Switching the averaging method to "moment" results in the MPPIC solver crashing at the first "evolve()", when the isotropy model is employed.
Steps To ReproduceJust run the GoldschmidtMPPIC tutorial with the moment averaging method.
Additional InformationERROR MESSAGE:

Solving 3-D cloud cloud

Cloud: cloud injector: model1
    Added 24750 new parcels

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 exp in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::exp(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#6 Foam::tmp<Foam::FieldField<Foam::Field, double> > Foam::exp<Foam::Field>(Foam::tmp<Foam::FieldField<Foam::Field, double> > const&) at ??:?
#7 Foam::IsotropyModels::Stochastic<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >::calculate() at ??:?
#8 void Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >::motion<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >(Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >&, Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> >::trackingData&) at ??:?
#9 void Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > >::solve<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >(Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >&, Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> >::trackingData&) at ??:?
#10 virtual thunk to Foam::ParcelCloud<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >::evolve() at ??:?
#11 Foam::parcelCloudList::evolve() at ??:?
#12 Foam::solvers::incompressibleDenseParticleFluid::prePredictor() at ??:?
#13 ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamRun"
#14 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#15 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16 ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamRun"
Floating point exception (core dumped)
TagsNo tags attached.

Activities

will

2024-01-17 16:05

manager   ~0013184

Last edited: 2024-01-17 16:06

The moment method is not reliable. It was a nice idea, but ultimately it did not solve the problem of how to construct lagrangian-to-lagrangian averages in such a way as to create a uniform distribution of parcels within the cells.

If someone had a need, and was able to fund development and maintenance work on this topic, then this issue could be revisited.

For the time being, to avoid confusion regarding the applicability of this functionality, I have removed the moment method.

https://github.com/OpenFOAM/OpenFOAM-dev/commit/a671b738c91280df42952e64a2b374b879062bc0

Issue History

Date Modified Username Field Change
2024-01-14 18:24 Federico Municchi New Issue
2024-01-17 16:05 will Assigned To => will
2024-01-17 16:05 will Status new => resolved
2024-01-17 16:05 will Resolution open => fixed
2024-01-17 16:05 will Fixed in Version => dev
2024-01-17 16:05 will Note Added: 0013184
2024-01-17 16:06 will Note Edited: 0013184