T-Junction using rhoPimpleFoam (Simulation blows and weird result)
 

Hi All,

I am trying to simulation T-Junction using rhoPimpleFoam (One Inlet, two Outlet). Mass Flow BC at both inlets (thats my requirement).

Simulation blows after some time and the results looks weird. Velocity goes to 2750 m/s. What might the problem.

I have attached the BC and problem setup, also vector plot.

Error Msg

-> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /home/umagopinath/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting


Thanks
KGN

I simulated the same case using "rhoCentralFoam" and the results converges. Attached Vector plot.

But for the same Co number the time step is significantly less in "rhoCentralFoam".

Co number 0.55

rhoCentralFoam = 3.2e-7
rhoPimpleFoam = 2.2e-06

Why such a difference?
But with rhoPimpleFoam the simulation blows after some time and also the results are weird.


I guess the boundary condition is creating the problem with "rhoPimpleFoam" solver. I know only mass inflow and mass outflow of the system plus total temperature at inlet.

Thanks
KGN

hello mecbe2002,


I am not experienced with rhoPimpleFoam or rhoPisoFoam but looking at your system I think it is over dimensioned. You can only set one Inlet and one outlet mass flow so the system can calculate the second outlet mass flow based on a fixed outlet pressure.


Best regards
Wouter

Hi wouter,

Thank for the reply.

Yes i too thinks so. But the problem is that, I know only mass flow rate/mass flux at inlet and mass flow rate at outlet.
For pressure based solver, one pressure and one velocity has to specified.

Is there is way to solve this problem?

Thanks
KGN

hello mecbe2002,
The easiest way I can think of is guessing a pressure and calculating with this to steady state and checking this with what you want and try another pressure and calculating again and so iterating till you got your solution.

hope this helps
Wouter