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

K and Epsilon Do Not Change

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 17, 2013, 12:54
Default K and Epsilon Do Not Change
  #1
New Member
 
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 4
mosquera is on a distinguished road
Greetings foamers,

I've been working on a 2d case and I feel making good progress improving my mesh, but I've run into a fairly substantial issue.

I'm using a realizable k-epsilon turbulence model and inputting experimental data as a boundary condition on one side of my 2d channel.

The solution converges, but in the solution k and epsilon both are their proper values on the BC inlet previously described and drop immediately to whatever I defined as the internal field in the /0 directory and remain that way.

I have a case I've been developing in 3d and it (though not yet converged sufficiently) it looks to have the same issue.

Does anyone have any experience with an issue like this?

Many thanks,
Matt
mosquera is offline   Reply With Quote

Old   July 18, 2013, 03:21
Default
  #2
New Member
 
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 4
mosquera is on a distinguished road
Some more information on the case:

I'm using timeVaryingMappedFixedValue as the inlet boundary condition for the channel. I have experimental data for the velocity at the inlet and use empirical relationships to find k and epsilon at these points.

Here's a schematic diagram of my case and boundary conditions (2D):

____B____
| |
| | <--C
A|_____ ___|
D

A: Inlet Boundary Condition
B: End of simulation wall
C: Outlet Boundary Condition
D: Test Chamber Wall

Boundary Conditions:
A:
U: timeVaryingMappedFixedValue
k: timeVaryingMappedFixedValue
epsilon: timeVaryingMappedFixedValue
p: zeroGradient
B:
U: symmetryPlane
k: symmetryPlane
epsilon: symmetryPlane
p: symmetryPlane
C:
U: zeroGradient
k: zeroGradient
epsilon: zeroGradient
p: fixedValue (0)
D:
U: fixedValue (0)
k: kLowReWallFunction
epsilon: epsilonLowReWallFunction
p: zeroGradient

Additionally, I checked and the solution converges more quickly using upwind instead of upwind linear, but the problem persists.

Thank you for any advice/tips on things to try,
Matt
mosquera is offline   Reply With Quote

Old   July 18, 2013, 17:47
Default
  #3
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
Hi Matt,

Which solver are you using?
Could you post the log-output of a few time steps?

Cheers,

L
Lieven is offline   Reply With Quote

Old   July 22, 2013, 10:24
Default
  #4
New Member
 
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 4
mosquera is on a distinguished road
Hi Lieven,

I'm using simpleFoam.

I still have solution instability, but I realized the reason for the zeros for k and epsilon was because of an error calculating them at the inlet.

I'm using timevaryingmappedfixedvalues to map experimental values to the inlet of my flow. I realized the problem is that I'm using empirical correlations to find k and epsilon at the inlet (I only have hotwire data). What's happening is that they're not exactly correct and so the solution is unstable in that region.

Do you have any experience with this boundary condition? Is there a way to "smooth" the boundary condition and make the solution more stable? I hope to expand this simulation to three dimensions, so a more stable solution would be incredibly helpful.
mosquera is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37
Calculation of k and epsilon freezes Nigirim OpenFOAM Running, Solving & CFD 1 November 14, 2012 08:52
epsilon and K blowing up. sivakumar OpenFOAM Running, Solving & CFD 1 October 25, 2012 04:50
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21


All times are GMT -4. The time now is 08:10.