|
[Sponsors] |
September 12, 2016, 10:48 |
Kolmogorov microscale for a check valve
|
#1 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
Hello!
I try to estimate the Kolmogorov microscale for a hydraulic check valve. I attached a schematic of the valve. The Kolmogorov microscale is defined as where is the average rate of dissipation of turbulence kinetic energy per unit mass and is the kinematic viscosity of the fluid. It is the smallest scale in turbulent flow in which viscosity dominates and the turbulent kinetic energy is dissipated into heat. The Kolmogorov microscale can be expressed with the Reynoldsnumber Re as with the chracteristic length scale l. Now my questions: What is the characteristic length scale in may case of the check valve? Can I use the hydraulic diameter of the valve? I look in serval book, but I never found a satisfactory answer. The Re-number of the valve is with the flow rate Q, diameter of the valve and the viscosity . Thank you for an answer! Regards, Christoph |
|
September 12, 2016, 10:55 |
|
#2 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
what about the magnitude of the velocity inlet?
|
|
September 12, 2016, 11:05 |
|
#3 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
||
September 12, 2016, 11:24 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
So you have Re=O(10^3). It's quite small to try working in DNS/LES formulation. The kolmogorov scale is about 10^-4m
|
|
September 12, 2016, 13:00 |
|
#5 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
Thank you for the answer!
But 0.1 mm is quite large for the small valve: 0.5 mm plate stroke and the diameter of the inlet is 2 mm. What is the right characteristic length for this valve? I took the hydraulic diameter which is 2*h (plate stroke) = 1 mm. I performed different simulations in OpenFoam (DNS, LES one equation eddy visosity, and RANS k-omega SST) with the assumption of an axisymmetric flow problem. I know with turbulence it is not exactly true, but it was only a try. The differences in the flow rate between the three simulation runs were small. |
|
September 12, 2016, 13:09 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
Quote:
A rapid estimation is obtained by considering that at the Kolmogorov lenght scale you can set Re_eta=1 -> eta=ni/U. You cannot assume the axisymmetry if you want to use DNS/LES |
||
September 12, 2016, 13:50 |
|
#7 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
ok, thanks.
So, it is not impossible to perform a DNS simulation maybe from a quarter of the check valve with the corresponding boundary conditions. With this rough estimation arount 500 000 to 3 000 000 cells are necessary for a DNS simulation. |
|
September 12, 2016, 14:19 |
|
#8 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
you cannot simulate a quarter of valve
|
|
September 12, 2016, 14:47 |
|
#9 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
If I have time (and the workstation), I will try a DNS for the whole, half and quarter geometry and compare these simulations against each other. Because I am not sure, how large the error is.
Thank you again for your advice! |
|
September 12, 2016, 15:10 |
|
#10 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
Quote:
you can use that only in RANS. |
||
September 14, 2016, 04:10 |
|
#11 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
I know that from literature.
But I found a lot of papers dealing with axisymmetric direct numerical simulations for high Re-numbers. I did not find any comment in these papers, why it is possible. So, I tried it for my check valve example. Simulation is still running - but the first results show nearly no difference between DNS in 3D and 2D. |
|
September 14, 2016, 04:31 |
|
#12 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
Quote:
I am not aware of any relevant publication in the DNS/LES community of an axisymmetric DNS study. I suspect: 1) the 3D case has not yet developed the fully developed unsteady flow. 2) your code has a lot of numerical viscosity 3) you are using a too coarse grid |
||
September 14, 2016, 04:46 |
|
#13 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
Thank you for the very fast answers!
One of these papers is: http://www.sciencedirect.com/science...96890412002981 The pressure difference through the valve is up to 80 MPa, so in my opinion with the given dimension a high trubulent flow. I am using about 500 000 cells, but I will increase the number further. Maybe it could be the numerical viscosity, but how can I check it? regards, Christoph |
|
September 14, 2016, 05:14 |
|
#14 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
Quote:
I don't want to be offensive but such a paper would not be published in any good fluid mechanics journal... |
||
September 14, 2016, 05:25 |
|
#15 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
Do you have suggestion for good Journals?
|
|
September 14, 2016, 05:34 |
|
#16 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
||
September 14, 2016, 08:04 |
|
#17 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
||
September 14, 2016, 08:32 |
|
#18 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
Thank you!
With the results in this paper it is quite clear that it is not recommendable to use DNS in 2d. In my case the results are not so impressive and clear. Now, I have an other question: I performed two 2d simulation of the valve. Once with the k-omega-SST turbulence model and the second time without using a turbulence model. I found only minor differences between both simulations. But I do not know why the differences so small; maybe I resolve the larges eddies of the flow (but it should be not possible in 2d)? regards, Christoph |
|
September 14, 2016, 12:03 |
|
#19 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,773
Rep Power: 71 |
First, You cannot compare the RANS solution directly to the DNS solution. This latter must be statistically averaged. But if you do a 2D DNS there is no physical meaning in the solution you get, so the statistics are not relevant.
If you have numerical viscosity that depends on the type of discretization of the convective term. Generally upwinding or flux-limiters are the sources such error. |
|
September 15, 2016, 01:21 |
|
#20 |
New Member
Christoph Gradl
Join Date: Feb 2014
Posts: 18
Rep Power: 12 |
I performed some additional analysis of the results. And now I see clearly that the axissymmetric assumption for DNS simulations is not correct. In the attached figure I depicted the pressure iso-surface and the velcoity magnitude in the background.
You helped me a lot to understand turbulence simulation a little bit better! Unbenannt.PNG |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Flow through throttle valve in NX flow | Ashwath | Main CFD Forum | 0 | August 10, 2016 11:43 |
select valve seat for a valve in IC Engine | tarun_aseri | FLUENT | 1 | September 12, 2014 23:57 |
Modelling an initially closed valve which a fluid opens | siw | CFX | 1 | April 26, 2014 06:32 |
Gate valve flow simulations... | nikesh | FloEFD, FloWorks & FloTHERM | 5 | January 28, 2014 01:31 |
CFD of tilting disk check valve | Greg Anderson | Main CFD Forum | 1 | July 13, 1999 17:07 |