
[Sponsors] 
March 19, 2019, 20:34 
Effect of coarse mesh on drag through valve

#1 
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Hi, something I have been wondering about is what would be the expected effect of using a mesh that is too coarse, particularly with regards to the 'drag' measured through a valve? Realizable kepsilon is the turbulence model that I am particularly interested in.
My thinking is that an overly coarse mesh would probably add in some artificial 'numerical viscosity', because it will model the boundary layer as being thicker than it really is. Also, it will probably give a bad prediction of boundary layer attachment as the flow turns round corners. Therefore, the drag would tend to be too high, and so it would underpredict the flow rate through a valve for a given pressure drop. Does that sound correct, or is there a flaw in my reasoning? 

March 20, 2019, 03:22 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,151
Rep Power: 66 
Quote:
The numerical viscosity is an issue appearing from the discretization of the convective terms rather then from the grid (even if it depends on the grid size). Actually, you can have also discretizations that do not produce numerical viscosity. As you are working with RANS formulation, the additional viscosity is produced by the turbulence model. The effect of the grid resolution on the viscous drag can be easily seen by considering the tangential stress that requires the discretization of the normal derivative of the streamwise velocity. If the noslip condition is prescribed the derivatives depends on the distance of the grid size from the wall. If you do not describe correctly the viscous sublayer such derivative is poorly computed. 

March 20, 2019, 07:57 

#3  
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Quote:
Ok. So, let's say I use a grid that is very coarse and doesn't adequately resolve the boundary layer (no inflation layers, very large wall y+ ...), does that mean the computed velocity gradient at the wall would be too low? I.e. it would actually underestimate the frictional drag? Would it also depend on what type of wall function I use? Most typically, I use standard or nonequilibrium wall functions in FLUENT for the valves I analyze. 

March 20, 2019, 08:05 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,151
Rep Power: 66 
Quote:
 yes, the stress should result lower than the actual  If you apply any sort of wall modelled BCs (instead of the natural noslip velocity BC) there is no further discussion about the stress at the wall: you are somehow already prescribing it and you do not describe the viscous sublayer. Thus, you cannot pretend to predict what you have already prescribed as BC. 

March 20, 2019, 08:28 

#5  
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Quote:
I'm not sure I fully understand what you are saying here:  Surely, even with a wall function, the nonslip condition at the wall still has to be respected?  Also, doesn't the nearwall velocity profile that is modeled by the function still have to be determined by the solver, based on the wall y+ and 'freestream' velocity? So, how is the wall shear stress being 'prescribed' in that situation, where a wall function is being used? 

March 20, 2019, 09:21 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,151
Rep Power: 66 
Quote:
Generally speaking, the wall modelled BCs do not prescribe the noslip velocity at the wall (y+=0)! Conversely, it is assumed you have a fully developed averaged velocity profile and considers the wall line like a fluid line at some y+ value. In other words, you cannot think to predict the drag when you are assuming that you already know a specific type of velocity profile. Of course, the flow behind a valve is very far to be similar to the assumed wall laws. I suggest to serach for similar posts in this forum 

March 20, 2019, 10:38 

#7  
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Quote:
Hmm. What you say is interesting and it differs from my understanding of how the kepsilon wall functions work. I will have to do some more background reading into the subject. My understanding is that the wall function is assuming a known normalized velocity profile in the boundary layer, but that the profile still has to be scaled during the solution to give an accurate representation of the boundary layer thickness. So, it seems to me that the predicted shear stress might still depend on the mesh resolution. However, perhaps I have some misconception. Out of interest, what type of RANS turbulence model/wall function would you recommend for modeling an industrial flow through a globetype control valve (which has a lot of swirling and separation)? 

March 20, 2019, 11:03 

#8  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,151
Rep Power: 66 
Quote:
When you assume that only a "scaling" must be determined on a know velocity profile what do you think are you prescribing? The key to understand is that a fully developed boundary layer profile is supposed and of course this is not true in a lot of wall turbulence flows. The flow around a valve is rich of separated flow region, a case in which DES/LES is suitable. 

March 21, 2019, 08:39 

#9  
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Quote:
As I understand it, the prescribed velocity profile has to 'match up' with the freestream velocity away from the wall, which is calculated during the solution. Otherwise, the velocity profile would be discontinuous and unrealistic. That is what I meant by 'scaling': the normalized velocity profile is fixed (u+ = y+ in the viscous sublayer), but u+ is normalized by the freestream velocity. Even in a fullydeveloped boundary layer, the wall shear stress will be higher, if the freestream velocity is higher. I see what you mean though, that a fullydeveloped boundary layer is assumed, which will not be accurate in some situations. So, you would not recommend using RANS at all for modeling flow through a valve? I agree that DES/LES would be ideal; however, I don't think they are really practical in an industrial setting. 

March 21, 2019, 09:59 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,151
Rep Power: 66 
Quote:
You can read about the wall functions and the limits in this approach in the book of Wilcox. Considering at present DES/LES for a flow around a valve should be considered almost mandatory at present, even in industrial applications. This study (https://www.researchgate.net/publica...utomobile_2003) was done 16 years ago on a notebook, almost as a joke, don't you think that now the computational power can do better ?? 

March 22, 2019, 18:31 

#11  
Senior Member
Lee Strobel
Join Date: Jun 2016
Posts: 133
Rep Power: 8 
Quote:
Thank you for linking that paper, it is very interesting. I am curious though, how many cells did the mesh contain (approximately)? Also, how finely was the nearwall region modeled (what was the wall y+)? (I didn't see those mentioned in the text). Comparing that example to the valves that I simulate on a daily/weekly basis, I have to say the valves I work with are quite a lot more complex than that geometrically. Typically, I use a halfsymmetric model with 510 inflation layers (maximum wall y+ about 200). The cell count is typically in the range of 56 million cells. I work for a large multinational and have access to a modern supercomputer. Using 40 cores, a steady, incompressible solution in FLUENT using RANS (kepsilon) takes about an hour to an hour and a half. The thing is, it's not just about the solving time. There is also time involved with setup, pre and postprocessing. I would expect LES to take longer to solve, as it is transient, especially if the boundary layer needs to be more finely modeled. But then, on top of that, there will be more processing time required to deal with transient results, averaging, etc. Plus, the large increase in the amount of data that will be generated and need to be stored/handled, compared to a simple, steady result. Of course, I accept that DES/LES would be more accurate than RANS in most situations. But, if it will increase the amount of time to setup and perform a solution by 34x, then that just isn't practical for me, given the time and resource constraints that I have to work under. I often need to do 5, 6 or more valve internal flow simulations in a week. That is why I was asking, if I have to use RANS due to my constraints (even though it is not as accurate), what turbulence model and wall modeling method would you recommend for an internal (incompressible) flow through a control valve? 

Tags 
drag, kepsilon, mesh, valve 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
IC engine valve dynamic mesh modelling  soogeng  FLUENT  0  April 14, 2018 10:27 
udf for valve closing a pipe using dynamic mesh  chem engineer  Fluent UDF and Scheme Programming  2  May 13, 2017 09:39 
[mesh manipulation] Importing Multiple Meshes  thomasnwalshiii  OpenFOAM Meshing & Mesh Conversion  18  December 19, 2015 18:57 
[ICEM] Generating Mesh for STL Car in Windtunnel Simulation  tommymoose  ANSYS Meshing & Geometry  48  April 15, 2013 04:24 
[Other] Mesh problem/ coarse OK  fine not OK  erichu  OpenFOAM Meshing & Mesh Conversion  10  April 10, 2013 12:29 