CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

[OpenFOAM v8] Why does Lun Savage radial function fail to converge?

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By edwardKGN

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 8, 2021, 01:23
Default [OpenFOAM v8] Why does Lun Savage radial function fail to converge?
  #1
New Member
 
Edward Ng Kian Guan
Join Date: Feb 2021
Posts: 5
Rep Power: 5
edwardKGN is on a distinguished road
I was testing usage of different radial distribution functions by modifying the tutorial case "fluidisedBed".

Originally the tutorial case used radial distribution of Sinclair Jackson.

The modified tutorial case uses radial distribution of Lun Savage, and had failed to converge at t = 1.41, yielding the error log shown at the bottom of this post.

From my understanding, the only difference between Sinclair Jackson's function and Lun Savage function, is the power of the function.

I am uncertain why did the simulation failed to converge.

The modified case file is at,https://github.com/edwardKGN/fluidis...age_radialDist.

Thank you for reading and helping!

/start of error log

Courant Number mean: 0.198841 max: 0.534878
deltaT = 0.002
Time = 1.408

PIMPLE: Iteration 1
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.105292, Final residual = 3.13208e-10, No Iterations 3
particles fraction, min, max = 0.274977 5.87739e-06 0.578727
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.673921, Final residual = 7.37451e-10, No Iterations 3
particles fraction, min, max = 0.274977 5.88012e-06 0.579236
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.740731, Final residual = 3.65256e-10, No Iterations 3
particles fraction, min, max = 0.274977 5.88293e-06 0.581966
Constructing momentum equations
particles min/max T 600 - 600
water min/max T 300 - 600
GAMG: Solving for p_rgh, Initial residual = 0.0370966, Final residual = 0.000304703, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.00313407, Final residual = 4.84552e-09, No Iterations 16
PIMPLE: Iteration 2
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0344012, Final residual = 6.5296e-11, No Iterations 5
particles fraction, min, max = 0.274977 5.8776e-06 0.577119
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0762682, Final residual = 3.26397e-10, No Iterations 5
particles fraction, min, max = 0.274977 5.88029e-06 0.578326
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.141742, Final residual = 1.92945e-10, No Iterations 5
particles fraction, min, max = 0.274977 5.88307e-06 0.580262
Constructing momentum equations
particles min/max T 600 - 600
water min/max T 300 - 600
GAMG: Solving for p_rgh, Initial residual = 0.0274588, Final residual = 0.000266116, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.00256743, Final residual = 8.92578e-09, No Iterations 15
PIMPLE: Iteration 3
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0729752, Final residual = 1.91868e-10, No Iterations 3
particles fraction, min, max = 0.274977 5.87765e-06 0.580288
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.727065, Final residual = 2.14024e-11, No Iterations 3
particles fraction, min, max = 0.274977 5.8804e-06 0.584208
MULES: Solving for alpha.particles
DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.772823, Final residual = 4.08709e-10, No Iterations 3
particles fraction, min, max = 0.274977 5.88322e-06 0.586477
Constructing momentum equations
particles min/max T 600 - 600
water min/max T 300 - 600
GAMG: Solving for p_rgh, Initial residual = 0.115123, Final residual = 0.00058304, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00681858, Final residual = 5.59614e-09, No Iterations 17
DILUPBiCGStab: Solving for Theta.particles, Initial residual = 0.855074, Final residual = 5.41886e-07, No Iterations 1
ExecutionTime = 149.94 s


Courant Number mean: 0.201101 max: 0.535999
deltaT = 0.002
Time = 1.41

PIMPLE: Iteration 1
MULES: Solving for alpha.particles
#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 powf64 in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
#6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at ??:?
#7 Foam::kineticTheoryModels::radialModels::LunSavage ::g0prime(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&, Foam::dimensioned<double> const&) const at ??:?
#8 Foam::RASModels::kineticTheoryModel::pPrime() const at ??:?
#9 Foam::MovingPhaseModel<Foam::ThermoPhaseModel<Foam ::phaseModel, Foam::rhoThermo> >::pPrime() const at ??:?
#10 Foam::MomentumTransferPhaseSystem<Foam::phaseSyste m>::DByAfs(Foam::PtrList<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) const at ??:?
#11 Foam::phaseSystem::solve(Foam::PtrList<Foam::Geome tricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) at ??:?
#12 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

/end of error log
mabrouka likes this.
edwardKGN is offline   Reply With Quote

Old   July 15, 2022, 09:56
Default Lun
  #2
New Member
 
Abdo
Join Date: Jan 2022
Location: Canada
Posts: 4
Rep Power: 4
mabrouka is on a distinguished road
Hi there,

Have you solved your problem, because I have same error
mabrouka is offline   Reply With Quote

Reply

Tags
multiphase, radial distribution

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 14:11
A girl fail to plot velocity profile when mesh changes + Wall function asherah STAR-CCM+ 0 February 19, 2010 17:45
Compilation errors in ThirdPartymallochoard feng_w OpenFOAM Installation 1 January 25, 2009 06:59
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
Droplet Evaporation Christian Main CFD Forum 2 February 27, 2007 06:27


All times are GMT -4. The time now is 23:49.