[silviliril]

Hi everyone,

I have been trying to compile a new solver phaseChangeHeatFoam with dynamic meshing capabilities. Which I call as phaseChangeHeatDyMFoam.

Now the problem is i have given end time as 20 sec and deltaT of 1e-5. still the code ends after solving for the first iteration. Its not running in loop. What can be the possible solution?

Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicRefineFvMesh
Reading field p_rgh

Reading field T

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Creating phaseChangeTwoPhaseMixture

Selecting phaseChange model Hardt
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar

Reading g

Reading hRef
Calculating field g.h

No finite volume options present

Creating field kinetic energy K


PIMPLE: Operating solver in PISO mode

time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMGPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0 max: 0
deltaT = 1.11111e-05
Time = 1.1111111e-05

Selected 368 cells for refinement out of 4000.
Refined from 4000 to 6576 cells.
Selected 0 split points out of a possible 368.
Execution time for mesh.update() = 0.03 s
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMGPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
MULES: Solving for alpha1
Liquid phase volume fraction = 0.05  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.05  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.05  Min(alpha1) = 0  Max(alpha1) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 2.31942e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 2.17407e-10, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 1.78713e-10, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0126993, No Iterations 1
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMG:  Solving for p_rgh, Initial residual = 0.00653561, Final residual = 0.000103673, No Iterations 3
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00134984, Final residual = 9.37753e-09, No Iterations 5
time step continuity errors : sum local = 0, global = 0, cumulative = 0
solve TEqn
DILUPBiCG:  Solving for T, Initial residual = 1, Final residual = 1.05455e-13, No Iterations 2
TAve = 373.15  Min(T) = 363.15  Max(T) = 373.15
ExecutionTime = 0.15 s  ClockTime = 0 s

End

[alexvaleije]

Hi,

You should give more information, but it is very likely that you have modified the main file of the solver (i guess it's called phaseChangeHeatDyMFoam.C) in a way that you commented or deleted the time loop that makes the solver go to the beginning.

Could u post that file over here? Or at least check if you have modified the main time loop

Regards,
Alex

[silviliril]

Alejandro Valeije,

Thanks. You were right. I misplaced bracket in the time loop. I corrected it. This solved the problem.