[ogloth]

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.

[rks171]

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% under-predicted). My calculated pressure drop for the bare section + obstruction was under-predicted 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 under-predict 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 k-omega, 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.

[Martin Hegedus]

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_cb-dp_eb)/dp_eb
-27% = 100*(dp_co-dp_eo)/dp_eo

Therefore the error of your wall friction of your bare section relative to the obstructed region is:

100*(dp_cb-dp_eb)/dp_eo = (1/3)*(100*(dp_cb-dp_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_cb-dp_eb)/dp_eo

Since that would mean,

100*(dp_cb-dp_eb)/dp_eb = -54%!!

[rks171]

This was right...

-18% = 100*(dp_cb-dp_eb)/dp_eb

That's how I calculated the error of the calculated bare section pressure drop.

This...

-18% = 100*(dp_cb-dp_eb)/dp_eo

...I did not do. I'm not sure what the reasoning would be for dividing by the pressure drop of the bare+obstruction section.

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.

[Martin Hegedus]

Originally Posted by rks171

...I did not do. I'm not sure what the reasoning would be for dividing by the pressure drop of the bare+obstruction section.

I agree with that! So the error of your obstruction is 21% rather than your claimed 8%.

Proof:
I assume this is how you calculated your overall loss

-27%=100*(dp_co-dp_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_co-dp_eo)/dp_eo =
100*((dp_coo+dp_cb)-(dp_eoo+dp_eb))/dp_eo =
100*((dp_coo-dp_eoo)+(dp_cb-dp_eb))/dp_eo=
100*(dp_coo-dp_eoo)/dp_eo + 100*(dp_cb-dp_eb)/dp_eo=
100*(dp_coo-dp_eoo)/dp_eo + 100*(dp_cb-dp_eb)/(3*dp_eb)=
100*(dp_coo-dp_eoo)/dp_eo + (1/3)*100*(dp_cb-dp_eb)/dp_eb

restating it
100*(dp_coo-dp_eoo)/dp_eo + (1/3)*100*(dp_cb-dp_eb)/dp_eb = -27%
therefore
100*(dp_coo-dp_eoo)/dp_eo = -27% - (1/3)*100*(dp_cb-dp_eb)/dp_eb
therefore
100*(dp_coo-dp_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%.

[rks171]

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.

[Martin Hegedus]

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.

[rks171]

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.

[Martin Hegedus]

Thanks for the clarification.

Originally Posted by rks171

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.

Yup, that was understood. That's why I focused on what goes into the pressure. Plus, the energy equation is not used for your incompressible unheated case.

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 k-epsilon and k-omega:

http://www.cfd-online.com/Forums/mai...ga-models.html

Good Luck!

[rks171]

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).

[sail]

Originally Posted by rks171

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).

i might have lost something in the thread, but you are running a fully resolved boundary layer simualtion (y+ < 1) with just three cells in the BL?

i'd say that 3 cells are not enough even when WF are used, but with a lo-Re 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?

[rks171]

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.5E-4 m, I believe, with the near-wall layer being about 1E-4 m. I have another case waiting where I increased that thickness.

[rks171]

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.5e-4 m, a refined case with 3 prism layers having total thickness of 1.9e-4 m, and a thicker case having 2 prism layers w/ total thickness of 5 e-4 m.
2. No significant difference in prediction of pressure loss was seen between any of the three cases
3. Cross-sectional 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.

[USU_CFDer]

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 Star-CCM+ 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+ ~30-100 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.

[USU_CFDer]

Also I forgot to ask you what level of convergence are you getting in your residuals?

[rks171]

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 cpu-seconds. 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 k-epsilon instead of realizable k-epsilon

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 cross-section became more uniform than configuration 1.

Using standard k-e, 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.

[USU_CFDer]

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 k-epsilon 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 k-epsilon 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 1e-3.

Sorry about the length, just some thougts.

[rks171]

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.

I just don't think you can ignore the need to resolve your bl

With people saying that I'm going to need something like 12 prism layers to resolve the boundary layer, I just don't think that's going to be feasible for my model to resolve the bl. I simply do not have access to the computational requirements for that kind of mesh. Therefore, I need to rely on wall functions and optimizing my mesh to work with them.

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 didn't only make the prism layer base size smaller, but I also reduced cell sizes throughout the entire mesh. This led to better predictions of velocity. Also, more energy was pulled away from the near wall region and made for a more uniform channel temperature gradient.

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.

Yes, I do all of this already. Particularly, I use a more refined mesh in the obstruction region. But thank you for pointing out these points.

One final thought is that for confidence your residuals really all need to be below 1e-3.

This last point, I'm really not so sure about. I've read in several places on here and in the Star-ccm+ manual that the actual values of the residuals aren't necessarily important and that you can still obtain convergence without having residuals reduced to some specific value. This is due to the fact that the residuals are normalized values. This is why I also watched actual flow field parameters and made sure they were settling into consistent values (if they should be behaving like that).

[USU_CFDer]

Two things, it really should be possible to generate a mesh the completely resolves the bl and is only 30-35% 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 CD-Adapco 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.

[rks171]

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.