# u plus vs y plus

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

 November 19, 2009, 05:21 u plus vs y plus #1 New Member   Jin Chen Join Date: Mar 2009 Location: sweden Posts: 27 Rep Power: 9 Dear all, I have a question about plotting u plus and y plus. I creat a line on floor with a certain height, define u plus and try to plot u plus vs y plus, but i can only get one pair data.(Based on y plus defination is the distance to the nearest wall, so i can only get one pair data) But if i want to see how does uplus vary with y plus, like the linear and the log relatoin, what should i do? Best regards Chne

 November 19, 2009, 09:45 #2 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 544 Blog Entries: 14 Rep Power: 18 What you want to plot is actually nothing different from the classical velocity profile in the direction normal to the wall, but in different coordinates. Hence, if you have the pairs u - y (for different values of y!), you have to plot u/u_tau vs. y * u_tau/(mu/rho), where y is the wall normal coordinate. You are probably making confusion between the y+, the transformed wall normal coordinate, and the dy+, the transformed wall normal coordinate of just the first point off the wall.

 November 20, 2009, 05:34 #3 New Member   Jin Chen Join Date: Mar 2009 Location: sweden Posts: 27 Rep Power: 9 i am not so clear about you mentioned dy+(do we need this variable when plotting u+ vs u+). I mean i ploted u+ vs y+ along a line which is perpendicular to the floor,but i can only get one pair value, seems like only for the neasrest node close to the floor. y plus is meaningful for the nearest node or for all the nodes ? If we look at its definetion , it's calcuated from the wall shear tress(U_tao=roor square(wall shear stree/density)), but wall shear stress happens on the wall. so how i can get the y plus for the rest of the nodes except for the neasrest one? or i misunderstand something?

 November 20, 2009, 07:41 #4 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 544 Blog Entries: 14 Rep Power: 18 Probably now i get the point: if you are using Fluent to do this, it will work as you described as it has y+ defined only at the near wall cells hence it cannot give you the values all along the specified lines. However this is just a Fluent feature (mostly due to the memory saving necessity) and of course y+ is meaningful all along the wall normal direction as, for example, the wall laws are defined in terms of y+. Actually, in fluent, you are just plotting y+ vs y which is dy+ vs y as it is given only in the near wall cell. The best thing you can do is to export your velocity profile (u vs y) in a file and then use, say, matlab to make the coordinate transformation with the wall stress computed by fluent. The wall stress is just the value at the wall and it is obviously the same all along the line Otherwise, i don't get why you should get just a point instead of the whole line you created Far and manukamin like this.

 November 20, 2009, 08:20 #5 New Member   Jin Chen Join Date: Mar 2009 Location: sweden Posts: 27 Rep Power: 9 Thanks very much for your advice! Now it's clear!

 December 10, 2014, 13:50 Problem solved! #6 Senior Member   Shamoon Jamshed Join Date: Apr 2009 Location: Karachi Posts: 160 Rep Power: 9 Sbaffini! Excellent reply of plotting the curve in Matlab or Excel. You made my day. The thing I had been looking for years is solved today. Hats of for you sbaffini likes this.

 February 22, 2016, 05:05 #7 Senior Member   Shamoon Jamshed Join Date: Apr 2009 Location: Karachi Posts: 160 Rep Power: 9 Sbaffini, how are you? I hope you will help me again. Here is the problem I am stuck at: " I have modeled a quarter pipe with its solid thickness modeled too. I have applied the flux over wall as 3500 W/m2. I have the formula for havg= flux/(Tavg,fluid-Tinside). The T inside is the temperature at the lower end of wall, the region that touches fluid. In experiment Tinside is measured via thermocouples that are installed across the tube length wherever the flux is applied. These TC however measures the inner temperature of the wall. Nusselt number is Nu=havg*Di/k. Di=inned dia, and k=thermal conductivity of fluid. My computed Nu is very low from experiment. I am using Fluent, and applied flux on top wall. and specify coupled BC on interface. This creates wall shadow.I do not use any heat value (flux or temp) except on top wall. I saw as well "via System coupling" written in these BC. do I need to turn it on too? In fluent I take Area weighted average for inlet and outlet (Ttotal) while for Thermocouples I made Isolines at TC positions and took Temp (total) with Area weighted average. The error is more than 10 % from experiment." Please help me as I am not able to find a solution yet.

February 22, 2016, 05:31
#8
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,303
Rep Power: 28
Quote:
 Originally Posted by swe704 Dear all, I have a question about plotting u plus and y plus. I creat a line on floor with a certain height, define u plus and try to plot u plus vs y plus, but i can only get one pair data.(Based on y plus defination is the distance to the nearest wall, so i can only get one pair data) But if i want to see how does uplus vary with y plus, like the linear and the log relatoin, what should i do? Best regards Chne

y* is the non-dimensional distance multiplied by Re_tau, this way you can plot u*(y*) in a log plot

 February 23, 2016, 09:36 #9 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 544 Blog Entries: 14 Rep Power: 18 Dear Shamoon, there might be several causes for your results. Unfortunately it's been a while since the last time i used Fluent. I suggest you to first check your setup (or a similar one) against a known laminar solution (for example, the CHT problem of a laminar poiseuille flow with an external solid tube can be solved analytically). An additional check might be to add an additional solid layer in your tube, made up of a different solid material, and adapting the b.c. to this new layer. If everything works properly, you should get the same solution in the original domain part. If everything is ok, the problems are either concerned to the postprocessing (which i experienced in the past for the CHT) or the turbulence modeling (which is unlikely, due to the very simple geometry). Concerning the postprocessing, the problems i had were related to the way the local heat flux was produced in output (both contours and files) for coupled walls. I finally decided to do this postprocessing by myself (actually a student i was supervising ), by having fluent write the temperature variables in celle centers on both sides of the surface and reconstructing the actual flux from the known interpolation scheme used by fluent. I don't know if the wrong behavior i experienced was determined by the fact that i was also using time statistics. Moreover, it was also fundamental to understand clearly what was the experimental output and how was it computed, in order to do the postprocessing in a consistent way. Hope this is helpful.

 February 23, 2016, 13:51 #10 Senior Member   Shamoon Jamshed Join Date: Apr 2009 Location: Karachi Posts: 160 Rep Power: 9 Dear Sbaffini I finally decided to take the interface temperature all along the pipe length and not the thermocouple locations. So error is observed at somewhat large Re . The thing you say, that do interpolation by yourself, may be I give it a try too. Strangely, the two interfaces, from which fluent creates wall and wall shadow have different temperatures. Please note that I modeled the wall thickness and meshed it as well and named the interface zones. I apply the flux on the outer side of the wall that has solid adjacent cell zone. Its very strange to me that the two interface are like two sides of a same coin and give different temperatures. 1. Could you tell me how Fluent does interpolation? Can I do it by myself? 2. Any particular case of pipe flow with wall thickness in your opinion? Best Regards,

February 24, 2016, 13:42
#12
Senior Member

Shamoon Jamshed
Join Date: Apr 2009
Location: Karachi
Posts: 160
Rep Power: 9
Dear sbaffini,

Thanks for your kind reply. I can try your method as well. The thing I mentioned yesterday about two different temperatures at interface was because I was taking the total temperature. However, static temp was same on both sides. Although the difference between static and total is not much.
Secondly, the formula in Nu is havg*D/k where havg=flux/(Tavg,w-Tavg,f)
Tavg,w = mean of the inner wall temperature
Tavg, f = mean of inlet and outlet

What I did was that first blindly followed the formula. I took area weighted average at 8 thrmocouple locations and noted the temperature at those stations (those were lines formed at interface).
I took mass weighted average btw inlet and outlet
Using this method gave me error range from 18-23%.
For second iteration, I took interface and took area weighted avg of total temperature. I took MWA value for second term just at pressure outlet. This improved my results. And now I am sticking to it because the experiment guys said themselves that the error in Nu is +- 18%. I am sending you some images from paper, that may be help you in understanding my prob.
Attached Images
 formula.png (37.0 KB, 6 views) formula2.png (67.6 KB, 6 views) pipe.png (13.4 KB, 7 views) uncertainities.png (8.1 KB, 5 views)

 February 25, 2016, 14:27 #13 Senior Member   Shamoon Jamshed Join Date: Apr 2009 Location: Karachi Posts: 160 Rep Power: 9 Ok, one more thing, I also have to compute friciton factor. FF depends upon pressure and I put the operating as 101325 and gauge press at outlet to zero. This gives me correct pressure at velocity inlet which is infact delta p in friciton factor formula. Problem is that even after 25000 iterations it is dropping down, but very slower rate, so at which point should I stop?

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

All times are GMT -4. The time now is 00:51.