
[Sponsors] 
January 30, 2012, 03:50 

#21 
Super Moderator
Oliver Gloth
Join Date: Mar 2009
Location: Todtnau, Germany
Posts: 116
Rep Power: 9 
The fact that the heat transfer is incorrectly predicted does speak in favour of a problem with the wall treatment.
To be honest, I have zero experience with the so called all y+ wall functions, but I have heard other people complain about it. You could try to test your wall function on a very simple flow (e.g. circular pipe); then try different values for y+ (1,5,10,30,100,...) and see how it influences the skin friction. If the all y+ wall functions are the problem you can then go for either low Re (y+ app. 1) or high Re (30 < y+ < 300) approach. From what I have read about your flow configuration, high Re is probably sufficient; I would try to aim for a y+ of maybe 50. 

January 30, 2012, 10:42 

#22 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
Of course, the form loss is a big contributor to the head loss. I am aware of this, but you do raise a good point. I checked the bare section pressure drop to the section with obstruction. A section having the obstruction has a pressure drop that is three times larger than a section without the obstruction (same lengths). This prompted me to investigate what is the culprit here. So I compared my calculated bare section pressure drop to experimental bare section pressure drop. It is still pretty far off (18% underpredicted). My calculated pressure drop for the bare section + obstruction was underpredicted by by ~27%. Just comparing the calculated and experimental pressure drop of the obstruction (by subtracting the pressure drop of the bare section from the section with the obstruction), I underpredict the pressure drop of the obstruction by ~8%.
This makes me think that the wall friction definitely has to be being predicted wrong. I tried running a case with komega, but it failed because I didn't properly select a free stream edge. That seems like opening up a whole other can of worms right now, so I'm going to put that on hold. I have a case running right now where I added a third prism layer and made the layers thicker to increase y+. I'm going to try what the other poster said and use a high y+ wall treatment and try and greatly increase my y+ values to about 30 or more. 

January 30, 2012, 12:25 

#23 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 469
Rep Power: 11 
What's the L/D for these sections?
Also, I am uncertain on your calculation of pressure drop percentage. Defininitions Experimental values: dp_eo: delta pressure (head loss) of the section with the obstruction of your experiment. dp_eb: delta pressure of a bare section of your experiment (a section without an obstruction) Calculated values dp_co: delta pressure (head loss) of the section with the obstruction of your calculation. dp_cb: delta pressure of a bare section of your calculation (a section without an obstruction) OK, so dp_eo = 3*dp_eb 18% = 100*(dp_cbdp_eb)/dp_eb 27% = 100*(dp_codp_eo)/dp_eo Therefore the error of your wall friction of your bare section relative to the obstructed region is: 100*(dp_cbdp_eb)/dp_eo = (1/3)*(100*(dp_cbdp_eb)/dp_eb) = 6% Therefore the wall friction accounts for only 6% of the experiment and the obstruction alone accounts for 21%. Or is my original understanding wrong? But, I find the following hard to believe 18% = 100*(dp_cbdp_eb)/dp_eo Since that would mean, 100*(dp_cbdp_eb)/dp_eb = 54%!! 

January 30, 2012, 13:53 

#24  
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
This was right...
Quote:
This... Quote:
The pressure drop of the bare section was measured. The pressure drop of the obstruction was not "directly" measured. It is calculated by subtracting the bare section pressure loss from the pressure loss of a bare+obstruction section... dp_eoo = dp_eo  dp_eb where dp_eoo is the pressure drop across only the obstruction and dp_eo is the pressure drop of both the bare section and the obstruction that is within the section. The L/D of a section (including the bare portion and the obstruction) is about 23. The L/D of the obstruction is less than 1. 

January 30, 2012, 14:23 

#25  
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 469
Rep Power: 11 
Quote:
Proof: I assume this is how you calculated your overall loss 27%=100*(dp_codp_eo)/dp_eo and I'm assuming dp_eo = 3*dp_eb using your definition dp_eoo = dp_eo  dp_eb and therefore dp_eo=dp_eoo+dp_eb and creating another definition dp_coo = dp_co  dp_cb and therefore dp_co = dp_coo+dp_cb therefore 27%= 100*(dp_codp_eo)/dp_eo = 100*((dp_coo+dp_cb)(dp_eoo+dp_eb))/dp_eo = 100*((dp_coodp_eoo)+(dp_cbdp_eb))/dp_eo= 100*(dp_coodp_eoo)/dp_eo + 100*(dp_cbdp_eb)/dp_eo= 100*(dp_coodp_eoo)/dp_eo + 100*(dp_cbdp_eb)/(3*dp_eb)= 100*(dp_coodp_eoo)/dp_eo + (1/3)*100*(dp_cbdp_eb)/dp_eb restating it 100*(dp_coodp_eoo)/dp_eo + (1/3)*100*(dp_cbdp_eb)/dp_eb = 27% therefore 100*(dp_coodp_eoo)/dp_eo = 27%  (1/3)*100*(dp_cbdp_eb)/dp_eb therefore 100*(dp_coodp_eoo)/dp_eo = 27%  (1/3)*(18%) = 21% So, Your loss due the obstruction for the obstruction section is 21% and not 8%. Your loss on the channel walls for the obstruction section is 6%. 

January 30, 2012, 15:46 

#26 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
I appreciate your effort on coming up with that proof, but please hear out my methodology. First, I must appologize because when I wrote that previous post saying that the total section pressure drop was 3 times that of the bare section, I was looking at the wrong pressure transducer. The total section pressure drop is actually only ~1.5 times the bare section.
However, what I did was I looked at the calculated pressure drop over a bare section (dp_cb) and I looked at the calculated pressure drop over a section having the bare portion and the obstruction (dp_cbo). I calculated the pressure drop of the obstruction by taking dp_cbo  dp_cb. I did the same for the experimental results (having a bare section measurement, dp_eb, and a bare + obstruction section, dp_ebo). The pressure drop of the obstruction is dp_ebo  dp_eb. Comparing the obstruction pressure drop for experimental and calculated, in this manner, gave me a difference of about 6%. Comparing the bare pressure drop for experimental and calculated gave me a difference of about 18%. I think, at the least, it's very clear that there is something wrong with the modeling of the bare region, as the difference between experiment and CFD that I'm seeing is an apples to apples comparison and the difference is quite large. I have an additional case submitted using a thicker prism mesh layer to increase the y+ values (and using the high y+ wall treatment). It is still waiting to start. I think I will have results tomorrow. I will post back with the effects of playing with the prism layer and the wall treatment. 

January 30, 2012, 17:11 

#27 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 469
Rep Power: 11 
Yes, please do post back! This is interesting.
LOL, I'm a man of pictures, equations, and plots. This painting with words is just not my cup of tea! Just for kicks I've attached the picture in my head. Sorry for horrible penmanship! Too many computers and spell checkers in my life. 

January 30, 2012, 17:33 

#28 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
The drawing you posted is a little off. The 'obstruction zone' that I was referring to is the area in the direct vicinity of the obstruction (red rectangle in the picture below). Everything else would be the bare flow area (blue area). The pressure measurement is made over bare section + obstruction. Note that there is also a pressure measurement made simply over a bare section (not shown), which is what I was using to determine how accurate the prediction of bare region pressure loss is (which should have only frictional loss, since there are no flow blockages present in that section).
As I believe I mentioned before, there was a case that was unheated and a case that was heated. Pressure prediction was bad in both cases. The walls aren't perfectly adiabatic (close enough though); however, for the unheated case that wouldn't make a difference anyway. 

January 30, 2012, 18:29 

#29  
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 469
Rep Power: 11 
Thanks for the clarification.
Quote:
For the picture I showed, and yours, there is a good chance for a recirculation region to occur both in front and behind the obstruction. This will create low (flow stagnating) or negative (circulation region) skin friction on the channel wall. You'll know by looking at the velocity contour plots. I didn't think wall functions where good at this. But I don't know the ins and outs either. I believe, to capture this, you will need to go with, at a minimum, SA or SST. That's why I'm say y+ < 1. That's also why I'm saying you need to be careful about capturing the obstruction and the channel region around it. I'm not sure what a "free stream edge" is, but that's probably because I'm an external aerodynamicst. SST may give you better results than SA, but SA is easy to use, at least from my perspective. Though, SST isn't that much harder. More of an issue of CPU time. But, SA does only have one variable, eddy viscosity (or it's form of it) so passing eddy viscosity from one turbulence model to SA might be easy. Not sure about SA to other turbulence models. And, I'm not sure how you are building up the velocity and eddy viscosity profiles for the input to your channel. Finally, maybe "wall functions" means something different to me than you. Discussion about kepsilon and komega: Use of kepsilon and komega Models Good Luck! 

January 31, 2012, 11:01 

#30 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
Thanks for the link. The first case that I ran, where I added a third prism layer, didn't solve nicely at all. After 2E6 CPU seconds, the solution only reached 230 iterations and crashed due to insufficient memory. It seems odd to me that simply adding another prism layer would cause such strange behavior (the total cell count is 6.7 million vs. 6.2 million the first time I ran this case (which was successful)). I can't remember changing anything else besides the mesh. Anyways, I'm still waiting to get results from the unheated case using the thicker prism layers (higher y+ values and high y+ wall treatment).


January 31, 2012, 17:29 

#31  
Senior Member
Vieri Abolaffio
Join Date: Jul 2010
Location: Always on the move.
Posts: 308
Rep Power: 9 
Quote:
i'd say that 3 cells are not enough even when WF are used, but with a loRe turbulence model you'll need at last 20 cells in the BL could you please post, for the sake of clarity, a screenshot of your BL mesh near the wall of the pipe? 

January 31, 2012, 17:50 

#32 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
As far as I know, no I'm not fully resolving the boundary layer. I selected the all y+ wall treatment. It's a wall function. And the thickness of those two prism layers (originally 2) was about 2.5E4 m, I believe, with the nearwall layer being about 1E4 m. I have another case waiting where I increased that thickness.


March 19, 2012, 18:03 

#33 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
Sorry for taking so long to reply back to this thread, but I got tied up in some other work and never ended up getting to that thicker wall prism layer case until now. To sum up the overall results I'm looking at so far:
1. I ran 3 cases, a base case with 2 prism layers, total thickness of layers was 2.5e4 m, a refined case with 3 prism layers having total thickness of 1.9e4 m, and a thicker case having 2 prism layers w/ total thickness of 5 e4 m. 2. No significant difference in prediction of pressure loss was seen between any of the three cases 3. Crosssectional averaged velocity in different locations of the test section were consistently significantly better using the refined mesh So it doesn't appear my meshing is affecting the pressure drop. Maybe it comes down to the turbulence model. 

March 23, 2012, 18:49 
Some thoughts on you problem

#34 
New Member
Aaron
Join Date: May 2011
Posts: 6
Rep Power: 7 
So I am currently working on my Master's Degree doind simulations on a similar setup thought it sounds like for a totally different application. Also it sounds like you are using StarCCM+ as your solver.
To clarify a little on the wall funtions the all y+ wall treatment is supposed to determine based on y+ whether a wall funtion is necessary or not. That is you should be able to have bl cells with a y+ that is valid for both no wall funtions and wall funtions and the solver will determine the appropriate handling. My personal experience however is that it is not very successful at doing this. It works better if you plan out your simulation to be either a fully resolved bl (y+ ~1 everywhere) or wall funcitons (y+ ~30100 everywhere). In practice either of these things can be very difficult to accomplish. That is you will likely have some bl cells that are not appropriate for wall functions of a fully resolved boundary layer solution. In these cases it has been my experience the the all y+ model works well. My recomendation would be to generate a mesh that fully resolves your boundary layer everywhere on the channel and the obstruction (average y+ ~1) and select the all y+ model. When you do this you will need to pay close attention to cell aspect ratio, try to keep it >= 0.1 as much as possible. Keep in mind thought that to much larger than that and you are probably wasting cells. From my expericence the trimmer meshing model provides the easiest controls to generate these meshes and trimmer meshes tend to solve fast as a bonus. Also if you make a mesh like this you will need at least 12 bl cells. I would agree with you that it is your wall treatment that is the source of your problems, and that it is likely your mesh that is the issue as this is the most determining factor on how well a boundary is simulated. I was having similar problems with my work and these practices have greatly improved the quality of my results. 

March 23, 2012, 18:56 

#35 
New Member
Aaron
Join Date: May 2011
Posts: 6
Rep Power: 7 
Also I forgot to ask you what level of convergence are you getting in your residuals?


March 23, 2012, 20:03 

#36 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
I'm going to have to get a lot better at making my meshes. Having 12 prism layers to resolve the BL is just not going to work for my model. It would be way too computationally expensive. Converging the models I have now take about 7 million cpuseconds. I tried to jack up the y+ values by increasing the prism layer thicknesses, but even after upping the thickness to a ridiculous amount, I still am left with a ton of y+ values below 1. I rationalize that it's due to the very low velocities in some regions of my models where there are crevices and due to the definition of the y+ parameter. So I don't know how else I can increase y+.
I finished compiling my results recently and I ended up with 3 different ways to model the same heated test: 1. a base mesh of 7 million cells 2. a refined mesh of 15 million cells (added a prism layer to make 3 layers total and reduced the thickness of them to get more refinement near the wall) 3. the same mesh as #1 but using standard kepsilon instead of realizable kepsilon Results show that using configuration 2 didn't change the prediction of the pressure drop. Nor did it improve prediction of the wall temperature. But the temperature gradient throughout the channel crosssection became more uniform than configuration 1. Using standard ke, though, had a significant impact on pressure drop prediction, improving it by 10%. Also, it greatly improved prediction of wall temperatures and led to a much more uniform temperature gradient. The residuals for my cases look pretty decent. The better indicators, I think, are the inlet/outlet mass balance, channel velocities w/ respect to iterations, and channel temperatures w/ respect to iteration. They all go to constant values, giving the indication of convergence. 

March 24, 2012, 17:02 

#37 
New Member
Aaron
Join Date: May 2011
Posts: 6
Rep Power: 7 
I would have to take a slightly different position than you on what your three cases mean. Your first two sound to me like they are handling the boundary region almost identically. Two or three prism layers are almost the same. And yes you did refine them but not in the manner that would make there handling of the near wall region significantly different. Thus it makes sense to me that the results were similar. As far as changing the turbulence model from the realizable to the standard kepsilon model I believe this actually supports the idea that it is your bl mesh. My understanding (and I didn't take the time to verify this) is that the difference between the two models is entirely in how the near wall region is handled. This tells me simply that for your meshing strategies the standard kepsilon modle is better suited for your boundary layers. The question then becomes are the results you are getting sufficiently accurate for your application. Only you can answer this but if they are not I am almost positive that it is only by significantly improving your mesh that you can get truly accurate results.
A few ideas that I have learned along the way are 1) when in doubt always assume the mesh is your problem, this is the single most influencable aspect of CFD and it has critical implications for accuracy. 2) Boudnaries are everything. If you think about it almost all CFD is 95% identical, you start out with the same set of equations. The 5% that is differnet are the boundary conditions (sometimes also initial conditions) these conditions contain all of the information that distinguishes flow over and airfoil from blood flow in an artery (minus constituitve models). Boundary condition information is the solution. If you are only moderately modelling your boundary conditions nothing you do will give you better than moderate results. And for the most part the success of modeling the boundary for a no slip wall depends enritely on the mesh. I just don't think you can ignore the need to resolve your bl, I am speaking from experience as my own research closely mimics your own. I am studying vortex generation as a means of increasing mixing in straight channels. Until i bucked down about my bls i wasn't getting very good results. Now that being said I was suprised that by adding one prism layer and refining the cell size you went from 7 to 15 million cells. I don't think that is necessary to get a really refined bl. True you will increase mesh size but you should be able to do it a lot more reasonably. If you are using Star then set your refence values for your prism layer boundary mesh seperately from those of the region. The prism mesh should be made based on absolute numbers and not relative to the base size. Also at a very minimum your obstruction and your channel walls need bl meshes designed specifically for them, you cannot make your prism layers for both with the same values. If you change the target surface size on you walls and the growth rate from the walls you can likely have a well refined bl without making your channel cells too small as well. Anyway just some thoughts. If your results are already good enough then all this is acedemicm, but if not I am almost positive that your mesh is what will make all the difference. One final thought is that for confidence your residuals really all need to be below 1e3. Sorry about the length, just some thougts. 

March 24, 2012, 18:55 

#38  
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
Thank you for the information. I agree that the mesh is a likely culprit considering the problems that I had with just obtaining convergence that were solved by mesh modification. At this point, I'm just not really sure how to model my near wall region any differently to expect any different wall drag and heat transfer results.
Quote:
Quote:
Quote:
Quote:


March 24, 2012, 20:33 

#39 
New Member
Aaron
Join Date: May 2011
Posts: 6
Rep Power: 7 
Two things, it really should be possible to generate a mesh the completely resolves the bl and is only 3035% larger than your current meshes. Which while larger shouldn't be completely prohibitive. That was going from no prism layer to 12 prism cells. When I switched from trying to have y+ of about 30 to y+ of 1 it was an increase on this order. I really think you need to play around with their meshing tools until you get a good balance having been there I know it is possible. Also I sould sacrifice some of the resolution in the domain in order to get better resolution at the wall if I had to.
In regards to the residuals I have heard those arguments before, in fact I will be joining CDAdapco as an application engineering in about a month. I know that for them they are big on the engineering parameter as a means of judging convergence. In my opinion you have to use both engineering parameters and residuals. You may already be aware of this but residuals are not simply the normalized change in solution from iteration to iteration. They are in fact a normalized represention of how completely the solution satisfies the governing equations. This is of critical importance because the governing equations are the only reason we have for attempting to model these phenomena at all. The bulk of Star users are on the industrial side of things and tend to be more satsified with CFD that is moderately accurate, after all this often times is all they require. Thus they focus more on the flow parameter of interest and don't bother with the rigeur required to converge the residuals. However, on the acedemic side this is not the case, and if it is your intention to publish your research I think you will find that journal reviewers care very much about the residuals. The number I stated is based of the minimum requirements for publication in Journal of Fluids Engineering and is only one of several requirements (they also have guidlines regarding the flow parameter of interest). Until your residuals are truly are converged satifactorily I don't think you should expect to get too much more accurate results, in industry + 10% is often acceptable, but it sounds like for you it is not. If I am just repeating information to you that you already know I apologize. Just trying to be helpful. 

March 25, 2012, 15:11 

#40 
Member
Join Date: Dec 2011
Location: State College, PA
Posts: 87
Rep Power: 6 
I was not aware that it would be possible to fully resolve my BL with such few cells. It's something I'm just going to play around with in the future. I won't have the time to do that for a while, but it's something I will consider when I get back to this. One thing that confuses me is, what percentage of my prism layer cells should have a y+ value of my target value? I don't see how it is possible to have every cell have a y+ value of 30, say. Because the prism layers need to shrink in many areas where surfaces get close together. Also, the velocity is just very low in some regions that causes y+ to go to values of like 0.01, even when I'm targeting y+=30.
I do keep the residuals in my mind when judging convergence. In fact, it's the very first thing I look for. And if they're oscillating a lot, even if my engineering parameters look reasonable, I'll try out something different with the mesh or my model to fix that behavior. It's just that I wasn't requiring all residuals to drop below a specified value so long as my engineering parameters look good. Thanks for all of your advice. I will keep these things in mind when I'm meshing up more cases. 

Tags 
pressure drop, wall drag 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Pipe Flow  Pressure Drop  Daniel L  FLOW3D  2  December 10, 2010 05:23 
Total and Static Pressure drop  Jaggu  FLUENT  3  December 9, 2008 07:03 
Estimating pressure drop from the shear stresses  K.Baker  FLUENT  2  July 26, 2007 03:56 
Periodic boundaries and pressure drop  MatejK  CFX  2  July 12, 2005 04:47 
Hydrostatic pressure in 2phase flow modeling (CFX4.2)  HB &DS  CFX  0  January 9, 2000 14:19 