|
[Sponsors] |
July 2, 2014, 12:00 |
Dealing with Source terms in turbulent flow
|
#1 |
Senior Member
Mehdi Babamehdi
Join Date: Jan 2011
Posts: 158
Rep Power: 15 |
Hi guys;
I applied a source term on the Non-zero pressure gradient flow over flat plate. the flow experiences natural transition to turbulence. The code worked for the flow without source term reasonably. By applying source term, the code diverges. the source term is synthetic jet on flat plate and not function of independent variables, hence they don't change during iterations. it has strong component toward -Y. by elimination of the Y component, the code works well, but after applying Y component, it diverges. I used slightly small under-relaxation numbers and different convetion discritization schemes (MUSCL, QUICK, TVD (mid-mode, superbee ,..)), but they did not make a big change (some of them postpone divergence). I want to know whether or not there is method can be used to handle this kind of source term. |
|
July 2, 2014, 12:39 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,832
Rep Power: 73 |
Quote:
The presence of a source term affects the stability constraint, have you checked for smaller parameters? Are you using DNS/LES/URANS? Furthermore, I suspect that it must be enforced in such a way to fulfill the continuity equation (divergence-free velocity) |
||
July 2, 2014, 22:29 |
|
#3 | |
Senior Member
Mehdi Babamehdi
Join Date: Jan 2011
Posts: 158
Rep Power: 15 |
Quote:
I used SST K-Omega (RANS) turbulence model. the developed code was evaluated for the case without source terms. I developed the code to simulate NZPG flow over a flat plate. And about your third statements (fulfilling the continuity equation), I can not understand what exactly you mean and what I should do to reach it. could you please explain more. |
||
July 3, 2014, 03:10 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,832
Rep Power: 73 |
I mean, is your source term in the momentum equation divergence-free?
Otherwise it must be included as source in the pressure equation |
|
July 4, 2014, 02:46 |
|
#5 | |
Senior Member
Mehdi Babamehdi
Join Date: Jan 2011
Posts: 158
Rep Power: 15 |
Quote:
the source term (included in "b" in momentum equations) impose on the momentum equation and indirectly (in U_hat and V_hat) affects the pressure equation. But as far as I know there is not any term in pressure equation on which source terms can be applied. Am I wrong? If I am wrong, could you please refer me to a source in which it is explained. |
||
July 4, 2014, 04:19 |
|
#6 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,832
Rep Power: 73 |
I cannot see without reading the code ...
but the general idea is: continuity: Div V = 0 momentum: dV/dt + conv + grad p = diff + b Therefore, taking the divergence of the momentum, the pressure equation is Div Grad p = source + Div b with proper bc.s including the term b.n on the boundaries |
|
July 5, 2014, 02:05 |
|
#7 | |
Senior Member
Mehdi Babamehdi
Join Date: Jan 2011
Posts: 158
Rep Power: 15 |
Quote:
I hope it will solve my problem. |
||
October 25, 2014, 22:23 |
|
#8 |
Senior Member
Mehdi Babamehdi
Join Date: Jan 2011
Posts: 158
Rep Power: 15 |
finally, I applied source term on the code I developed based on changing in changing in face velocities resulted from MIM.
It needs to decrease underrelaxation factor (in my case about 0.1), but the results is not promising, particularly pressure. pressure on the surface oscillates, and it produces wrong velocity. the problem is more severe upstream of the point where source term applies, and it is confusing because source term should not have any effects in those places. does anyone knows why it produces such oscillatory pressure? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] swak4foam building problem | GGerber | OpenFOAM Community Contributions | 54 | April 24, 2015 16:02 |
OpenFOAM without MPI | kokizzu | OpenFOAM Installation | 4 | May 26, 2014 09:17 |
centOS 5.6 : paraFoam not working | yossi | OpenFOAM Installation | 2 | October 9, 2013 01:41 |
pisoFoam compiling error with OF 1.7.1 on MAC OSX | Greg Givogue | OpenFOAM Programming & Development | 3 | March 4, 2011 17:18 |
Adding source terms to turbulent models | makaveli_lcf | OpenFOAM Running, Solving & CFD | 0 | June 8, 2009 09:34 |