|
[Sponsors] |
July 17, 2013, 12:54 |
K and Epsilon Do Not Change
|
#1 |
New Member
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 12 |
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 |
|
July 18, 2013, 03:21 |
|
#2 |
New Member
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 12 |
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 |
|
July 18, 2013, 17:47 |
|
#3 |
Senior Member
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 299
Rep Power: 22 |
Hi Matt,
Which solver are you using? Could you post the log-output of a few time steps? Cheers, L |
|
July 22, 2013, 10:24 |
|
#4 |
New Member
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 12 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 16 | March 4, 2017 08:30 |
calculation stops after few time steps | sivakumar | OpenFOAM Running, Solving & CFD | 7 | March 17, 2013 06:37 |
Calculation of k and epsilon freezes | Nigirim | OpenFOAM Running, Solving & CFD | 1 | November 14, 2012 07: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 20:21 |