CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Minimum PDF enthalpy exceeded during injection, without physical explanation

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 3 Post By Daveo643

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 23, 2015, 13:54
Default Minimum PDF enthalpy exceeded during injection, without physical explanation
  #1
New Member
 
Join Date: Mar 2013
Location: Canada
Posts: 22
Rep Power: 13
Daveo643 is on a distinguished road
Hello everyone,

I have searched solutions to this problem and have not found anything to help, so I hope someone would be able to address this issue anew.

I'm doing a 2D transient simulation of hydrogen injected into a chamber. The injector is modeled simply as 2 jets with a pressure inlet boundary condition. Air simultaneously enters through a separate sets of channels also with profile-defined pressure inlet boundary condition that mimics an in-cylinder compression following a polytropic.

Here is a picture of the setup.



My problem is that a some timestep point of the simulation during injection, there is a spike in the residuals - with no good physical reason occurring at that moment - and after that, I get messages "Minimum PDF enthalpy exceeded in xxx cells" and my temperature monitor shows a zone where the temperature has dipped to 173 K - an impossible phenomenon. This zone at the minimum enthalpy eventually spreads until almost the entire area of the modeled chamber.

Strange thing is, on a very similar geometry but with identical boundary and initial conditions, I don't have this problem at all!!

Here are some settings and other info that might be germane to the problem:
- partially-premixed combustion model
- Automated Grid Refinement enabled
- There is an over-expanded injection jet of hydrogen using a pressure inlet BC from 10 bar to ~1 bar, thus the PDF minimum temperature is such that during injection, this remains below the calculated static temperature.
- Coupled P-V solver
- Body force-weighted pressure discretization; second-order upwind for everything else; high-order term relaxation enabled. Second-order blending ( solve > set > numerics = 0.75)
- F-cycle multigrid with BCGSTAB stabilization for all variables.

Based on search results for solutions, I have tried:
- Gradually reducing underrelaxation as the simulation runs
- define > models > species > partially-premixed-combustion-expert > Maximum Heat Loss Fraction [0.667] = 0.333
- One suggestion is to change the polynomial coefficients for the Thermal Conductivity and Viscosity, but I have no idea where the suggested values come from and I'd rather not change these from the defaults if they work for other cases.

Any insights and assistance would be greatly appreciated!!

Last edited by Daveo643; July 23, 2015 at 15:03.
Daveo643 is offline   Reply With Quote

Old   November 2, 2015, 12:47
Default
  #2
New Member
 
Join Date: Mar 2013
Location: Canada
Posts: 22
Rep Power: 13
Daveo643 is on a distinguished road
FYI: Fixed with remeshing. There was an area poor mesh quality (that wasn't spotted because overall mesh quality statistics were excellent) where the "minimum PDF enthalpy exceeded" warning initiates and the error appears to propagate into the zone.
Far, Zari and bikooo3878 like this.
Daveo643 is offline   Reply With Quote

Reply


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
Minimum PDF table enthalpy exceeded suchy_t FLUENT 2 July 27, 2022 07:13
Minimum PDF Table Temperature jamesmaca FLUENT 15 April 8, 2022 14:30
Minimum PDF table enthalpy exceed in xxx cells charismatrix FLUENT 10 April 24, 2019 04:24
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
injection problem Mark New FLUENT 0 August 4, 2013 02:30


All times are GMT -4. The time now is 17:40.